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A  METHOD  OF  ESTIMATING  THE  INFLUENCE  OF  SHIP'S  NOISE 
ON  AMBIENT  NOISE  MEASUREMENTS 

by 

Horst  Bendig 
Fried.  Krupp  GmbH 
Krupp  Atlas-Elektronik  Bremen 
D2800  Bremen,  West  Germany 


ABSTRACT 

Ambient  noise  measurements  are  often  adversely  affected  by  noise  from  pas¬ 
sing  ships.  Here,  a  method  is  presented  which  allows  the  estimation  of 
this  influence.  The  method  is  especially  suitable  for  shallow  water  appli¬ 
cations.  It  is  based  on  mode  interference.  The  broad-band  ship-generated 
noise  is  transmitted  to  the  measuring  sensor  by  modes.  Because  of  the  in¬ 
terference  of  these  modes,  the  spectrogram  of  the  noise  exhibits  inter¬ 
ference  patterns.  Disturbance  of  ambient  noise  measurements  occurs  when¬ 
ever  interference  patterns  appear  in  the  spectrogram.  The  degree  of  the 
disturbing  effect  is  specified  by  a  mode  interference  modulation  factor 
which  is  defined  here.  In  addition,  the  number  of  ships  and  their  distance 
from  the  measurement  sensor  can  be  estimated  from  the  interference  pat¬ 
terns.  Measurement  results  are  shown  which  have  been  evaluated  with  this 
method. 


1.  INTRODUCTION 

Ambient  noise  measurements  are  often  disturbed  by  noise  from  passing 
ships.  This  is  especially  the  case  in  sea  areas  which  carry  a  large  amount 
of  shipping,  such  as  the  North  Sea  and  the  Baltic. 

The  aim  of  this  investigation  is  to  Identify  the  disturbances  and  to  esti¬ 
mate  the  degree  of  disturbance.  Level  observations  alone  are  not  adequate 
for  solving  this  problem.  Ir  areas  mentioned  above,  which  are 

shallow-water  areas,  the  ambient  noise  level  itself  fluctuates  considerably 
U,  2],  so  that  it  is  difficult  to  distinguish  between  fluctuations  in 
ambient  noise  level  and  changes  in  disturbance  level.  For  this  reason,  the 
disturbance  is  to  be  characterised  in  a  different  way.  A  method  of  doing 
this  is  presented  here.  The  method  is  based  on  mode  propagation. 

The  structure  of  this  paper  is  as  follows:  In  Chapter  2,  some  basic  matters 
are  discussed.  In  Chapter  3,  the  mode  interference  modulation  factor  is  de¬ 
fined.  Chapter  4  deals  with  the  estimation  of  the  number  of  disturbance 
sources  and  their  distance  from  the  measurement  point*  In  Chapter  5, 
l  measurement  results  are  discussed.  Finally,  Chapter  6  gives  conclusions. 

i5:. 

!• 

I 

5,  I  * 

f:'  |  SACLANTCEN  CP-32 


13-1 


BENDIG:  Influence  of  ship's  noise 


2.  BASICS 

The  acoustic  energy  emitted  by  the  disturbing  ship  is  carried  to  the  mea¬ 
surement  point  with  the  aid  of  normal  modes.  This  is  a  well-known  fact 
(e.  g.  [3],  [43 •  [5]  ).  Because  of  the  mode  propagation,  the  spectrograms 
(defined  in  [6],  for  example)  of  the  measured  underwater  sound  contain 
interference  patterns  [7,  8]. 

The  spectrogram  is  the  power  P  versus  frequency  f  and  measurement  time  t. 
According  to  this,  the  function  for  the  spectrogram  can  be  written  as 
P(t,f).  In  cases  where  motion  is  involved,  the  time  may  be  replaced  by 
the  distance. 

An  example  of  a  measured  interference  pattern  is  shown  in  the  lower  part 
of  Fig.  1;  the  distance  is  plotted  against  the  frequency.  The  degree  of 
blackness  corresponds  to  the  acoustic  intensity.  The  associated  measure¬ 
ment  situation  is  sketched  in  the  upper  part  of  the  illustration.  If  the 
disturbing  source  of  noise  is  moving  away  from  the  measurement  point  at 
constant  speed,  the  modulation  of  the  power  in  each  narrow  analysis  band 
(e.g.  at  frequency  F  )  is  harmonic. 

Because  of  its  orderly  structure,  this  interference  pattern  is  easily  dis¬ 
tinguished  from  the  ambient  noise,  and  is  therefore  very  suitable  for  the 
detection  of  disturbance  sources.  This  is  substantiated  by  Fig.  2.  This 
illustration  shows  a  longer  spectrogram.  Here,  a  frequency  band  from  0  to 
450  Hz  is  shown.  The  ordinate  scale  is  not  in  distance  units  like  the 
scale  in  Fig.  1;  instead,  it  is  in  time  units  (the  usual  procedure  in 
spectrograms).  The  spectrogram  contains  an  interference  pattern  of  a  pas¬ 
sing  ship.  Shortly  after  the  starting  instant  t  =  0  of  the  spectrogram  re¬ 
cording,  the  measured  noise  had  not  yet  been  affected  by  the  ship  which 
would  later  cause  disturbance.  No  mode  interference  pattern  is  apparent 
yet  in  this  period  of  time.  The  pattern  does  not  begin  until  later;  it 
then  becomes  more  and  more  pronounced  as  the  distance  between  the  distur¬ 
bance  source  and  the  measurement  point  becomes  smaller.  After  the  closest 
£oint  of  approach  (CPA)  has  been  reached,  the  pattern  becomes  fainter  again. 

The  clear  distinction  between  spectrogram  regions  with  and  without  a  mode 
interference  pattern,  as  illustrated  in  Fig.  2,  causes  us  to  make  the 
following  statement:  whenever  interference  patterns  are  present  in  the 
spectrogram,  it  should  always  be  assumed  that  a  source  of  disturbance  is 
present.  The  degree  of  disturbance  is  to  be  described  by  a  modulation  fac¬ 
tor,  which  is  referred  to  here  as  the  "mode  interference  modulation  factor". 

In  addition  to  the  degree  of  acoustic  disturbance  of  the  ambient  noise 
measurement,  a  rough  estimate  of  the  number  of  disturbance  sources  and  their 
distance  from  the  measurement  point  is  often  of  interest.  Mode  interference 
offers  possibilities  for  this.  A  glance  at  the  spectrogram  in  Fig.  2  will 
verify  this.  The  interference  lines  which  charaterise  the  fan-like  struc¬ 
ture  of  the  interference  pattern  have  different  gradients.  At  a  given  fre¬ 
quency,  the  absolute  value  of  the  gradient  of  the  interference  lines  in¬ 
creases  with  the  distance  between  the  source  and  the  receiver.  This  fact 
forms  the  basis  of  the  distance  estimate. 
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Let  us  mention  the  following  boundary  conditions  for  the  method  proposed 
here  for  estimating  the  influence  of  disturbance  sources  on  the  ambient 
noise  measurement. 

The  sound  propagation  channel  must  be  such  as  to  ensure  that  mode  propa¬ 
gation  can  take  place.  This  is  so  in  shallow  water  areas  especially.  The 
disturbance  source  should  radiate  acoustic  energy  over  a  wide  band  of  fre¬ 
quencies.  This  condition  is  fulfilled  for  the  broad-band  cavitation  and 
flow  noise  produced  by  moving  ships.  Here,  the  frequency  band  from  a  few 
tens  of  Hz  to  about  500  Hz  is  investigated.  The  underwater  sound  is  to  be 
measured  omnidirectionally. 


3 *  DETERMINATION  OF  THE  MODE  INTERFERENCE  MODULATION  FACTOR 

The  mode  interference  modulation  factor  is  intended  to  describe  the  degree 
of  modulation  of  a  periodogram  component.  It  is  thus  a  measure  of  the  mea¬ 
sured  time-dependent  variation  of  the  power  in  the  narrow  analysis  bands. 

To  be  able  to  compute  the  degree  of  modulation  formally,  we  make  the  fol¬ 
lowing  assumptions.  We  assume  that  each  periodogram  component  can  be  des¬ 
cribed  approximately  by  a  stochastic  process  Z(t): 

Z  (t)"XU)-AU)  (3.1) 

Equation  (3.1)  is  valid  for  a  fixed  frequency.  The  frequency  is  suppressed 
in  the  notation.  According  to  the  function  for  the  power  P ( t , f )  in  the 
spectrogram,  our  function  Z(t)  approximates  the  power  P(t,  f  =  const.)  at 
constant  frequency. 

The  product  process  Z ( t }  is  to  consist  of  a  carrier  process  X ( t )  and  a  mo¬ 
dulation  process  A(t),  which  are  both  assumed  to  be  stationary. 

We  shall  assume  that  the  carrier  process  and  the  modulation  process  are  in¬ 
dependent  of  each  other.  We  shall  also  assume  that  short-time  power  spectra 
computed  one  after  another  are  likewise  independent  of  each  other. 

In  the  harmonic  modulation,  the  modulation  factor  is_defined  as  the  quotient 
of  the  amplitude  and  mean  value;  the  amplitude  isV2~times  the  standard 
deviation  Or  . 

We  shall  adopt  this  definition  and  describe  the  mode  interference  modulation 
factor  M  of  a  periodogram  component  by: 

0.2) 


whore  <3*  is  the  standard  deviation  of  the  modulation  process  and  ^  is  its 
mean  value. 
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Since  the  product  process  only  (and  not  the  modulation  process)  is  measur¬ 
able,  the  quotient  fyl/Uj  cannot  be  determined  directly.  For  this  reason, 
the  modulation  process  will  be  described  by  the  obtainable  relative  stan¬ 
dard  deviations  of  the  product  process  Z(t)  and  of  the  carrier  process 
X(t).  We  obtain  the  following: 


(3.3) 


We  now  go  over  to  the  discrete  case,  since  discrete  periodograms  are  com¬ 
puted  in  practice,  e.  g.  with  the  aid  of  the  fast  Fourier  Transform.  The 
modulation  factor  is  thus  computed  for  discrete  frequency  values  f.  only. 
We  also  assume  that,  for  estimating  the  modulation  factor,  only  time-in¬ 
tervals  of  duration  T  are  used  in  the  spectrogram.  The  position  of  an  in¬ 
terval  with  respect  to  time  is  assigned  to  the  discrete  instant  t.-.  The 
modulation  factor  will  therefore  be  symbolised  by  My(t.,f.j)  .  Averaging 
in  the  frequency  direction  over  K  modulation  factors  ’'of  adjacent  spec¬ 
trogram-components  gives  the  mean  mode  interference  modulation  factor  for 
the  instant  tj: 


_  ,  ,  UK-1 

MrW  “  IT  S.  ,  (3.4) 

i  »  L 


where  fj  «*  i  •  Af  and  Af  is  the  frequency  resolution  cell  of  the  discrete 
periodogram.  During  the  ambient  noise  measurement,  this  mean  value  Rr  is 
continually  updated.  Depending  on  the  magnitude  of  RT,  a  decision  is  then 
made  as  to  whether  the  disturbance  is  too  great  or  is  still  acceptable. 

To  estimate  the  modulation  factor  in  accordance  with  equation  (3.3),  we 
need  not  only  the  relative  standard  deviation  (which  can  be 

estimated  directly  from  the  spectrogram)  but  also  the  relative  standard 
deviation  for  the  carrier  process.  (In  the  following,  estimated 

values  are  indicated  by  a  "hat14,  *  ,)This  quotient  Is  generally  not  known. 
The  required  quotient  6"*  //<*  is  a  constant  and  can  be  computed  after 
a  special  procedure  has  been  applied.  This  procedure  is  called  “normalisa¬ 
tion4*  [9].  This  normalisation  is  applied  to  each  individual  periodogram 
of  the  spectrogram  before  determination  of  the  modulation  factor. 

To  summarise,  the  expected  value  of  the  power  density  in  the  neighbourhood 
of  the  individual  frequency  points  f .  of  the  original  periodogram  is  esti¬ 
mated  first.  Next,  the  values  of  Unoriginal  periodogram  are  divided  by 
the  estimated  expected  values.  An  effect  of  this  normalisation  is,  for 
example,  that  broad-band  pink  noise  turns  into  white  noise  (more  or  less!. 
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Furthermore  -  and  this  is  important  for  us  -  in  the  unmodulated  case  all 
moments  for  the  components  of  the  normalised  peridogram  can  be  stated,  and 
thus  also  the  required  values^  and  of  our  carrier  process. 

We  expect  the  mean  modulation  factor  My  for  an  unmodulated  spectrogram  to 
be  around  zero.  In  this  case  especially,  it  is  clear  that  estimated  rela¬ 
tive  standard  deviations  of  individual  periodogram-components  may 

be  smaller  than  the  calculated  quotient  I Mk  of  the  carrier  process. 

The  modulation  factor  according  to  equation  (3.3)  would  become  imaginary. 
However,  this  result  is  not  meaningful.  For  this  reason,  equation  (3.3)  is 
modified  somewhat: 


(3.5) 


in  equation  (3.5),  the  radical  expression  is  now  unsigned.  The  sign  has  been 
moved  to  a  place  in  front  of  the  root,  and  is  determined  by  the  operation 
"sign{.j  The  following  applies: 


sign  (  c}  ■ 


1,  c  >  0 
0,  c  -  0 
.-1,  c  <  0 


(3.6) 


According  to  equation  (3.5),  the  estimate  may  yield  negative  modulation  fac¬ 
tors.  This  is  unimportant  as  only  the  mean  value  of  the  modulation  factor 
must  be  zero  or  positive. 


4 •  ESTIMATION  or  THE  NUMBER  OF  SHIPS  ANQ  THUB  DISTANCE 

By  means  of  the  interference  patterns  in  the  spectrogram,  it  is  possible  to 
estimate  the  number  of  disturbing  ships  and  to  make  a  rough  estimate  of  the 
distance.  However,  this  requires  the  spectrogram  to  be  available  in  a  form 
which  will  allow  visual  evaluation,  e.  g.  as  a  display  on  a  grey-shade  recor¬ 
der  or  on  a  TV  screen. 

If  the  ambient  noise  measurement  is  being  disturbed  by  more  than  one  ship, 
the  corresponding  interference  patterns  are  superimposed  on  one  another.  In 
the  spectrogram,  this  super imposition  is  clearly  recognisable,  especially 
for  the  practised  eye.  An  example  is  shown  in  Fig.  3.  There,  two  principal 
sources  of  disturbance  are  involved.  The  superimposition  is  clearly  visible 
in  the  time  interval  from  15:10  to  15:17. 

The  interference  pattern  in  Fig.  2  had  shown  that  the  pattern  depends  on  dis¬ 
tance.  The  magnitude  of  the  gradient  of  the  "interference  lines"  is  propor- 
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tional  to  the  distance  between  the  disturbance  source  and  the  measurement 
point.  This  fact  is  used  to  make  a  passive  estimate  of  distance.  In  addi¬ 
tion  to  the  modulation  factor,  we  thus  also  obtain  an  important  parameter 
for  situation  assessment  in  ambient  noise  measurement. 

The  method  described  here  works  with  relatively  simple  aids.  It  requires 
that  the  sources  of  disturbance  are  approximately  on  a  collision  course 
with  respect  to  the  measurement  point.  Since  this  requirement  is  only  more 
or  less  fulfilled  in  practice,  the  estimate  can  only  produce  rough  values 
for  the  distance. 

Fig.  4  summarises  the  main  steps  in  the  evaluation.  In  addition  to  the 
measured  interference  pattern,  results  from  the  model  calculation  are  re¬ 
quired:  for  the  evaluation  frequency  Fg,  we  require  the  interference  wave¬ 
length  IWl  (also  called  the  "horizontal  interaction  distance")  A  of  the 
mode  pair  which  is  dominant  in  the  measured  interference  pattern.  Second¬ 
ly,  the  frequency-dependent  changes  of  the  IWL,  3A/3f  ,  must  be  deter¬ 
mined.  These  tasks  can  be  performed  by  computer  programs  such  as  those 
described  in  [10]  or  [11]  or  by  much  simpler  procedures  [12]. 

Two  things  must  be  determined  from  the  measured  spectrogram  -  firstly,  the 
number  k  of  interference  wavelengths  passed  through  per  unit  time  by  the 
disturbing  ship  (k  does  not  have  to  be  an  integer),  and  secondly  the  gra¬ 
dient  of  the  interference  lines. 

The  expression  k  *A/T  represents  the  speed  of  approach  or  departure,  v. 
The  distance  can  be  estimated  from  the  determined  values  according  to  the 
following  formula: 


(4.1) 


In  this  formula,  the  interference  line  gradient  H  /8f  is  the  most  im¬ 
portant  parameter.  The  other  parameters  act  as  scaling  factors  to  some 
extent. 


It  is  immediately  apparent  that  k  can  be  estimated  with  sufficient  reliabi¬ 
lity  if  and  only  if  the  interference  maxima  along  the  line  f  *  P0  are  equi¬ 
distant.  However,  this  is  so  only  if  the  speed  of  approach  or  departure  of 
the  disturbance  source  (relative  to  the  measurement  point)  is  constant  with¬ 
in  the  evaluation  time  T. 


5.  DISCUSSION  OF  THE  MEASUREMENT  RESULTS 

The  evaluation  results  of  many  experiments  at  sea  show  that  the  proposed 
method  for  estimating  the  influence  of  disturbance  noise  is  practicable. 


******  rm 
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The  measured  data  come  from  North  Sea  areas  with  a  water  depth  of  about 
43  m.  The  measurements  were  performed  in  the  spring  of  two  successive 
years.  An  omnidirectional  hydrophone  was  used  to  receive  the  sound. 

The  mean  mode  interference  modulation  factor  Ry  was  evaluated  in  the  fre¬ 
quency  band  from  50  Hz  to  450  Hz.  The  duration  T  of  the  evaluation  in¬ 
terval  was  fixed  at  128  sec.  The  normalisation  produced  a  computed  value 
of  1.054  for  Gxl/Uf. 

All  evaluations  gave  similar  results,  so  we  can  restrict  ourselves  here 
to  a  discussion  of  a  few  examples.  In  Figs.  5,  6  and  7,  the  modulation 
factor  Ry  is  plotted  against  time.  The  shaded  area  symmetrical  to  the 
modulation  factor  Ry  indicates  the  standard  deviation  of  the  K  =  401 
individual  modulation- factors.  A  new  estimation  of  the  modulation  fac¬ 
tor  was  carried  out  every  20  sec,  so  there  is  an  overlap  of  108  periodo- 
grams. 

The  modulation  factor  My  in  Fig.  5  shows  the  typical  variation  with  time 
as  a  ship  passes  by:  at  the  beginning  of  the  experiment,  the  disturbing 
ship  is  far  away  and  the  modulation  factor  is  small.  As  the  distance 
between  the  ship  and  the  measurement  point  decreases,  the  modulation 
increases.  After  the  ship  has  passed  the  closest  point  of  approach  (CPA), 
the  modulation  decreases  again.  The  drop  in  the  modulation  factor  in 
the  region  of  the  CPA  is  due  to  the  particular  form  of  the  interference 
pattern  in  this  region.  A  glance  at  the  spectrogram  in  Fig.  2  which  re¬ 
fers  to  this  experiment  will  make  this  clear. 

At  ‘N*  beginning  and  end  of  the  experiment,  the  modulation  factor  in 
Fig.  5  stays  at  a  value  of  about  0.2,  However,  the  passing  ship  did  not 
have  any  disturbing  influence  here.  Reduction  of  the  modulation  factor 
to  zero  was  prevented  here  by  slight  disturbances  of  a  non-definable 
type.  These  disturbances  can  be  recognised  in  the  original  plot  of  the 
spectrogram.  (Owing  to  the  copying  process,  the  grey  shades  are  unfortu¬ 
nately  not  reproduced  very  well,  so  that  it  is  no  longer  possible  to  see 
all  of  the  details  in  the  reproduced  spectrograms. ) 

If  sources  of  disturbance  are  present,  the  model  is  an  appropriate  one. 
This  has  been  shown  by  the  measurement  results.  It  has  not  yet  been  pos¬ 
sible  to  carry  out  the  test  for  the  case  in  which  disturbance  sources 
are  absent,  since  there  were  always  slight  disturbances  in  the  sea  areas 
at  all  available  measurement  times. 

The  noise  for  Fig.  $  had  been  measured  with  a  hydrophone  situated  about 
1.5  m  above  the  bottom.  Fig,  6  shows  an  example  from  a  measurement  5  m 
below  the  surface  of  the  water.  Here  too,  a  ship  passed  by.  The  princi¬ 
pal  difference  between  the  two  examples  is  to  be  found  in  the  time  period 
after  the  instant  of  CPA:  the  decrease  in  the  modulation  factor  with 
time  is  much  steeper  in  Fig.  6. 
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Whereas  in  Fig.  5  the  increase  and  decrease  of  the  modulation  factor  be¬ 
fore  and  after  the  CPA  are  roughly  symmetrical,  the  variation  in  Fig.  6 
is  very  asymmetrical.  This  effect  was  found  generally.  It  is  probably 
due  to  a  strong  screening  effect  produced  by  the  wakes  of  the  ships. 

Fig.  7  shows  the  modulation  factor  My  as  two  sources  of  disturbance  pass 
by.  In  the  first  25  minutes  of  the  experiment,  the  first  disturbance  source 
is  dominant,  and  then  the  second  one  is  dominant.  A  corresponding  excerpt 
of  the  spectrogram  for  the  time  period  from  10  minutes  to  15  minutes  af¬ 
ter  the  start  of  the  experiment  is  shown  in  Fig.  3. 


No  representative  results  can  be  offered  for  the  distance  estimate  since 
the  accuracy  of  the  results  depends  to  a  very  large  extent  on  the  amount 
of  practice  which  the  estimator  has  had  -  especially  in  the  case  of  com¬ 
plicated  interference  patterns.  We  want  to  show,  by  means  of  an  example, 
that  the  evaluation  method  works  successfully. 


We  shall  evaluate  a  part  of  the  experiment  which  gave  the  spectra  shown 
in  Fig.  2.  For  this  experiment,  the  model  calculation  gives  an  interfe¬ 
rence  wavelength  780  m  and  a  gradient  PA  &  3.55  m  /  Hz, 

determined  at  an  evaluation  frequency  F0  =  200  Hz.  In  the  spectrogram, 
the  number  of  interference  wavelengths  which  have  been  passed  thro1  gh  in 
5  minutes  is  estimated  as  3.9.  Because  the  pattern  is  very  complicated, 
a  "ten  point  divider"  is  required  for  this.  With  the  aid  of  this  instru¬ 
ment,  the  equidistant  interference-maxima  can  be  found  quickly.  The  dis¬ 
tances  estimated  after  the  dt/df  evaluation  are  listed  for  various 
times  in  Table  I . 


TIME  SINCE 
START  OF 
EXPERIMENT 

(minutes) 


ESTIMATED 

DISTANCE 

(nautical 

miles) 


DISTANCE 
MEASURED 
BY  RADAR 
(nautical 
miles) 


9 

3.4 

3.55 

10 

3.0 

3.25 

11 

2.6 

2.95 

12 

2.4 

2.55 

13 

2.1 

2.25 

14 

1.8 

1.8 

15 

1.5 

1.5 

TABLE  I:  EXAMPLE  OF  DISTANCE  ESTIMATION 

For  the  sake  of  comparison,  the  radar  measurements  are  also  listed  here. 
They  are  not  normally  available  during  application  of  the  proposed  method. 
Their  purpose  here  is  to  act  as  reference  values.  The  comparison  shows 
that  the  estimate  gives  acceptable  values. 
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6.  CONCLUSIONS 

A  method  was  to  be  developed  with  which  the  influence  of  disturbing  ship's 
noise  on  the  ambient  noise  measurement  could  be  estimated.  This  objective 
has  been  achieved  by  the  defining  of  a  mode  interference  modulation  factor. 
The  theoretical  model  selected  for  this  can  be  used  in  practice.  This  has 
been  demonstrated  with  the  aid  of  measurement  results. 

The  additional  information  gained  from  mode  interference,  namely  informa¬ 
tion  about  the  distance  of  the  disturbance  sources,  extends  the  basis  for 
situation  assessment  during  ambient  noise  measurement. 

The  mode  interference  modulation  factor  offers  the  possibility  of  automa¬ 
ting  the  monitoring  of  the  measurement  process;  the  idea  is  to  use  a  thres¬ 
hold  decision. 

In  the  evaluation  of  the  mode  interference  modulation  factor  carried  out  so 
far,  all  periodogram  components  within  a  specified  frequency  band  have  been 
taken  into  account.  However,  it  may  be  better  to  evaluate  only  some  of  the 
periodogram  components  within  a  specified  frequency  band,  namely  those  with 
the  largest  modulation  factors.  In  this  way,  frequency  bands  involving  poor 
propagation  behaviour  would  be  excluded  from  the  evaluation  process.  This 
would  lead  to  an  improvement  in  the  estimation  result. 
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OPTIMAL  DETECTION  AND  TRACKING  OF  ACOUSTICAL  NOISE  SOURCES 
IN  A  TIME-VARYING  ENVIRONMENT 

by 

Hendrik  Van  Asselt 
SACLANT  ASW  Research  Centre 
La  Spezia,  Italy 


ABSTRACT 


A  dynamical  detector  estimator  is  derived  for  nonlinear  nongaussian  signal 
processes  based  on  Bucy's  representation  theorem  which  is  the  fundamental 
result  in  stochastic  filtering.  A  quantized  version  of  two-dimensional 
conditional  probability  density  functions  is  generated  for  bearing  and 
bearing- rate  as  well  as  for  frequency  and  frequency-rate.  Results  of  the 
analysis  of  real  data  for  bearing  and  bearing  rate  are  presented  and 
discussed  using  optimal  nonlinear  filtering  techniques. 


INTRODUCTION 


Target  detection  and  tracking  from  passive  sensors  is  a  nonlinear 
estimation  problem  and  has  practical  significance.  This  paper  discusses 
this  problem  and  presents  an  analysis  of  real  target  data  measured  in 
shallow  water. 

The  present  state  of  the  art  makes  it  feasible  to  develop  detection  and 
estimation  algorithms  in  discrete  time,  since  nonlinear  processing 
algorithms  based  on  Bayes  law  are  developed  in  discrete  time. 

In  this  paper  we  have  restricted  ourselves  to  hardware  that  is  available 
today,  which  means  that  the  closeness  of  approximations  to  optimality 
depends  largely  on  the  dimensionality  of  the  signal  processes. 

We  have  divided  our  system  into  two  parts,  bearing  and  frequency 
estimators,  which  each  turn  out  to  be  two-dimensional  dynamical  processes 
having  nonlinear  system  functions  and  linear  observation  functions. 
Because  of  the  lower  dimensionality  of  bearing  and  frequency,  when  the 
state  spaces  are  decoupled,  one  can  use  a  more  accurate  approximation  to 
optimality  than  for  the  coupled  state/space  representation  if  one  thinks  In 
terms  of  the  number  of  computations  required  and  of  the  real-time 
implementation. 

For  bearing  and  frequency  estimation  we  are  able  to  generate  quantized 
versions  of  two-dimensional  conditional  probability  density  functions  and 
represent  these  as  a  set  of  point  masses  on  a  moving  mesh.  A  detailed 
discussion  Is  given  in  <1>. 


SACLANTCEN  CP- 32 


VAN  ASSELT:  Detection  and  tracking  of  noise  sources 


Our  motivation  for  this  type  of  data-processing  scheme  is  that  the  nature 
of  most  signals  relating  to  passive  sonar  detection  and  tracking  are 
multidimensional  signal  processes  and  are  commonly  dynamical,  nonlinear, 
and  nongaussian  and,  further,  may  have  time-varying  statistics. 
"Classical"  signal-processing  schemes  based  on  the  stationary  steady-state 
filtering  theory  therefore  often  give  undesired  performance,  especially  in 
cases  of  low  signal-to-noise  ratios. 


1  SYSTEMS  CONSIDERATIONS 

Kalman  filtering  is  often  applied  for  designing  suboptimal  nonlinear 
estimators.  The  procedure  is  that  the  models  are  linearized  and  then  the 
linear  theory  is  applied.  This  is  one  of  the  simplest  methods  of  deriving 
suboptimal  data-processing  schemes  for  nonlinear  dynamical  signal 
processes.  Often,  however,  these  techniques  lead  to  undesirable 
performance,  as  we  will  illustrate  in  a  simple  example. 

Figure  1  compares  a  recursive  realization  of  the  kalman  filter  with  the 
discrete  Bayes-Law  estimator  for  the  case  of  truncated  gaussian  densities. 
The  a-priori  conditional  probability  density  function  is  indicated  by  pn_^ 

having  a  variance  P  ^  =  0.6.  The  probability  density  function  of  the 

observations  Z  has  a  variance  R  =1.2  and  the  system  noise  has  a 
n  n 

probability  density  function  m  and  variance  Qn  =  0.696. 

If  we  calculate  the  a  posteriori  conditional  probability  density  function, 
then  both  the  kalman  estimator  and  the  Bayes-Law  estimator  provide  almost 
the  same  output  variance:  Pg  =  0.6  and  P^  =  0.62.  The  difference  results 

from  the  fact  that  Bayes  Law  computes  on  the  real  truncated  gaussian 
densities  whereas  the  kalman  filter  assumes  non-truncated  densities. 

Figure  2  makes  the  same  comparison.  Consider  now  the  real  a-priori 

conditional  probability  density  function  is  nongaussian,  whereas  the  kalman 
filter  assumes  It  is  gaussian,  both  having  the  same  variance.  Also  the 
real  probability  density  function  Z  of  the  input  data  is  nongaussian  and 
again  the  kalman  filter  assumes  it  is  gaussian. 

Comparing  the  results  of  this  example  we  note  that  a  kalman  filter  severely 
underestimates  the  conditional  a-posteriori  variance.  As  a  consequence  it 
will  also  severely  underestimate  the  filter  gain  matrix.  This  causes  too 
little  information  to  be  extracted  from  the  measurements ,  leading 
frequently  to  divergence.  This,  of  course,  is  not  to  blame  the  kalman 
filter,  It  is  a  result  of  its  wrong  application. 


1*1  Implementation  of  a  Bayes-Law  System 

Figure  3  shows  the  configuration  of  a  Bayes-Law  system.  The  input  is 
obtained  from  a  line  array  of  20  hydrophones  approximately  2  m  apart.  This 
is  followed  by  a  conventional  beaaform er  in  order  to  sample  the  noise  field 
passively.  The  beam-axes  are  1°  apart  and  the  beaawidth  is  5°  at  the  3  dB 
level. 
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All  the  beams  do  not  necessarily  carry  the  same  levels  of  energy  and  the 
energy  levels  per  beam  fluctuate  within  a  range  of  approximately  10  dB. 

Along  the  array  we  observe  the  signal,  which  is  "more  or  less"  coherent. 
We  cannot  specify  a  measure  of  coherence,  because  it  is  not  a  constant. 
After  beamforming  we  could  observe  this  as  presented  in  Fig.  4.  The  shaded 
zone  shows  the  collection  of  2000  beamformer  output  samples,  where  we  have 
eliminated  time  dependence.  If  this  shadowed  zone  were  to  collapse  to  a 
broken  line  we  would  receive  a  fully  coherent  signal;  so  a  measure  of 
noncoherence  is  indicated  by  the  width  of  the  shaded  area.  Note  also  that 
the  noise  field  is  nonhomogeneous. 

Given  these  sample  functions  we  then  evaluate  the  likelihood  function  in 
order  to  generate  a  probabilistic  structure  in  space. 

In  Fig.  3  these  observations  are  indicated  as  Z[h(£),  G(£)],  where  h(£)  is 
the  input  function  and  G(£)  the  array  gain  in  a  given  direction.  We  then 
continue  by  two-dimensional ly  convolving  the  probabilistic  structure 
described  above  with  the  product  of  the  probability  density  function 
ro(y,  f(£).rl  describing  the  target  trajectory  noise  and  the  a-priori 

conditional  probability  density  function  D  ,((,)  (see  also  Ch.  4  of  <2>, 
especially  pp  43-59).  1 

We  have  m  modelled  as  a  gaussian  density  function,  the  arguments  y  (and 
dummy  £)  describing  the  target-motion  variables  in  state/space  (bearing, 
bearing  rate  and  frequency,  frequency  rate)  and  the  target-motion  model 
represented  by  f(£),  where  r  Is  the  average  manoeuvring  time.  During  the 
detection  phase  the  target  motion  is  modelled  as  a  linear  trajectory  in 
cartesian  coordinates  having  random  accelerations  (the  constant-speed 
model).  During  the  tracking  phase  we  use  a  constant-acceleration  model  in 
order  to  allow  target  manoeuvring.  This  can  be  achieved  by  changing  the 
autocorrelation  time. 

After  we  have  generated  the  conditional  probability  density  function,  it  is 
simple  to  compute  mean  values  and  variances,  as  is  indicated  in  Fig.  3. 
Figure  5  presents  the  estimated  variances  of  a)  bearing  and  d)  bearing- 
rate  and  shows  the  potential  of  the  nonlinear  optimal  estimator. 

Figure  6  shows  a  typical  example,  where  the  dynamics  behave  reasonably 
stationary.  We  cut  the  density  function  into  two  halves  along  the 
correlation  axis  In  ordt‘r  to  study  its  transient  during  the  detection  phase 
(Fig.  4).  At  the  start,  t  ,  we  have  no  information,  which  is  represented 

by  a  uniform  initial  probability  density  0< t0>.  Then  at  intervals 

of  2  seconds  we  update  the  estimated  density  D(  ),  e.g.  tA  refers  to 
10  seconds  processing  time.  We  stop  the  detection  phase  after  t*.  which  is 
where  the  detector  has  reached  the  steady-state  solution.  Table  1  presents 
the  detector  performance  numerically,  with  reference  to  the  highest,  level, 
which  is  at  91°  and  zero  degrees  per  second.  After  the  detection  phase  the 
autocorrelation  time  of  the  target  accelerations  is  Increased  so  that  the 
transient  is  observed  from  detection  to  tracking  phase.  Figure  7  shows  a 
sequence  of  pictures  showing  both  transients.  Figure  8  is  another  example, 
In  which  the  estimators  have  been  left  in  the  detection  mode  so  as  to  show 
the  detectors  capability  of  handling  dynamical  processes  (a  relative  fast 
ship  moving  through  a  fixed  set  of  beams). 


H 


SACLANTCEN  CP-32 


14-3 


VAN  ASSELT:  Detection  and  tracking  of  noise  sources 


A  sequence  of  2000  sequentially  estimated  conditional  probability  densities 
has  been  filmed  to  show  in  more  detail  the  systems  dynamics.  This 
pictorial  display  demonstrates  why  techniques  based  on  stationary 
steady-state  filtering  theory  cannot  work  properly.  This  is  because  often 
averaging  time  is  too  long  with  respect  to  the  non-stationarity  of  the 
underlying  signal-process  <3>. 


2  DETECTOR  DERIVATION 

The  stochastical  differential  equation  representing  a  diffusion  process  is: 

dx^  =  f(xt>  t)dt  +  a(x^,  t)du^  ,  (Eq.  1) 

where  xt,  u^  and  f(xt>  t)  are  n-dimensional  vectors 
and  t)  is  an  n  x  n  matrix. 


Assume  that  a  solution  exists  and  is  unique  and  that  the  domain  of  the 
solution  will  be  the  bounded  time  interval  (a,  b).  The  process 

T 

{u^}  ={(u^  U£  t . un  t^  ^  defines  the  vector  of  n  independent  Brownian 

motions.  The  components  of  f (x^, t)  and  o(xt,t)  are  continuous  In  t  and 
lipschitz-continuous  in  xt- 

For  the  m-dimensional  observation  process  we  use  the  differential  equations 


dzt  »  h(xt,  t)dt  +  y(t)dvt  signal  present 

dzt  a  v(t)dvt  signal  absent 

whore  dxt  is  defined  in  Eq.  1, 

x(a)  =  c, 

i( a)  “  0, 

vt  is  an  m-dimensional  Brownian  motion,  independent  of  u^. 


(Eq.  2) 


(Eq.  3) 


Then,  following  <4>,  consider 

dxt  a  f(xt,  t)dt  ♦  o(xt,  t)dut  (Eq.  4) 

dzt  a  (f(xt,  t)  ♦  a(iv  t)  y(zt,  t))dt  +  o(zt,  t)dut  ,  (Eq.  5) 

where  t<g{a,  b]  x(a)  =  z(a) 

f,  y  are  n-vectors  and  o  is  an  n  x  n  matrix 
f,  o  and  y  are  measurable  in  both  arguments 

t  2 

/ 1  y(xt»t)  |  dt  <  *  w.p.  i 

Jy(x.,t)  |  <  yjlx.  I)  ;  yn  being  a  nondecreasing  function  of  a 
I  realWiatle?  *  0 

1 

i 
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Then  the  measures  m  and  m  induced  on  the  space  of  all  continuous 

functions  in  R  ,C  [a,bj  by  the  processes  {x.}  and  {z.}  are  mutually 
absolute  continuous?  The  Radon-Nikodym  derivative  dm^/dmz  is  then  given  by 


(Eq.  6) 

The  nonlinear  detection  and  estimation  problem  is  described  <2>  by: 

dx^  =  f(x^,t)dt  +  o(xt,t)du^  (Eq.  7) 

dz^  =  h(t)xtdt  +  y(t)dvt  (Eq.  8) 


dm. 


=  exp 


/  YT(xs,s)dvs  -  |  J*  |y(xs,s)|2  ds 


The  conditional  mean  for  this  estimation  problem  is  given  by 


* 

El*t  I  zr;a  <  »'  <  t]  =  xt  = 


(Eq.  9) 


The  expectation  symbol  E  defines  integration  with  respect  to  the  function 
space  measure  mx  generated*  by  solving  Eq.  7  and 


«t  s  flt^xr»zPi  a  <  t  <  r)  * 
exp 


/  xphplvjvrl  *d2r  *  H  /  ^hJtY^r'hpXj.dr 

A  B 


Recalling  the  detection  problem  of  Eq.  2  and  defining  the  process  {P^J  as 
Pt  B  **•  E«JQt*  * 

which  satisfies  the  stochastic  differential  equation 

dPt  *  »Jh[hJytr1  dzt  -  ^  (Eq.  10) 

where  is  given  by  Eq.  9. 

This  can  be  proved  because  IL-ftl*]  is  positive  and  finite  w.p.  1;  therefore 
Pt  is  well  defined.  If  w2r  apply  the  stochastic  differential  rule  to 
**  ob*a*n 

w 


,  EwJKtht°tl  *YtYt*  E«x^htxtflt^ 
* - - 
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and  we  find  the  likelihood  function 

dpt  =  xtht(Vt1_1  dzt  + 

■  **  ht  *t dt 

This  concludes  our  detection  process  for  a  diffusion  process  in  white 
noise. 


CONCLUSIONS 

A  dynamical  detector/estimator  is  presented  for  nonlinear  nongaussian 
signal  processers.  We  have  shown  the  performance  of  the  bearing/bearing 
rate  estimator  where  the  input  was  real  data. 
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TRUNCATED  GAUSSIAN  EXAMPLE 
FOR  LINEAR  SCALAR  SYSTEM 


«*n-l 


SHIFT  m  UNTIL  k.  5 
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—  REAL  pdt  — GAUSS  P„  ,  »  0.6 

I  I 

OBSERVATION  DENSITIES 

—  REAL  pdl  —GAUSS  R„*  1.2 


SYSTEM  NOISE  DENSITIES 
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SYSTEM  CONFIGURATION 
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QUANTIZED  CONDITIONAL  PROBABILITY  DENSITY  FUNCTION 

NUMERICAL  VALUES 


-106 

-131 

-134 

-147 

-164 

-194 

-203 

-220 

-302 

-  co 

-  00 

2.5 

-  81 

-  98 

-  70 

-106 

-127 

-130 

-138 

-179 

-187 

-281 

**  00 

2 

-  75 

-  73 

-  69 

-  65 

-  65 

-  97 

-115 

BH 

-122 

-171 

-281 

1.5 

-  74 

-  48 

-  62 

-  37 

-  61 

-  33 

-  80 

-  96 

-106 

-151 

1 

-  66 

-  58 

-  51 

-  42 

-  33 

-  29 

-  33 

-  66 

-  77 

-  84 

-  85 

0.5 

-  94 

-  89 

-  78 

-  73 

-  62 

... 

-  32 

n 

-  47 

-  65 

-  63 

-  72 

0 

-107 

-  98 

-  90 

-  83 

• 

-  74 

-  63 

-  33 

-  29 

-  36 

-  45 

-  53 

-0.5 

-172 

-137 

-128 

-119 

-117 

-  94 

-  65 

-  32 

-  50 

-  39 

-  55 

-1 

-302 

-201 

-151 

-136 

-135 

-122 

-112 

-  94 

-  65 

-  67 

-  71 

-1.5 

-  CO 

-303 

-216 

-201 

-  84 

-165 

-157 

-145 

-121 

-100 

■  72 

•2 

33 

-  * 

-333 

-240 

-234 

-225 

-194 

-175 

-162 

-150 

-134 

-2.5 

85  | 
_ L 

86 

87 

68 

89 

90 

91 

92 

93 

94 

95 

V.  • 
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-  SPATIAL  BINS  (degrees) 


VALUES  IN  dB 
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INFLUENCE  OF  BACKGROUND  NOISE  SPATIAL  COHERENCE 
ON  HIGH  RESOLUTION  PASSIVE  METHOD 

by 

Georges  BIENVENU  and  Laurent  KOPP 

THOMSON-CSF,  ASM  Division,  BP  53,  06801  CAGNES-sur-MER  Cedex,  FRANCE 


ABSTRACT 


New  signal  processing  methods  for  passive  listening  have  appeared  recently  : 
they j  are  called  "high  resolution"  methods  because  they  have  theoretically  a 
better  resolving  power  than  conventional  or  adaptive  beamforming.  But  they 
need  additional  knowledge  on  background  noise  structure,  and  particularly 
on  ambient  noise  for  low  frequency.  These  methods  are  based  on  the  estima¬ 
tion  of  the  spectral  density  matrix  of  the  signal  received  on  the  array  and 
on  the  utilization  of  its  eigenvectors  and  eigenvalues.  This  matrix  contains 
the  background  noise  spectral  density  matrix  for  which  a  good  parametrized 
model  is  needed.  Performances  estimated  by  simulations  show  on  one  hand  the 
interest  of  high  resolution  methods  compared  to  conventional  methods,  and 
on  the  other  hand  their  sensitivity  to  the  background  noise  model  which  is 
used. 


1  -  INTRODUCTION 

One  of  the  main  functions  of  an  underwater  passive  listening  system  is  the 
determination  of  the  number  of  sources  present  in  the  medium,  as  well  as 
their  characteristic  parameters.  For  that,  one  uses  noises  radiated  by  the 
sources  which  are  received  on  the  sensors  of  an  array.  The  basic  function 
is  array  processing  which  is  traditionally  carried  out  by  conventional 
beamforming.  Effort  towards  performance  improvement  has  led  first  to  adap¬ 
tive  array  processing.  It  provides  a  gain  which  is  limited  by  the  signal  to 
noise  ratio  of  the  sources.  Recently,  studies  have  been  done  on  new  methods 
which  are  more  powerful  0-2-3]  :  they  are  called  high  resolution  methods. 
But  they  need  an  additional  hypothesis  on  the  medium. 

Adaptive  array  processing,  as  conventional  beamforming,  needs  hypothesis 
only  on  the  sources  :  they  are  point  like,  perfectly  spatially  coherent, 
and  the  wavefront  shape  from  a  source  is  a  known  function  of  the  source 
position  (the  transfer  function  of  the  sensors  is  also  supposed  known).  Let 
r(t)  be  the  vector  representing  the  signal  received  on  the  K  sensors  of  the 
array.  The  correlation  matrix  of  the  received  signal  is  : 

C(t)  -  r(t)  r+(t+x)3  (I) 

where  E()  denotes  expectation  and  r  is  the  conjugate  transpose  of  r.  With 
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the  above  hypotheses,  the  spectral  density  matrix,  which  is  the  Fourier 
transform  of  C(x),  of  a  source  alone  is  equal  to  : 


r<f)  =  y(f)  D(f)  D+(f) 


(2) 


where  y(f)  is  the  spectral  density  of  the  received  signal  and  D(f)  the 
source  position  vector,  composed  of  the  transfer  functions  between  the 
source  and  each  sensor,  normalized  by  the  transfer  function  between  the 
source  and  a  reference  point  of  the  array.  This  matrix  is  rank  one,  which 
is  characteristic  of  its  perfect  spatial  coherence. 

The  additional  hypothesis  needed  by  the  high  resolution  methods  concerns 
the  spatial  coherence  of  the  background  noise.  In  general,  background  noise 
is  supposed  to  be  statistically  independent  between  the  sensors  ff|  or 
suppressed  (jf|  .  High  resolution  methods  have  better  performances  than  those 
of  adaptive  array  processing  thanks  to  that  hypothesis.  They  suppose  also 
that  the  noisefield  can  be  resolved  :  the  number  N  of  sources  is  less  than 
the  number  K  of  sensors.  The  spectral  density  matrix  of  the  received  si¬ 
gnals  is  equal  to  (sources  and  background  noise  are  statistically  indepen¬ 
dent)  : 

N 

r (f >  -  r  (f)  +  rc(f)  =  o(f)  i  +  l  y.(f)  1(f)  b.(f)  (3) 

B  s  i«i  1  1  1 

Fg(f)  and  r  (f)  are  respectively  the  spectral  density  matrices  of  the 
background  noise  and  of  the  sources,  I  is  the  spatial  coherence  matrix  of 
the  background  noise  which  is  the  identity  matrix  and  o(f)  its  spectral 
density,  D.(f)  and  y.(f)  are  respectively  the  position  vector  and  the 
spectral  density  of  t/»e  i*h  source. 


2  -  HIGH  RESOLUTION  METHOD  PRINCIPLES. 


High  resolution  methods  are  based  on  eigenvalue-eigenvector  decomposition 
of  the  spectral  density  matrix  P(f).  An  eigenvector  V(£)  and  its  correspon¬ 
ding  eigenvalue  X(f),  are  defined  by  the  relation  : 

N 

r<f>  V(f)  -  3(f)  v(f)  ♦  l  y .  ( f )  5.<f)  r  8  .(f)  v(f)l  -  x(f)  v(f )  <4) 

i-i  1  1  1 

It  can  be  shown  that  F(f)  has  : 


a)  N  eigenvectors  V,(f)  (i  8  0 ,Nj )  that  are  the  eigenvectors  of  the  sour¬ 
ces  alone  raatrixr^f)  corresponding  to  the  N  non-zero  eigenvalues  X  ,(f) 
of  I*(f)  the*  rank  of  which  is  equal  to  N  }  the  corresponding  eigenvalues 
are  equal  to  s  X.(f)  *  X  .(f)  ♦  o(f)  }  these  N  eigenvectors  are  an  ortho¬ 
gonal  basis  of  trie  N  dimensional  source  subspace  spanned  by  the  N  posi¬ 
tion  vectors  D. (f)  of  the  sources  ;  then  : 

N  N 

l  yAi)  1(f)  1(f)  «  l  \  At)  V.(f)  1(f)  (5) 

i-1  1  1  1  i«l  al  1  L 


b)  (K-N)  eigenvectors  V.(f)  orthogonal  to  the  preceding 
to  each  position  vector  1(f)  i 


ones  and  therefore 
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vT(f)  1(f)  =  0 


(6) 


They  are  a  basis  of  the  (K-N)  dimensional  subspace  orthogonal  to  the  source 
subspace.  The  corresponding  (K-N)  eigenvalues  are  all  equal  to  o(f),  and 
therefore  smaller  than  the  previous  ones. 


From  the  above  analysis,  it  is  deduced  that  : 


a)  The  number  K  of  sensors  minus  the  number  of  minimum  and  equal  eigenvalues 
gives  the  number  of  sources. 


b)  The  position  of  the  sources  can  be  determined  by  using  : 


-  either  the  source  subspace  and  relation  (5)  as  proposed  in  ff]  or  (Y| 

-  or  the  orthogonal  subspace  as  proposed  in  | 3]  :  as  the  eigenvectors  of 
this  subspace  are  orthogonal  to  each  source  position  vector,  if  D(f ,0) 
is  a  position  vector  according  to  the  wavefront  shape  of  the  model 
(including  propagation  and  array)  and  corresponding  to  a  scarce  ^he  , 
position  of  which  is  denoted  by  3,  the  scalar  products  :  V.(f)  D(f,0) 
are  equal  to  zero  for  each  value  of  0  corresponding  to  the  position  of 
a  source  in  the  medium.  For  stability  and  ambiguity  reasons  Q)|,  the 
orthogonal  subspace  method  uses  the  function  : 


K-N 

G(f,6)  -  l 

j-i 


i  vT ( t )  D(f ,0) I  2 


(7) 


It  is  that  separation  into  two  subspaces,  one  which  contains  the  sources, 
the  source  subspace,  and  the  other  which  contains  the  background  noise  only, 
the  orthogonal  subspace,  which  gives  to  high  resolution  methods  their  main 
property  .•  their  asymptotic  (infinite  observation  time)  resolving  power  is 
infinite,  that  is  to  say  two  sources  can  be  resolved  even  if  they  are  very 
close  together  and  very  weak  ;  it  is  no  longer  limited  by-  the  input  signal 
to  noise  ratio  y^(£)/o(f)  as  it  is  for  adaptive  array  processing. 


The  above  properties  are  valid  only  for  an  infinite  observation  time.  In 
practice,  the  observation  time  is  limited  ant  ^nly  an  estimation  P(f)  of  t'(f) 
is  available.  Therefore,  only  estimation  of  the  source  parameters  can  be 
obtained.  The  minimum  eigenvalues  of  I* (f )  are  not  exactly  equal  but  show  a 
spread  chat  depends  upon  the  observation  time.  Thus  the  determination  of  the 
number  of  sources  is  an  actual  detection  test  with  detection  and  false  alarm 
probabilities  :  it  must  be  noticed  that  this  detection  test  of  the  number  of 
sources  does  not  need  the  wavefront  shape  knowledge. 


Eigenvectors  are  also  estimates  of  the  actual  ones  and  only  estimates  of  the 
source  positions  are  obtained.  In  particular,  orthogonal  subspace  eigenvec¬ 
tors  are  not  exactly  orthogonal  to  the  source  position  vectors  and  t!(f,0) 
(relation  (7))  exhibits  minimum  instead  of  zeroes.  Fig.  I  shows  an  example 
obtained  by  simulation.  Receiving  array  is  linear,  composed  of  12  equi spa¬ 
ced  sensors.  There  are  two  sources  at  infinity  in  the  medium  with  bearings 
0*  and-5.3°  and  spatially  incoherent  background  noise.  Sources  are  in  the 
same  plane  as  array.  The  figure  shows  versus  bearing  the  output ^spectral 
density  of  adoptive  array  processing  y.  *«d  the  inverse  of  G(f,0),  which 
exhibits  spikes  at  the  source  position  instead  of  minimum.  The  spectral 
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density  matrix  is  estimated  by  : 

P 

r (f )  =  l  l  x. (f)  x!(f) 

i=1 


(8) 


where  X.(f)  is  the  discrete  Fourier  transform  of  the  vector  signal  r(t) 
computed  on  a  time  duration  equal  to  T.  When  the  parameter  P,  that  is  to 
say  the  observation  time,  increases,  the  adaptive  array  response  remains 
about  the  same  while  the  response  of  G(f,6)  is  improved  :  the  spikes  become 
more  narrow. 

The  physical  limitations  to  the  asymptotically  infinite  resolving  power  of 
high  resolution  methods  are  given  by  the  imperfect  knowledge  and  the  fluc¬ 
tuations  of  the  background  noise  spatial  coherence  and  of  the  source  wave- 
front  shape.  Influence  of  the  background  noise  spatial  coherence  is  exami¬ 
ned  now. 


3  -  INFLUENCE  OF  BACKGROUND  NOISE  SPATIAL  COHERENCE  |T| 

In  fact,  at  sea,  the  background  noise  has  not  a  spatial  coherence  matrix 
equal  to  an  identity  matrix.  It  is  composed  of  several  components  :  the  sea 
noise,  generated  by  the  wind  at  the  surface,  the  flow  noise,  generated  along 
the  array,  and  the  traffic  noise  at  the  low  frequencies.  Generally,  the 
electronic  noise  of  the  hydrophone  channels  is  negligible. 

So  in  the  general  case,  the  spectral  density  matrix  of  the  background  noise 
can  be  written  : 

rB(f)  -  o(f)  J(f)  (9) 

where  o(f)  is  the  spectral  density  and  J(f)  the  spatial  coherence  matrix  of 
the  background  noise. 

In  that  case,  the  properties  of  the  eigenvectors  and  eigenvalues  of  the 
spectral  density  matrix  stated  in  the  preceding  section,  are  no  more  valid. 
For  example,  in  order  to  illustrate  the  phenomenon,  the  very  simple  ease 
where  the  background  noise  is  reduced  to  surface  noise  only  QQ  can  be  consi¬ 
dered.  The  cross-spectral  density  between  two  sensors  distant~*of  d  is  equal 
to  : 

V£)  *  o(£>  Jm(2,£d'c>  <,0) 

where  c  is  the  sound  velocity  in  the  sea,  J  <)  the  Bessel  function  of  order 
m,  and  m  a  modelling  parameter.  m 

The  signals  received  on  a  linear  array  of  12  equally  spaced  sensors  have 
been  simulated.  Bearing  only  is  considered.  Table  1  shows  the  eigenvalues 
obtained  for  a  noiseficld  composed  of  background  noise  only  with  a  parame¬ 
ter  m  equal  to  zero,  and  for  a  length  fc  betwen  two  adjacent  sensors  equal 
to  0.3  or  0.7  vavolength.  If  the  spatial  coherence  matrix  of  the  background 
noise  was  an  identity  matrix  according  to  the  hypothesis,  all  the  eigen- 
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values  would  be  equal  :  Table  1  shows  significant:  differences 


£=0.31 

2.45 

2.26 

1.34 

1.29 

1.13 

1.10 

1.05 

1.00 

0.32 

0.05 

0.004 

-4 

2,10 

4-0.71 

2.67 

2.22 

1.38 

1.11 

0.76 

0.61 

0.53 

0.  55 

0.54 

0. 53 

0.52 

0.52 

Table  1  :  Eigenvalues  of  spatially  correlated  background  noise 

Fig.  2  presents  the  diagram  versus  bearing  0  of  the  eigenvector  correspon¬ 
ding  to  the  maximum  eigenvalue  of  the  spectral  density  matrix  estimate  P(f). 
The  noisefield  is  composed  of  background  noise  and  of  one  source  with 
bearing  5.4°  and  signal  to  noise  ratio  :  -15  dB.  The  length  between  t^o 
adjacent  sensors  is  half  a  wavelength.  In  this  case,  the  eigenvector  V^Cf) 
corresponding  to  the  maximum  eigenvalue  of  the  theoretical  spectral  density 
matrix  T(f)  is  equal  to  the  source  direction  vector.  Therefore,  its  diagram 
defined  by  : 

D(e)  «  | D+(f ,0)  VM ( f ) | 2 

is  equal  to  the  classical  beamforming  diagram.  Three  diagrams  are  drawn  on 
Fig.  2  :  they  are  obtained  for  background  noise  the  spatial  coherence  para¬ 
meter  m  of  which  is  equal  to  0.5,  0.2  and  0.8.  For  0.5  the  spatial  coheren¬ 
ce  matrix  of  the  background  noise  is  equal  to  the  identity  matrix  according 
to  the  hypothesis.  The  diagram  is  effectively  good  for  ma0.5.  Diagrams  show 
weak  differences  for  ma0.2  and  a  completely  different  shape  for  111=0.8. 

Fig.  3  presents  results  obtained  with  the  same  noisefield  as  for  Fig.  2,  but 
for  the  orthogonal  subspace  method  ;  expression  (7).  G  (f,0)  is  drawn 
instead  of  G(fA(p  as  in  Fig.  1 ,  in  order  to  see  spikes  instead  of  minimum. 
The  diagrams  G  (f,(l)  has  been  drawn  for  the  same  background  noise  spatial 
coherence  parameters  as  before.  The  best  result  is  obtained  for  m®0.5,  that 
is  to  say  when  the  spatial  coherence  of  the  aetuul  background  noise  is  equal 
to  the  spatial  coherence  of  the  model. 

Fig.  4  presents  results  obtained  also  with  the  orthogonal  subspace  method. 
The  noisefield  is  composed  of  background  noise  and  of  two  sources  with  bea¬ 
rings  :  0°  and  4®,  and  signal  to  noise  ratio  !  -10  dB.  Three  diagrams  G“kf)» 
are  drawn  versus  bearing  0  for  the  same  background  noise  spatial  coherence 
parameters  :  0.5,  0.2  and  0.8.  There  also  the  best  result  is  obtained  for 
m-0.5. 


The  few  results  given  in  that  section  show  clearly  that  the  background 
noise  spatial  coherence  has  to  be  know  as  exactly  as  possible. 


4  -  ADAPTIVITY  TO  BACKGROUND  NOISE  SPAT I At.  COHBKKNCR. 

When  the  background  noise  spatial  coherence  matrix  is  not  equal  to  the 
identity  matrix  : 

-  o(f)  J(f) 

it  can  be  show  Qfj  that  if  J(f)  is  know,  the  problem  can  be  solved.  As 
j'(f)  is  hermit ian,  positive  definite,  it  is  possible  to  define  a  matrix  C(f) 
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such  that  : 

C(f)  J(f)  C+(f)  =  I  (11) 

Therefore,  the  spectral  density  matrix  T(f)  transformed  by  C(f)  is  equal  to 

N 

F.(f)  =  C(f)  r (f )  C+(f)  =  o (f )  C(f)  J(f)  C+(f)  +  l  y.(f)C(f)D.(f)Dt(f)C+(f) 

t  .=1  i  i  i 

<»2> 

As  :  C(f)D.(f)  is  a  vector,  which  can  be  written  D  .(f),  the  transformed 
matrix  F  ,<*)  is  equal  to  : 

C  N 

r  (f)  -  o (f )  I  +  l  y.(f)  D  .(f)  D+.(f)  (13) 

C  X  Cl  Cl 


Therefore,  the  background  noise  has  been  spatially  whitened,  and  the  eigen¬ 
vectors  and  the  eigenvalues  of  the  transformed  spectral  density  matrix  F  (f) 
have  the  properties  stated  in  the  theory  (section  2),  Of  course,  £o  Recover 
the  source  positions,  a  transformed  source  position  vector  model  D  (0,f) 
must  be  used, in  expression  (7)  in  particular  : 


Dc(6,f)  -  C(f)  D ( 8 , f )  (U) 

So,  the  problem  is  solved  when  J(f)  is  known. 

In  fact,  it  has  been  shown  [jb]  that  it  is  not  necessary  to  know  exactly  the 
spatial  coherence  matrix  of  the  background  noise  :  it  is  only  necessary 
that  it  may  be  modal ized  as  a  function  of  unknown  parameters.  In  this  case, 
the  background  noise  spectral  density  matrix  model  can  be  written  s 

rB(f,m)  ■  o(f )  J(f,m)  (15) 

where  o(f)  is  unknown  and  m  represents  the  unknown  parameter. 

Let  C(f,m)  be  the  matrix  such  that  : 

C(f ,m)  J(f ,m)  C+(f,m)  -  I  (16) 


Let  be  the  parameter  values  of  the  actual  background  noise.  If  the  eigen¬ 
values  of  the  transformed  received  signal  spectral  density  matrix  : 

r  (f*o)  «  C(f,o)  T(f)  C*(f,»)  (17) 

C 

4  4 

are  computed  versus  m,  when  m  is  equal  to  m  ,  (K-N)  eigenvalues  become 
equal  and  the  other  N  are  greater.  In  pract?  ce,  of  course,  only  an  estima¬ 
tion  r(f)  of  r(f)  is  available.  If  the  eigenvalues  diagram  of  the  matrix  : 

r  (f,m)  *  c(f,m)  r(f)  e+(f,m)  (is) 

c 

is  drawn  versus  m,  when  m  is  negr  a  ,  a  focusing  spot  is  observed  which  is 
as  much  sharp  as  the  estimate  r(f)°is  better.  The  value  of  m  corresponding 
to  that  focusing  spot  gives  estimates  of  the  background  noise  parameter  m 
and  of  the  number  of  sources  together.  Eigenvectors  and  eigenvalues  corres¬ 
ponding  to  that  value  of  m  arc  used  for  source  parameters  estimation. 
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In  order  to  show  the  ability  of  the  method,  simulations  have  been  conducted', 
using  the  same  array  and  noisefield  as  in  the  preceding  section*,  the  back¬ 
ground  noise  spatial  coherence  model  (expression  (10))  has  one  unknown  para¬ 
meter  :  m.  The  noisefield  is  composed  of  two  sources  with  bearings  :  0°  and 
4°,  and  signal  to  noise  ratio  :  -10  dB  (as  for  Fig.  4),  and  the  background 
noise  spatial  coherence  parameter  is  equal  to  0.5.  Fig.  5  presents  the  plots 
versus  m  of  the  eigenvalues  of  the  theoretical  spectral  density  matrix  (a), 
and  of  the  matrix  ?(f)  estimated  with  P  =  80  (expression  (8))  (b) ,  and 
P=400  (c)  (P  is  the  number  of  instantaneous  spectra  integrated  to  estimate 
?(f)).  The  focusing  spot  happens  for  m=0.5  :  it  becomes  sharper  when  P,  that 
is  to  say  the  observation  time,  is  increasing.  The  two  higher  eigenvalues 
which  indicate  the  presence  of  the  two  sources  are  clearly  seen. 

This  eigenvalue  diagram  is  therefore  an  interesting  mean  to  estimate  the 
background  noise  spatial  coherence  parameters. 


5  -  CONCLUSION 

The  interest  of  high  resolution  methods  is  clearly  demonstrated  by  their 
asymptotical  performances.  But  in  order  to  be  used  in  practice,  it  is  abso¬ 
lutely  necessary  to  have  a  good  parametrized  model  for  the  spatial  coheren¬ 
ce  of  the  background  noise  and  in  particular  of  the  ambient  noise. 
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FIG.  2  ADAPTIVE  ARRAY  OUTPUT  POWER  y  AND  INVERSE  ORTHOGONAL  SUBSPACE 

-7 

HIGH  RESOLUTION  METHOD  RESPONSE  G  (f , 6)  VERSUS  BEARING  0,  FOR 
SPATIALLY  INCOHERENT  BACKGROUND  NOISE  AND  TWO  SOURCES  WITH 
BEARINGS:  0°  AND  -5.3°,  AND  SIGNAL  TO  NOISE  RATIO:  0  dB 


FIG.  2  MAXIMUM  EIGENVECTOR  DIAGRAMS  VERSUS  BEARING  0  FOR  SPATIALLY 
COHERENT  BACKGROUND  NOISE  WITH  SPATIAL  COHERENCE  PARAMETER  m 
AND  ONE  SOURCE  WITH  BEARING,  5.4*,  AND  SIGNAL  TO  NOISE  RATIO: 
•IS  dB.  For  m  »  0.5#  background  noise  is  spatially  incoherent. 
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FIG.  5  EIGENVALUE  DIAGRAMS  FOR  THE  SAME  NOISE  FIELD  AS  IN  FIG.  4  FOR 
THE  THEORETICAL  SPECTRAL  DENSITY  MATRIX:  a) ,  AND  FOR  THE 
SPECTRAL  DENSITY  MATRIX  ESTIMATE  FOR  TWO  VALUES  OF  THE 
OBSERVATION  TINE  (proportional  to  P) :  b)  and  c) . 
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PERFORMANCE  OF  THRtE  AVERAGING  METHODS, 
FOR  VARIOUS  DISTRIBUTIONS 

by 

Albert  H.  Nuttall 
Naval  Underwater  Systems  Center 
New  London,  CT  Oo32U  USA 


ABSTRACT 

The  performance  of  three  averaging  methods,  namely  the  sample  median,  the 
sample  arithmetic  mean,  and  the  sample  geometric  mean,  are  analyzed  in  terms 
of  their  bias,  variance,  and  mean  square  error.  The  bias  and  variance  are 
numerically  evaluated  for  various  parent  distributions  and  plotted  versus  the 
number,  N,  of  data  points  employed  in  the  sample  statistics.  Also,  the 
limiting  behaviors,  as  N  increases  without  limit,  are  derived.  It  is  found 
that  the  best  averaging  method  is  very  dependent  upon  the  distribution  of  the 
data,  with  the  sample  median  being  favorud  for  data  with  occasional  large 
out-1  iers. 


1N1RUUULT1UN 

Estimation  of  average  properties,  such  as  the  average  power  in  a  particular 
angular  sector  and/or  frequency  bin  is  often  accomplished  by  taking  N 
Independent  measurements  of  such  data  and  calculating  a  simple  arithmetic 
average.  However,  when  the  desired  process  is  subject  to  random  fade-outs  or 
occasional  large  out-1 iers,  this  sample  arithmetic  mean  (SAM)  is  severely 
perturbed,  and  alternative  averaging  methods  should  be  considered.  Two 
possible  candidates  are  the  sample  median  (SMU)  and  the  sample  geometric  mean 
(SUM);  these  nonlinear  processors  of  the  available  data  Have  the  potential  of 
suppressing  the  deleterious  effects  mentioned  above.  Here  we  Investigate*  the 
performance  of  all  three  of  these  averaging  methods  in  terms  of  the  number,  N , 
of  independent  data  points  employed  in  the  pertinent  average,  and  the  parent 
distribution  of  the  data.  A  wide  variety  of  distributions  are  considered, 
some  with  parameters  which  allow  for  significantly  different  character  and 
shapes  of  the  governing  probability  functions. 


♦The  basic  analysis,  derivations,  and  programs  are  given  in  Kef.  1. 
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DtFINiTlUNS 


We  have  available  N  statistically-independent  identically-distributed  samples 
(random  variables)  x^,  x>,  . ...  xu  from  some  parent  population  with 
cumulative  distribution  function  P(u)  =  Prob{x<u}  and  probability  density 
function,  PDF,  p(u)  =  P'(u).  The  SAM  of  the  available  measurements  is 

a(N)  =  ^  (xL  +  x2  +  ...  +  xN);  (1) 


the  SMU  is 


qlN)  5  middle  value  of  {xj,  X£,  ....  XjJ,  for  N  odd; 
and  the  SUM  is  (for  non-negative  random  variables) 

glN)  «  (xA  x2  ...  xN)i/N  «  exp - g - 


(2) 


=  a 


^  A  l ogB (x^+  ...  +  A  1  og^ 


for  any  base  B  >  0 
and  scaling  A. 


13) 


The  last  form  in  U)  for  base  B  »  10  goes  under  the  name  of  dil  averaging. 

As  N  tends  to  infinity,  the  sample  quantities  above  tend  to  definite 
(non-random)  limits,  in  particular,  as  N -*•*, 

a(N)— »  arithmetic  mean  °  Jdu  u  p(u> ; 
q(N)-*  median  ■  u^  ,  where  P(ui/2)  *  ^  * 

g(M)  —♦  geometric  mean  *  exp  /Tdu  ln(u)  p(u));  (4) 


where  we  drop  the  prefix  ‘sample'  for  these  deterministic  quantities.  The 
last  result  in  (4)  follows  from  the  exponential  form  of  the  $GN  in  (3).  If 
the  1/2  in  the  median  definition  is  replaced  by  r,  we  have  for  q(tO  the  sample 
quantile  of  order  r  (Ref.  2,  page  181). 

The  three  limiting  quantities  in  (4)  will  generally  r.ot  be  equal.  For 
example,  for  an  exponential  parent  PUP 

p(u)  *  ;*•  exp^  for  u  >  0,  (5) 
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we  have 


arithmetic  mean  =  m; 

median  =  m  1  n  (2)  =  m  . t>y 3 ;  (6) 

geometric  mean  =  me“*  =  m  .662. 

Thus  we  define  the  bias  of  each  of  the  sample  statistics  (l)-(3)  as  the 
difference  between  tneir  mean  value  and  their  asymptotic  value: 


bias  {Sam}  -  a(N]  -  arithmetic  mean; 

Dias  (SML)]  =  qj¥7  -  median; 

\ 

bias  {SUM]  =  gTffy  -  geometric  mean.  ,  •  (7) 

By  virtue  of  tin's  definition,  all  three  biases  will  tend  to  zero  as  ; 

that  is,  all  three  estimators,  (l)-(3j,  are  asymptotically  unbiased,  each  with 
respect  to  its  desired  value  as  given  by  (4),  respectively. 

It  is  then  convenient  to  define- a  normalized  bias,  NB,  for  each  sample 
statistic  as 

NB  IN)  =  H  ,  ‘  (8) 

0 

wliere  o  is  the  standard  deviation  of  parent  PDF  p(u).  Tirt5  scale  factor  of  N 
leads  to  a  non-zero  value  of  the  normalized  bias  for  large  N,  while  the  scale 
factor  of  o  is  convenient  in  that  it  eliminates  the  dependence  of  the 
normalized  bias  on  the  absolute  scale  of  the  input  data.  For  large  N,  (8) 
yields 


bias  -  o  — as  N  ;  ••  (9) 


thus  MB(co)  is  an  Important  measure  of  quality 'of  the  particular  sample 
statistic  under  consideration. 


The  variances  of  sample  statistics  { 1 )- ( 3)  are  defined  as 

var  {Sam]  = 

a*(N)  - 

2  ..  ' 

a(N)  ; 

var{SMD]  = 

q2(U)  - 

_ 2 

q(N)  ; 

II 

S- 

'O 

g*(N)  - 

_ 2 

y(N)  . 

(10) 
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Again,  since  these  quantities  tend  to  zero  for  large  N,  it  is  more  convenient 
to  define  a  normalized  variance,  NV,  as 

NVCN)  lanance  .  (11) 


Then  we  can  state  that 

variance  ~  —  as  (12) 


thus  NV(oo  )  is  also  an  important  measure  of  the  quality  of  a  particular  sample 
statistic. 

We  present  results  here  for  NB(N)  and  NV(N),  along  with  their  asymptotic 
values  at  N  *  <*  ,  for  a  variety  of  parent  distributions  P(u).  Additional 
results  for  the  sample  quantile  with  r  =  .75  and  .9,  and  for  the  PDF, 
cumulative  distribution  function,  characteristic  function,  cumulants,  and 
moments  of  the  various  sample  statistics  are  available  in  Ref.  1. 

RESULTS 

The  first  case  we  consider  is  the  Gaussian  PDF  with  arithmetic  mean  m  and 
variance  a2.  Since  this  random  variable  can  go  negative,  the  SGM  is  undefined. 

The  SAM  and  SMD  are  unbiased  for  all  N;  thus  NB(N)  =  0  for  all  N,  for  this  example. 
Results  for  the  normalized  variance  are  presented  in  Fig.  1  for  the  number  of 
samples,  N,  between  1  and  51,  for  both  the  SAM  and  the  SMD.  The  normalized 
variance  for  the  SMD  is  computed  only  at  odd  values  of  N,  Indicated  by  an  X, 
and  straight  lines  drawn  between  these  points  for  ease  of  association  of 
values.  It  is  seen  that  the  variance  of  the  SMD  is  always  greater  than  that 
for  the  SAM,  the  limiting  value,  NV(oo),  being  */2  for  the  SMD;  see  also  Ref. 

2,  page  369.  Observe  that  the  parameters  m  and  a  of  this  PDF  have  dropped  out 
of  this  normalized  plot. 

It  is  worth  pointing  out  here,  and  for  similar  results  to  follow,  that 
although  the  curve  for  the  SMD  increases  with  N,  that  does  not  mean  that  the 
variance  Increases  with  N;  rather,  the  normalizing  factor  of  N  In  definition 
(U)  causes  this  behavior.  The  actual  (unnonmali zed)  variance  decreases 
monotonlcally  with  N,  eventually  being  of  order  1/N. 

For  a  Rayleigh  randora  variable,  the  results  for  the  normalized  bias  are  given 
In  Fig.  2.  The  SAM  Is  unbiased  for  all  N,  whereas  the  SGM  and  SMD  are,  of 
course,  only  asymptotically  unbiased.  The  limiting  values,  NB(*>),  for  both 
of  these  latter  sample  statistics  are  given  by  analytically  complicated 
expressions  and  are  not  repeated  here,  for  the  sake  of  brevity;  they  are 
Indicated  numerically  by  horizontal  lines  at  the  right  edge  of  the  figure. 

The  corresponding  results  for  the  normalized  variance  are  given  in  Fig.  3. 

They  indicate  that  whereas  the  SGM  has  about  the  same  stability  as  the  SAM, 
the  variance  for  the  SMD  is  about  65% greater. 

For  an  exponential  PDF  (as  given  by  (5)),  the  normalized  bias  and  variance 
results  are  presented  in  Figs.  4  and  5  respectively.  The  biases  of  the  SGM 
and  SMU  are  comparable,  but  we  observe  that  the  variance  for  the  SGM  is  twice 
as  small  as  that  for  the  SAM  and  the  SMD. 
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For  a  log-normal  PDF,  u 


for  u  >  0 


(13) 


the  NB(N)  and  MViM)  results  are  independent  of  location  parameter  my,  but 
they  do  depend  on  spread  factor  oy.  This  may  be  anticipated  by  plotting  the 
PDF  (13)  for  various  values  of  oy  and  observing  that  the  shape  changes  as 
ay  does.  Since  NB(N)  and  NV(N)  depend  upon  the  shape  of  the  PDF  (rather 
than  upon  absolute  location  and  scale),  results  will  depend  on  the  particular 
value  of  ay  selected.  An  example  of  NB(N)  and  NV ( N )  for  ay  =1  is 
presented  In  Figs.  6  and  7.  Now  we  observe  that  the  variance  of  the  SMD  is  3 
times  better,  and  that  of  the  SGM  4.6  times  better,  than  for  the  SAM,  at  least 
for  laryer  values  of  N.  However,  as  oy-»0,  the  log-normal  PDF  in  (13) 
approaches  a  Gaussian  PDF  about  the  point  u  *  exp(my),  and  the  behaviors  would 
revert  back  to  Fig.  1  then. 


The  next  example  is  the  Rice  PDF;  physically,  this  corresponds  to  the 
squared-envelope  of  the  sum*  of  a  sine  wave  and  a  centered  narrowband  Gaussian 
noise  process.  That  is,  the  PDF  is 


p(u)  =  — 

T 


2 

I 


exp 


uJJlA 
2<rf  / 


1 


0 


for  u  >  0 


(14) 


where  A  is  the  sine  wave  amplitude  and  cr,  is  the  noise  standard  deviation. 

Once  ayain,  the  shape  of  the  PDF  depends  on  a  parameter,  namely  A/<r, .  Results 
for  A/<r,  =  1  are  given  in  Flys,  8  and  9.  The  SMD  has  147*  greater  variance 
than  the  SAM,  but  the  SGM  has  about  607*  of  the  SAM  variance.  As  A/oj-*0,  the 
exponential  PDF  results  are  obtained,  whereas  as  A/a;-*«>,  the  Gaussian  case  is 
realized.  Thus  (14)  represents  a  transition  case  between  these  extremes. 

The  last  example  we  consider  here  Is  an  exponential  PDF  with  out-liers.  That 
is,  each  sample  or  measurement  {xj^  in  (l)-(3)  is  given  by 

X|c  B  xa  +  xp,  (IS) 

where  xa  has  an  exponential  PDF, 

pd(u)  °  ~  exp  )  for  u  >  0,  (16) 


*lf  each  of  the  observed  random  variables  xi,  x;>,  ...,  xn  in  (l)-(3)  is 
obtained  by  first  summing  up  the  envelope-squared  outputs  of  M  narrowband 
filters,  as  for  example  In  diversity  reception,  the  PDF  in  (14)  is  replaced  by 
the  Q*  distribution.  Results  for  this  case  are  available  in  Ref.  1,  but  are 
not  given  here,  for  sake  of  brevity. 
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and  disturbance  xb  is  a  random  variable  which  is  zero  most  of  the  time,  but 
occasionally  takes  on  a  large  value  (out-lier)  L.  That  is,  its  PDF  is 

Pb(u)  =  (1  -  QJ  S(u)  +  Q  $(u  -  L), 

where  Q  is  the  probability  of  an  out-lier.  Then  the  parent  PDF  of  observation 
(sample)  x^  is  the  convolution  of  (16)  and  (17): 

P(u)  exp  (-£-)  ll(«)  +  jj-  exp  (-  U(u  -  L),  (18) 

a  ^  a  a  \  a  ' 

where  unit  step 

f 1  for  t  >  (f 
U(t)  =  < 

l 0  for  t  <  0_ 

The  important  parameter  now  is  L/ma,  which  obviously 
PDF  (18), 

Now  however,  before  we  get  into  the  detailed  bias  and  variance  results, 
another  consideration  is  of  paramount  importance.  Our  sample  statistics, 

(l)-(3),  will  no  longer  extract  (estimate)  the  arithmetic  mean,  median,  and 
geometric  mean,  respectively,  of  the  (disturbance-free)  exponential  PDF  (16), 
but  perforce,  the  corresponding  statistics  of  the  measurement  PDF  (18).  If, 
however,  we  are  really  interested  in  the  parameters  of  (16),  then  we  must 
Inquire  into  the  quantitative  disturbance  caused  by  the  out-llers  described  in 
(17).  Here  we  merely  cite  the  results  for  one  numerical  case;  additional 
results  are  given  In  Ref.  1. 

For  probability  Q  .05,  and  out-lier  value  L/m  =  6,  we  find  that  the  ratio  of 
arithmetic  means,  for  (18)  with  respect  to  (16),  Is  1.3.  The  corresponding  ratio  of 
medians  is  only  1.08,  whereas  the  ratio  of  geometric  means  is  1.13,  Thus  the 
SMD  and  SUM  are  more  resistant  to  the  presence  of  Infrequent  out-llers, 

Insofar  as  their  effects  on  the  particular  parameters  of  median  and  geometric 
mean. 

The  results  for  the  normalized  bias  and  variance  are  given  in  Figs.  10  and 
11.  The  variances  of  the  SMD  and  SUM  are  again  smaller  than  that  for  the 
SAM.  The  bias  of  the  SMD  and  SUM  are  comparable. 

If  we  define  a  mean-square  error  as  the  average  value  of  the  squared  > 
difference  between  a  sample  statistic  s  and  a  desired  parameter  da  of  the 
disturbance-free  PDF,  we  can  develop  it  as  follows: 


mse 


» +  (?  -  day  .  (a) 


(19) 

affects  the  shape  of 
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But  the  first  term  is  the  variance  of  the  sample  statistic,  and  the  second 
term  can  be  expressed  as 

s  '•  da  =  is  -  d)  +  (d  -  da)  =  bias  +  deflection  in  desired  parameter,  (21) 


where  d  is  the  modified  value  of  the  desired  parameter  da,  due  to  the 
disturbance.  Thus 

MSE  =  variance  +  (bias  +  deflection)^  .  (22) 

Now  the  bias  and  variance  are  0(N  ^ )  for  large  N,  whereas  the  deflection  of  the 
desired  parameter  does  no*-  decay  with  N  at  all;  in  fact,  it  is  independent  of  N. 
Thus  the  considerationsTB’ove,  whether  for  the  ratio  of  arithmetic  means  or 
medians  or  geometric  means,  are  very  important,  since  they  dominate  the 
magnitude  of  the  mean-square  error  for  very  many  samples  available. 

U1SCUSSI0N 


The  ability  of  the  SMD  and  SGM  to  suppress  deleterious  effects  due  to  occasional 
large  interferences  is  very  pronounced  for  some  probability  density  functions. 
Not  only  is  the  deflection  of  the  desired  parameter  (arithmetic  mean  or 
median  or  geometric  mean)  decreased,  but  the  bias  and  variance  of  the  estimate 
can  bo  markedly  reduced  in  some  cases.  The  exact  amounts  depend  on  the 
magnitude  and  frequency  of  the  i nterference. 

Another  possible  approach  to  alleviate  the  effects  of  additive  large  out-1 iers 

is  to  subject  the  available  samples  xi,  x^- . .  to  a  nonlinear 

transformation  such  as  saturation,  in  order  to  suppress  the  large  contri¬ 
butions,  prior  to  evaluating  the  SAM  or  SGM  or  SMD.  Knowledge  of  relative 
levels  (such  as  L/ma  for  the  above  example)  would  be  required  for  optimal 
adjustment  of  the  saturation  level,  but  performance  could  be  markedly  improved. 
The  nonlinear  transformation  would  reduce  the  deflection,  while  the  averaging 
of  N  samples  would  reduce  the  bias  and  variance.  This  possibility  has  not  yet 
been  pursued. 
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FIG.  9  NORMALIZED  VARIANCE ;  RICE  PDF  (A/a  f  1) 
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IS  POWER  AVERAGING  THE  BEST  ESTIMATOR 
FOR  UNDERSEA  ACOUSTIC  DATA? 


By 

R.A.  Wagstaff  and  J.L.  Berrou 
SACLANT  ASW  Research  Centre 
La  Spezia,  Italy 


ABSTRACT 


The  power  average  is  the  statistic  which  is  commonly  used  in  all  types  of 
acoustic  data  processing.  It  is  an  excellent  estimator  for  "well  behaved" 
Gaussian  distributions  when  the  central  or  average  value  is  desired. 
However,  when  there  are  "outliers"  and  errors  or  when  the  distributions  are 
not  Gaussian  there  are  other  statistics  which  are  better  estimators.  In 
the  case  of  ambient  noise,  the  distributions  are  seldom  Gaussian  or  free 
from  "outliers"  and  errors.  For  this  case  the  median  and  the  geometric 
mean  power  level  (dB  average)  are  less  sensitive  to  interference  in  the 
environment  and  to  some  of  the  minor  system  faults  and  are  more  powerful 
statistics. 


INTRODUCTION 


It  is  not  a  simple  task  to  measure  the  undersea  ambient  noise  and  to 
process  and  report  tne  results  with  a  high  level  of  confidence.  The 
measurement  itself  is  plagued  with  many  potential  hazards.  The  equipment 
is  seldom  in  perfect  condition  and  the  noise  environment  is  usually  not 
cooperative.  Nearby  ships,  seismic  exploration,  and  countless  other 
interferences  have  ruined  many  hours  of  otherwise,  good  data.  The 
measurement  system  has  its  own  problems  which  incre-.se  the  difficulty  of 
achieving  good  results.  Hydrophones  can  go  bad  as  well  as  can  the 
electronic  components  in  amplifiers,  preamplifier*-.,  and  filters  which 
condition  the  acoustic  signals.  If  the  system  i<  n  vjwed  array  there  is 
the  additional  problem  of  towship  noise.  In  some  cases,  the  noise  from  the 
towship  is  considerably  more  than  the  ambient  which  is  to  be  measured. 

Past  experience  in  the  measurement  of  ambient  noise  in  many  different 
areas,  with  different  environmental  conditions  and  noise  source 
distributions,  and  with  different  measurement  systems  indicates  that 
ambient  noise  data  are  generally  not  Gaussian  distributed.  The 
distributions  tend  more  toward  the  shape  of  a  log-normal  or  Rayleigh 
distribution.  The  distribution  function  in  levels  has  a  longer  tail  on  the 
low  noise  level  end  than  on  the  high  noise  level  end.  Host  of  the  usual 
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interferences,  however,  are  distributed  more  toward  the  high-level  end. 
This  is  intuitively  obvious,  since  an  interference  of  low  level  noise 
cannot  reasonably  interfere  with  the  measurement. 

The  average  power  level  is  a  commonly  used  statistic  in  ambient  noise  data 
processing.  Because  of  the  extremely  high  level  of  the  interference,  the 
average  power  level  can  be  dominated  even  when  the  interference  constitutes 
a  relatively  small  percentage  of  the  data.  If  a  "persistant  background" 
noise  is  desired  from  the  measurement,  the  average  power  will  not  provide 
it.  This  is  also  true  if  trying  to  detect  a  target  during  times  of  high 
level  noise  interference,  which  is  illustrated  by  the  time  series  in 
Fig.  1.  Two  explosions  were  received  from  seismic  prospecting  during  the 
measurement  period.  Three  statistics  have  been  calculated:  the  average 
power  level,  the  decibel  average  (dB  AVG),  and  the  median  level  (MEDIAN). 
In  this  case,  the  average  power  level,  being  biased  upward  about  5  dB,  is 
not  a  good  estimator  for  what  might  be  considered  representative  of  the 
data.  The  other  two  seem  reasonable.  A  target  line  would  be  masked  if  the 
averaging  period  contained  one  or  more  of  these  transient  noise  "spikes". 

The  selection  of  a  best  statistic  is  based  on  a  very  simple  concept:  for  a 
source  to  interfere  it  must  be  of  reletively  high  level.  A  lower  level 
source  does  not  interfere.  The  ideal  statistic,  then,  is  one  which 
discriminates  against  high  levels.  This  would  not  be  the  average  power 
level.  The  average  power  can  easily  be  captured  and  biased  upward  as  in 
Fig.  1  when  the  high  level  interference  constitutes  a  small  percentage  of 
the  data.  However,  it  is  this  statistic  which  is  most  commonly  used  in 
ambient  noise  processing.  The  geometric  mean  power  level  (dB  average)  and 
median  levels,  on  the  other  hand,  are  relatively  insensitive  to 
interferences  which  constitute  a  small  percentage  of  the  data.  A  simple 
example  illustrates  this  point.  Consider  a  set  of  50  observations.  Let 
45  be  of  60  dB  level  and  5  of  100  dB  level.  The  average  power  level  of  the 
50  samples  is  about  90  dB.  The  dB  average  is  64  dB  and  the  median  is  60 
dB.  The  differences  between  these  latter  two  and  the  average  power  level 
are  about  26  and  30  dB.  If  It  is  accepted  that  the  five  samples  of  100  dB 
are  a  result  of  interference,  then  the  dB  average  and  median  are  better 
statistics  than  the  average  power  level.  These  statistics  have  the  obvious 
advantage  that  they  Qive  very  good  estimates  without  the  need  for  deleting 
suspect  data.  This  is  important  for  automatic  data  processing.  The 
advantage  of  course,  decreases  as  the  percentage  of  the  contaminated  data 
Increases.  Some  examples  will  be  given  from  measured  data  which  illustrate 
the  "power"  of  these  two  statistics  for  ambient  noise  measurement.  The 
same  concepts  apply  to  other  processes  in  acoustic  measurement  and  signal 
processing.  It  will  be  left  as  a  challenge  to  the  reader  to  discover  In 
which  phases  of  analysis  this  simple  concept  can  be  applied  to  his 
situation. 


RESULTS 

During  one  ambient  noise  measurement  exercise  with  the  towed  array, 
interference  from  seismic  prospecting  wes  received  during  the  entire 
measurements.  The  range  to  the  source  was  not  known,  only  the  azimuth. 
The  level  of  the  interference  at  the  array  was  'bout  95  dB  at  480  Hz.  The 
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repetition  rate  was  between  9  and  12  seconds  and  each  explosion  lasted 
about  2  seconds.  The  480  Hz  ambient  noise  without  the  interference  was 
approximately  66  dB.  The  effects  of  these  explosive  sources  on  the  data 
processing  are  illustrated  in  Fig.  2.  There  are  three  curves  in  each  plot 
of  Fig.  2.  The  dashed  curves  at  the  top  are  the  average  power  levels  for 
50  spectral  estimates.  The  next  two  curves  are  the  for  the  median  (solid) 
and  the  dB  average  (dashed)  for  the  same  50  data  points. 

The  explosions  cause  an  interference  pattern  across  the  array  in  the 
hydrophone  noise  power  average  (top  dashed  curve).  The  mean  level 
decreases  about  5  dB  across  the  array  and  has  oscillations  with  an 
amplitude  of  about  one  decibel.  Such  effects  could  cause  problems  to  many 
types  of  signal  processing.  The  curves  of  the  median  (solid  curve)  and  dB 
average  (bottom  dashed  curve)  register  considerably  different  behavior. 
The  median  has  about  a  1.5  dB  slope  across  the  array  with  very  little 
evidence  of  oscillations  and  the  dB  average  is  flat  with  no  oscillations. 


The  reasons  the  median  and  the  dB  average  give  such  excellent  results  are 
the  following.  The  data  points  with  the  explosions  are  all  above  the  median 
level.  About  half  of  them  might  have  been  below  the  median  if  there  had 
been  no  explosions.  Therefore,  the  median  level  is  shifted  in  the  ordering 
of  the  data  by  an  amount  equal  to  about  half  the  number  of  explosions.  If 
4  explosions  were  received  in  the  50  samples  the  shift  would  be  about  two 
numbers  but  not  more  than  four.  The  resulting  change  in  the  level  of  the 
median  would  be  small  and  possibly  none  at  all.  As  the  number  of 
occurrences  increases  in  proportion  to  the  total  data,  the  affect  on  the 
median  increases.  The  same  is  true  for  the  dB  average.  However,  the 
discrimination  capability  of  the  dB  average,  in  addition,  is  affected  by 
the  relative  magnitude  of  the  difference  between  the  good  data  and  the 
interference.  Estimates  of  the  beam  noise  level  and  the  noise  field 
horizontal  directionality  could  be  made  from  the  median  and  dB  average 
beam  noise  levels;  it  could  not  be  done  from  the  average  power  levels 
without  being  biased  by  the  explosions.  The  results  in  the  beam  noise  plot 
Fig  2b  indicate  the  explosion,  near  beam  number  13,  which  dominates  the 
curve  for  the  average  power  levels  (top  dashed  line),  but  is  not  evident  in 
the  curves  for  the  dB  average  and  median.  The  discrimination  level  of  the 
latter  two  is  approximately  34  dB  against  the  explosion. 

Figure  3  illustrates  the  three  different  statistics  on  data  from  two  other 
time  periods.  In  each  case  the  data  base  consists  of  sets  of  50  beam  noise 
spectral  estimates.  Figure  3a  is  for  1460  Hz  when  interference  was  being 
received  from  a  source  which  was  believed  t.o  be  of  biological  origin, 
possibly  snapping  shrimp.  The  plot  of  average  beam  power  levels  shows  the 
obvious  effects  of  the  interference.  Judging  from  the  other  two  curves 
(median  and  dB  average)  in  the  same  plot,  the  average  power  level  has  been 
biased  up  by  as  much  as  15  dB.  The  curves  for  the  median  (solid)  and  dB 
average  (lower  dashed  curve)  appear  unaffected,  with  the  anisotropy  due  to 
shipping  noise  being  clearly  evident. 


The  data  for  Fig.  3b  were  obtained  when  a  different  type  of  biological 
source  was  nearby.  The  noises  received  on  the  array  sounded  like 
carpenters  hammering  on  a  roof.  The  sources  of  the  noise  were  most  likely 
carpenter  fish  (Sperm  whale).  The  average  power  levels  (top  dashed  curve) 
are  again  "captured",  being  biased  upward  as  much  as  17  dB,  while  the 
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RESULTS  DURING  SEISMIC  PROSPECTING 


FIG.  2  NOISE  LEVEL  VERSUS  HYDROPHONE  NUMBER  PLOT  (2a)  AND  BEAM  NOISE  LEVEL 
VERSUS  BEAM  NUMBER  PLOT  (2b)  FOR  A  PERIOD  WHEN  NOISE  FROM  SEISMIC 
PROSPECTING  WAS  BEING  RECEIVED.  POWER  AVERAGE  (top  dashed  curve) , 
MEDIAN  (solid  curve),  AND  dB  AVERAGE  (bottom  dashed  curve) 


BIOLOGICAL  NOISES 


FIG. 


3  NOISE  \BVEL  VERSUS  BEAM  NUMBER  FOR  TIMES  OF  NOISE  INTERFERENCE 
FROM  <:*)  SNAPPING  SHRIMP ,  AND  (3b)  CARPENTER  FISH  (Sperm  whale). 
POWER  AVERAGE  (dashed  curve) ,  AND  dB  AVERAGE  (bottom  dashed  curve) 
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curves  for  the  medians  and  dB  averages  do  not  register  any  significant 
effects.  Hence,  in  both  sets  of  data  in  Fig.  3  the  background  shipping 
no'se  could  be  measured  by  the  median  and  the  dB  average  during  the  periods 
of  interference,  while  it  could  not  be  measured  with  the  power  average. 
This  might  also  have  been  true  if  instead  the  objective  was  to  detect  a 
target.  The  gain  against  the  interference  is  the  level  differences  between 
the  power  average  and  the  other  two  statistics.  The  potential  improvement 
is  substantial . 

The  previous  examples  demonstrate  that  the  median  and  the  dB  average  can 
effectively  discriminate  against  high  level  transient  sources  in  the 
environment.  Unfortunately,  the  environment  is  not  the  only  source  of  bad 
data.  The  measurement  system  can  also  be  a  problem.  Most  of  the  present 
systems  contain  many  components  and  are  highly  sophisticated.  As  the  number 
of  components  increases,  the  probability  of  experiencing  a  mechanical  or 
electronic  failure  also  increases.  In  fact,  past  experience  indicates 
there  is  always  a  defect  in  the  measurement  system.  It  is  a  challenge  to 
minimize  the  deleterious  effects  of  the  defect  and  produce  high  quality 
results. 

The  median  and  the  dB  average  can  also  be  used  to  improve  the  data  quality 
for  some  types  of  system  faults.  Consider,  for  example,  the  two  plots  in 
Fig.  4.  The  plot  on  the  left  is  a  strip  chart  recording  of  a  hydrophone 
output  from  a  towed  array  during  a  time  of  acoustic  measurements.  Noise 
"spikes"  of  about  18  dB  above  the  background  occur  about  every  10  seconds. 
The  spikes  are  believed  to  be  caused  by  the  relative  motion  of  the 
hydrophone  in  its  fasteners.  The  affect  the  transient  noise  spikes  in  Fig. 
4a  have  on  data  quality  is  illustrated  in  the  beam  noise  plot  of  Fig.  4b. 
The  beams  from  22  to  42  are  steered  into  real  space  and  receive  acoustic 
energy  on  their  main  lobes.  The  beams  outside  this  range  correspond  to 
phase  shifts  greater  than  for  endfire  beams.  They  are  virtual  beams  and 
receive  acoustic  energy  only  through  the  sidelobes  which  extend  into 
acoustic  space  (see  reference  1  for  more  detail).  When  the  sidelobes  are 
good  the  levels  of  these  virtual  beams  are  measures  of  the  self  noise  of 
the  system.  For  the  case  in  Fig.  4b  the  sidelobes  were  in  the 
neighborhood  of  40  dB  down  at  the  time  of  this  measurement  but  were  masked 
by  the  “system"  noises  in  Fig.  4a. 

The  self  noise  resulting  from  the  noises  illustrated  in  Fig.  4a  are  given 
by  the  levels  of  the  virtual  beams.  Relative  to  the  plot,  the  self  noise 
is  about  38  dB  for  the  power  average  and  about  28  dB  for  the  dB  average. 
The  median  is  a  few  decibels  below  the  dB  average.  These  results  indicate 
that  for  the  power  average  the  self  noise  of  the  system  is  about  4  to  6  dB 
below  the  acoustic  levels  on  the  beams.  The  lowest  average  power  beam 
levels  would  most  certainly  be  contaminated  by  the  hydrophone  noise  spikes. 
The  median  and  dB  average  beam  levels,  on  the  other  hand,  have  a  greater 
safety  range.  The  differences  between  the  acoustic  and  virtual  beams  is 
generally  in  excess  of  12  dB  for  these  latter  two  statistics. 

The  dB  average  and  median  can  be  used  in  the  spatial  domain  as  effectively 
as  In  the  time  domain.  For  example,  consider  the  problem  of  estimating  the 
omnidirectional  level  and  horizontal  directionality  of  the  ambient  noise 
from  the  beam  outputs  of  a  towed- line  array.  There  are  two  major  problems 
which  must  be  solved  to  do  this.  The  first  is  that  the  line  array  has  an 
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FIG.  4  STRIP-CHART  RECORDING  OF  NOISES  ON  ONE  HYDROPHONE  CHANNEL  (4a)  AND 
THE  EFFECTS  ON  THE  BEAMFORMED  OUTPUT  (4b).  POWER  AVERAGE  (top 
dashed  curve) ,  MEDIAN  (solid  curve) ,  AND  dB  AVERAGE  (bottom  dashed 
curve) 
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inherent  left-right  ambiguity  which  must  be  resolved.  This  requires  that 
measurements  be  made  on  different  array  headings.  Unfortunately,  the 
noise-field  changes  over  the  time  required  to  do  this.  Hence,  an  exact 
solution  is  generally  impossible.  The  second  problem  is  caused  by  the 
towship.  The  noise  received  from  it  interferes  with  the  measurement  of  the 
ambient.  It  can  be  as  much  as  30  dB  above  the  ambient  for  the  entire  set 
of  measurements. 

The  capability  of  the  dB  average  and  the  median  to  discriminate  against  the 
high  level  transient  noise  can  also  be  used  to  discriminate  against  the 
towship  noise.  When  the  measurements  are  performed  on  different  array 
headings,  the  towship  noise  appears  on  a  different  true  azimuth  each 
measurement  period.  This  constitutes  a  spatial  transient.  Flence,  the 
transient  discrimination  capability  of  the  median  and  the  dB  average  can  be 
utilized  to  discriminate  against  the  towship  in  the  estimation  of  the 
horizontal  directionality  of  the  noise  field.  Furthermore,  the 
ambiguities  in  the  beam  noise  measurements  will  be  resolvea  in  the  same 
way.  The  true  azimuth  of  a  fixed  source  is  not  transient.  The  ambiguous 
azimuth,  however,  will  vary  with  the  array  heading  and  this  constitutes  a 
spatial  transient.  It  will,  therefore,  be  discriminated  against  by  the 
median  and  the  dB  average  performed  in  the  spatial  domain.  V/hen  the 
source  does  not  remain  on  a  fixed  azimuth  for  the  entire  set  of 
measurements,  all  true  and  ambiguous  azimuths  will  be  spatial  transients. 
The  median  and  dB  averages  will  arrive  at  values  which  will  be  between  the 
extremes.  They  will  be  good  estimators  of  the  noise  field  most  of  the 
time.  The  omnidirectional  level  can  then  be  obtained  by  integrating  the 
horizontal  directionality  pattern.  Because  of  the  discrimination 
capability  of  the  median  and  the  dB  average,  the  omnidirectional  level  of 
the  ambient  noise  can  be  estimated  exceptionally  well  in  the  presence  of 
masking  by  the  towship  and  other  high  level  moving  or  transient  sources. 

The  spatial  discrimination  capability  of  the  median  and  the  dB  average  can 
be  illustrated  by  estimating  the  horizontal  directionality  and 
omnidirectional  level  of  the  ambient  noise  for  various  levels  of 
interference  from  the  towship.  This  was  done  for  simulated  measurements  in 
tho  Alboran  Sea  south  of  the  shipping  lane.  The  shipping  surveillance 
data,  obtained  during  the  actual  measurements  were  used  in  the  RANDI  II 
noise  model  <2>  and  the  positions  dead-recknoned  from  one  measurement 
period  to  another.  The  measurements  were  made  each  hour  for  nine  hours  on 
nine  array  headings  which  were  incremented  80  deg.  each  hour.  Noise  was 
added  to  the  beams  normally  contaminated  by  the  towship  and  the  measurement 
was  simulated  again.  This  was  done  for  towship  noise  levels  received  at 
the  array  equal  to  74  dB  and  again  for  114  dB.  This  latter  level  was 
about  40  dB  above  the  ambient  and  70  dB  above  the  lowest  beam  noise  levels 
(10  kt  of  wind  at  480  Hz).  The  algorithm  in  reference  3  was  used  to  obtain 
estimates  of  the  horizontal  directionality  and  the  omnidirectional  level  of 
the  ambient  noise  from  the  beam  noise  data  for  each  level  of  towship  noise 
contamination.  This  algorithm  utilizes  the  dB  average  in  the  deconvolution 
of  the  beam  data  by  an  iterative  approach  and  either  the  dB  average  or  the 
median  in  space  to  resolve  the  ambiguities  and  discriminate  against  the 
towship. 

Example  polar  plots  of  the  beam  data  for  three  of  the  nine  headings  are 
given  in  Fig.  ba.  Only  a  few  of  the  forward  beams  are  affected  by  the 
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towship  noise  as  is  evident  by  the  superposition  of  the  beam  noise  plots 
for  the  cases  of  no  towship  noise,  74  dB  towship  noise,  and  114  dB  towship 
noise.  The  sidelobes  have  been  reduced  sufficiently  to  eliminate  the 
system  as  a  source  of  degraded  results  in  order  to  concentrate  only  on  the 
discrimination  capabilities  of  the  dB  average  and  the  median. 

The  directionality  estimates  and  the  corresponding  omnidirectional  levels 
are  given  in  Fig.  5b.  The  dashed  curve  in  each  plot  was  obtained  directly 
from  the  model.  It  is  the  unambiguous  estimate  (ideal)  of  the  time 
averaged  and  spatially  smothed  noise  field.  It  is  not  degraded  by  the 
towship  noise  or  the  ambiguities  of  the  array.  It  is  the  quantity  that  the 
algorithm  <3>  is  attempting  to  estimate.  The  solid  curve  in  each  plot  is 
the  estimated  directionality  pattern  using  the  median.  The  dotted  curve 
corresponds  to  the  dB  average  processing.  The  plots  from  left  to  right  are 
for  no  noise  from  the  towship,  74  dB,  and  114  dB  from  the  towship.  The 
omni-directional  levels  corresponding  to  each  pattern  are  also  included 
beside  each  plot. 

In  each  of  the  three  cases  illustrated  in  Fig.  5b  the  median  estimate 
agrees  well  with  the  ideal  pattern  except  in  the  90  deg.  sector  south  of 
southeast.  In  this  region  the  estimated  levels  are  biased  up  by  about  7 
dB.  The  azimuthal  extent  of  this  region  increases  with  the  level  of  the 
contamination.  The  omnidirectional  level  of  the  median  estimate  is  biased 
up  about  1.4  dB  when  the  contamination  is  41  dB  above  the  ambient.  This  is 
rather  impressive  suppression  of  the  contamination.  The  dB  average  does 
not  perform  as  well  but  still  gives  impressive  results.  The 
omnidirectional  levels  of  the  dB  average  estimates  are  biased  up  about  5  dB 
for  the  114  dB  contamination  and  1.5  dB  for  74  dB  contamination.  The 
excellent  agreement  between  the  patterns,  especially  for  the  first  two 
plots  of  Fig.  5b,  indicates  that  the  dB  average  and  the  median  are  very 
effective  in  resolving  the  ambiguities  and  estimating  the  omnidirectional 
levels  and  the  directionalities.  The  last  two  plots  indicate  that  these 
statistics  are  also  very  good  at  discriminating  against  spatial  transients, 
in  this  case  the  towship. 

There  is  always  a  question  of  the  validity  of  modeled  results.  The 
previous  results  are  supported  by  measurements  of  a  similar  nature  obtained 
in  the  Alboran  Sea.  The  level  of  the  towship  noise  was  increased  by 
switching  on  the  active  rudder  of  the  towship.  This  increases  the  noise 
received  from  the  towship  by  up  to  30  dB.  The  exact  amount  depends  on 
frequency.  It  is  negligible  at  100  Hz  and  30  dB  at  1000  Hz.  Ambient  noise 
measurements  were  made  on  nine  different  headings.  Data  were  acquired  on 
each  heading  with  the  active  rudder  off  and  another  set  with  it  on.  The 
data  were  processed  as  two  separate  nine-sided  polygons,  one  with  the 
active  rudder  noise  and  one  without. 

Figure  6  illustrates  the  results  for  750  Hz.  The  top  row  of  plots  are  beam 
noise  measured  during  three  of  the  nine  legs.  The  effective 
omnidirectional  level  was  about  73  dB.  The  noise  from  the  towship  exceeds 
the  ambient  beam  noise  from  5  to  20  dB  depending  on  array  orientation.  The 
second  row  of  plots  are  for  the  same  orientations  but  the  active  rudder  is 
switched  on.  In  these  cases  the  noise  from  the  towship  exceeds  the  beam 
noise  from  the  ambient  by  30  to  40  dB.  The  omnidirectional  levels  were 
about  96  dB,  23  dB  higher  than  without  active  rudder.  The  active  rudder 
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measurements  for  each  leg  were  immediately  preceeding  or  following  the 
measurements  without  the  active  rudder  on. 

The  final  row  of  plots  are  the  ambient  noise  directionality  estimates  from 
the  beam  noise  data  for  9  legs.  The  pattern  on  the  left-hand  side  was 
obtained  from  the  data  sets  without  the  active  rudder  and  the  beams 
containing  the  towship  noise  were  eliminated  from  the  processing.  Hence, 
this  pattern  is  completely  free  of  contamination  by  the  towship.  It  will 
be  used  to  judge  the  effectiveness  of  the  dB  average  and  the  median  in 
discriminating  against  the  towship  noise.  The  omnidirectional  level, 
printed  below  the  plot,  is  70.7  dB.  The  center  plot  is  the  directionality 
estimate  from  the  active  rudder  data  when  dB  averaging  s  used.  The 
omnidirectional  level  is  biased  up  7  dB  and  the  low  levels  in  the  southern 
half  space  have  not  been  well  reproduced.  The  results  for  the  median  in 
the  third  plot  are  much  better.  The  omnidirectional  level  is  biased  up  only 
3  dB  and  the  low  levels  have  been  reproduced  much  better.  The  median  is 
clearly  a  more  effective  discrimination  than  the  dB  average. 

The  lowest  levels  achieved  by  the  dB  average  and  the  median  are  not  as  low 
as  when  the  active  rudder  is  not  used  and  the  towship  noise  is  eliminated 
in  the  processing.  The  main  reason  this  is  so  is  that  the  sidelobes  of  the 
array  were  not  sufficiently  suppressed  during  the  measurements.  The 
differences  between  the  active  rudder  noise  and  the  lowest  beams  were  about 
46  dB  (96-50).  This  level  of  sidelobe  suppression  was  not  always  achieved 
during  the  measurements. 

As  an  example  of  discrimination  in  both  the  temporal  and  spatial  domains 
consider  the  ambient  noise  horizontal  directionality  results  which  were 
obtained  by  a  towed  array  at  one  site  in  the  northeast  Atlantic.  During 
most  of  the  measurement  periods  interference  was  being  received  every  9  to 
XI  seconds  at  levels  in  excess  of  30  dB  above  the  normal  beam  levels.  The 
hydrophone  time  series  data  in  Fig.  i  and  the  hydrophone  and  beam  noise 
plots  in  Fig.  2  are  from  the  same  measurement  time  period  but  for 
different  frequencies.  The  difficulty  of  estimating  the  ambient  from  the 
beam  data  was  increased  by  a  sea  state  of  about  ij.  In  such  a  case,  the 
towship  noise  is  the  dominant  source,  not  considering  the  seismic 
prospecting  noise.  It  is  much  easier  to  estimate  the  ambient  at  higher  sea 
states  when  the  towship  noise  is  a  smaller  part  of  the  total  and  less 
discrimination  is  required. 

This  case  is  one  In  which  there  was  interference  in  the  time  domain  and 
spatial  domain.  The  time  domain  interference  consisted  of  explosive 
sources  being  received  about  every  9  to  12  seconds  from  an  azimuth  of 
approximately  333°.  The  spatial  domain  interference  consisted  of  the 
towship  on  the  forward  beams.  these  two  sources  dominated  the 
measurements.  The  seismic  prospecting  noises  were  discriminated  against  in 
the  time  domain  by  using  the  median  beam  outputs  and  the  towship  noise  was 
discriminated  against  by  using  the  dii  average  in  the  spatial  domain.  The 

resulting  directionality  patterns  for  500  to  1500  Hz  are  included  in 

Fig.  7.  These  patterns  do  not  show  evidence  of  either  the  seismic 

prospecting  noise  or  the  towship  noise.  They  are  obviously  dominated  by 

distant  shipping.  The  two  peaks  in  the  northeast  quadrant  are  along 
azimuths  to  a  Port  (72°)  and  a  shipping  lane  (52°).  Without  having 
employed  the  median  and/or  dB  average  in  both  spatial  and  temporal  domains, 
it  is  unlikely  that  such  excellent  results  could  have  been  achieved. 


SACLANTCEN  CP- 32 


17-12 


i 


WAGSTAFF  &  BERROU:  Power  averaging  the  best  estimator? 


In  general  the  towed  array  has  a  variable  geometry.  As  it  is  towed  along 
it  snakes  through  the  water.  The  deviation  from  linearity  often  increases 
with  increasing  sea  state  and  wind  conditions.  One  of  the  effects  of  this 
on  the  acoustic  performance  is  a  degradation  in  the  beamforming  capability 
of  the  system.  One  measure  of  this  degradaton  is  the  deviation  of  the 
measured  beamwidth  from  theoretical.  Figure  8a  illustrates  the  effective 
beam  response  of  a  towed  array  system  to  the  noise  from  a  source  on  the 
towship.  The  arrival  angle  was  about  6  deg.  from  forward  endfire.  The 
results  for  the  three  statistics  indicate  the  wiggling  of  the  array  over 
the  13  minute  time  period  required  to  collect  the  data  has  had  no  adverse 
effects.  All  three  statistics  produce  beam  responses  near  theoretical. 
This  is  not  the  case  near  broadside  where  the  effects  of  the  wiggling  are 
maximum.  The  bottom  reflected  arrival  for  the  same  time  period  was  received 
70  deg.  from  forward  endfire.  The  results  in  Fig.  8b  indicate  that  the 
effective  beam  response  obtained  by  all  three  statistics  has  increased.  The 
theoretical  half-power  beamwidth  should  be  about  2.15  beam  spacings.  It  is 
2.4  for  the  median,  2.66  for  the  dB  average  and  2.95  for  the  power  average. 
Hence,  the  median  and  the  dB  average  have,  in  effect,  “ straightened”  the 
array  compared  to  the  results  for  power  averaging.  The  median,  however,  is 
the  best  statistic  for  this  case. 


CONCLUSION 


It  has  been  shown  by  example  that  the  power  average  is  not  always  the  best 
estimator  for  underwater  acoustic  data.  The  median  and  the  dB  average  are 
very  effective  in  eliminating  or  discriminating  against  some  of  the  sources 
of  data  degradation  in  the  environment  and  in  the  data  measurement  and 
processing  systems  and  give  improved  result  over  the  power  average. 
Tradition  is  probably  the  main  reason  the  power  average  has  been  used  so 
extensively  by  the  underwater  acousticians.  If  the  power  average  is  not 
vital  to  the  processing  technique  the  dB  average,  median,  or  some  other 
statistic  may  provide  significant  improvement  in  the  results.  This  has 
been  the  experience  of  the  authors  in  the  measurement  and  analysis  of 
towed  array  data.  The  technique  is  sufficiently  general  that  it  could  be 
applied  to  other  situations,  such  as  detection  or  signal  processing,  with 
similar  results. 


SACIANTCEN  CP* 32 


17-14 


WAGSTAFF  &  BERROU:  Power  averaging  the  best  estimator.' 


(  3  )  DIRECT  arrival  from  tow  SHIP  H*) 


BEAM  NUMBER 


(b) 


BOTTOM  REFLECTED  NOISE  FROM  TOWSHIP  HO*) 


« 


2 


20  21  22  2J  24  2$ 

BEAM  NUMBER 


PIG.  8  BEAM  RESPONSE  TO  A  SOURCE  (a)  6  DEG.  FROM  FORMARD  ENDFIRE  AMD 
(b)  70  DEG.  FROM  FORNARD  ENDFIRE 


SACLANTCEN  CP-32 


17-15 


WAGSTAFF  &  BERROU:  Power  averaging  the  best  estimator? 


REFERENCES 


1.  Berrou,  J.L.  and  Wagstaff,  R.A.  “Virtual  beams  from  an  FFT  beamformer 
and  their  use  to  assess  the  quality  of  a  towed-array  system", 
Proceedings  of  the  IEEE  ICASSP  82  conference  in  Paris,  May  82. 

2.  Hamson,  R.M.  and  Wagstaff,  R.A.  "RANDI  II  -  An  ambient  noise  model 
that  includes  coherent  hydrophone  summation  for  sonar  systems  in 
shallow,  intermediate,  and  deep  water".  Presented  at  and  published  in 
the  proceedings  of  the  Ambient  Noise  Conference  at  SACLANTCEN,  May  82. 

3.  Wagstaff,  R.A.  “Iterative  techniques  for  ambient-noise  horizontal 
directionality  estimation  from  towed  line-array  data",  J.  Acoust.  Soc. 
Am.  63,  1978,  863-869. 


SACLANTCEN  CP-32  17-16 


FRASER:  Properties  of  low-frequency 
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ABSTRACT 

Models  for  ship-generated  ambient  noise  vary  considerably  in  the  detail 
required  of  their  input  data  on  sources  and  the  environment.  This  paper  describes  a 
model  that  predicts  ambient-noise  levels  and  cumulative  distribution  functions  for 
Poisson-distributed  shipping.  Representing  statisticaiiy  both  the  locations  end  source 
levels  of  the  ships  serves  to  minimize  the  input  requirements.  It  then  proves 
relatively  easy  to  study  parametrically  the  sensitivity  of  the  predictions  to 
propagation  conditions,  shipping  densities,  source-level  distributions,  and  array 
directivity.  The  model  is  also  well-matched  to  the  limited  environmental  information 
often  available  for  sea  trials.  Predictions  will  be  compared  with  experimental  data  to 
illustrate  the  utility  of  the  model. 


1  INTRODUCTION 

Many  acoustic  sea  trials  collect  ambient-noise  data  In  the  frequency  range 
of  10  to  260  Hz  where  ship-radiated  noise  Is  the  dominant  source  mechanism. 
However,  a  very  considerable  effort  is  then  required  if  enough  information  on  the  ship 
sources  and  on  environmental  effects  is  to  be  collected  ;o  permit  the  a.nbiant-nc'ss 
field  to  be  modelled  deterministically.  Time  and  resources  are  very  seldom  available 
to  do  this.  Then  we  must  rely  on  statistical  estimates  of  source  levels  and  locations 
of  ships,  and  upon  transmission-loss  model  predictions,  to  provide  the  necessary  input 
data  for  the  ambient-noise  model. 

It  is  essential  to  'ava  a  thorough  understanding  of  the  limitations  of  the 
predictive  capability  of  such  models.  How  dependent  is  each  aspect  of  the  statistical 
pradiction  upon  each  of  the  input  parameters?  How  much  confidence  can  be  placed  In 
the  model  if  environmental  or  shipping  conditions  change?  Are  all  of  the  major 
influences  en  ambient-noise  statistics  accounted  for  correctly  in  the  model?  it  is 
hoped  that  this  paper  and  the  ambient-noise  model  FLAN  described  herein  will  make 
significant  progress  towards  answering  these  questions. 

The  model  la  similar  to  a  few  others’*3  in  that,  it  divide*  the  ocean  into 
arbitrary  areas  within  which  transmission  loss  to  a  sensor  is  opnstant.  The  number  pf 
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ships  in  each  of  the  areas  is  assumed  to  be  a  random  variable  (abbreviated  rv) 
governed  by  Poisson  statistics.  Based  on  these  assumptions  and  statistically 
distributed  source  levels  (SL),  we  can  estimate  the  zeroth  and  first  order  ambient 
noise  statistics,  i.e.  mean,  standard  deviation  (sd),  and  (cumulative)  distribution 
function  (cdf),  where  the  cdf  F(x)  is  related  to  the  probability  density  function  (pdf) 
f(y)  by  the  following  expression: 


F(x)  = 


f(y)  dy  . 


FLAN  differs  from  the  other  models1"3  to  some  extent  in  its  ease  of 
specifying  input  parameter  values.  It  also  uses  a  different  procedure  for  specifying 
the  areas  governed  by  each  independent  Poisson  distribution,  and  for  convolving  the 
contribution  of  each  area  to  the  total  noise  field. 

A  set  of  predictions  is  made  in  Section  3  to  illustrate  the  influence  of  each 
adjustable  parameter  upon  the  noise  statistics.  To  indicate  that  the  model,  using 
plausible  input  values,  is  able  to  reproduce  experimental  data,  some  of  the  examples 
are  compared  with  data  collected  by  DREA  during  sea  trials.  Section  4  briefly 
describes  the  the  influence  of  horizontal  directivity  of  an  array  of  sensors  upon  the 
ambient-noise  statistics.  A  summary  of  points  covered  and  the  main  conclusions  that 
stem  from  this  study  appear  in  Section  5. 


2  THE  MODEL 


2.1  ASSUMPTIONS 

In  the  FLAN  model  for  ship-generated  ambient  noise,  the  emphasis  Is  on 
making  the  input  specification  as  simple  as  possible.  In  deference  to  tills  aim  several 
assumptions  are  made.  These  assumptions  are  more  easily  visualized  if  reference  is 
made  to  the  map  in  Fig.  1,  on  which  are  defined  some  of  the  model  parameters.  The 
principal  assumptions  are  as  follows: 

(1)  the  ocean  area  of  interest  can  be  divided  into  a  number  of  subareas, 

each  with  constant  mean  transmission  loss  to  the  receiving 
hydrophone; 

(2)  transmission  loss  (TL),  or  more  specifically,  transmission  efficiency 

(TE)  defined  by  TE  »  1O’^TL/10\  is  independent  of  azimuth  within  a 
specified  sector  (In  this  paper  360°  will  be  assumed  for  simplicity); 

(3)  points  (1)  and  (2)  above  suggest  that  the  subareas  should  be  annular 

In  shape,  as  illustrated  In  Fig.  1 ,  with  the  sensor  or  array  located  at 
their  center; 

(4)  the  density  of  shipping  (number  per  unit  area  Q)  Is  constant  between 

arbitrary  minimum  and  maximum  range  limits,  Rm|n  and  R^, 
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(5)  the  number  of  ships  in  a  given  annulus  is  an  rv  nh  drawn  from  a  Poisson 

distribution,  with  a  mean  ship  count  of  bh  In  the  h-th  annulus; 

(6)  each  contributing  class  of  ships  is  distributed  uniformly  throughout  the 

area  of  interest  so  that  the  composite  source-level  distribution  will 
prevail  everywhere; 

(7)  transmission  loss  will  normally  be  assumed  to  increase  as  20  log  R  up 

to  transition  range  R0  and  as  1 0  log  R  thereafter,  with  an  additional 

attenuation  proportional  to  range  and  specified  by  a  coefficient  or 
(there  is  also  provision  for  using  TL  calculated  with  a  ray-trace 
model); 

(8)  virtually  any  source-level  distribution  can  be  mocked  up  by  combining 

m  relatively  narrow  distributions  of  arbitrary  mean  source  level  and 
weighting  factor; 

(9)  the  statistical  predictions  are  representative  of  an  ensemble  average 

over  a  large  number  of  statistically  independent  data  samples. 


2.2  MODEL  DEVELOPMENT 

The  FLAN  model  is  called  upon  principally  to  relate  ship  and  environmental 
information  to  cumulative  distribution  functions  of  ambient-noise  level.  Therefore 
convolution  of  the  density  functions  of  Input  parameters  will  be  expected  to  play  an 

important  role  in  the  calculation.  The  USI  model1  uses  Ingenious  techniques  to  avoid 
traditional  convolution  formulations.  However,  in  the  present  case  computational 
efficiency  is  not  of  primary  concern.  Rather,  the  emphasis  is  on  minimizing  the  number 
of  user-input  parameters  and  easing  the  burden  of  specifying  their  values. 

The  discussion  will  first  focus  on  the  model  development.  Numerical 
procedures  will  be  described  in  Sectic  2.3. 

The  method  of  dividing  the  ocean  intu  annuli,  each  with  constant 
transmission  efficiency,  was  illustrated  in  Fig.  1.  The  model  development  proceeds 
from  this  point  as  shown  schematically  In  Fig.  2.  Four  factors  must  be  taken  Into 
account:  the  sources,  the  ocean  medium  for  conducting  the  sound,  the  hydrophone 
sensor  or  array,  and  the  data-analysls  procedures. 

The  source  power  for  each  ship  is  an  rv  y  drawn  from  the  same  density 
f(y),  as  Implied  by  Fig.  2.  The  selection  of  a  form  for  the  source-power  pdf  is 
governed  by  two  requirements:  (1)  that  Its  characteristic  funotlon  be  easily 
evaluated,  and  (2)  that  it  be  easily  transformed  into  the  log  domain  to  facilitate 
comparison  with  published  source-level  distributions.  Perhaps  the  simplest  pdf  to 
fulfill  these  conditions  is  the  gamma  density.  In  order  to  fit  the  great  variety  of 
possible  8L  density-function  widths  and  shapes,  it  is  necessary  to  provide  for  a 
weighted  sum  of  gamma  densities  with  different  median  values.  The  souroe-power  pdf 
then  has  the  following  form: 


f(y).Wmill  yFe’cy, 
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where  in  the  present  use  p  is  arbitrarily  chosen  as  2,  Wm  is  the  weight  for  the  m-th  SL 
density,  and 


„  _  o  exp(-SLom/g) 

c“3  te;  * 

Here  SLom  is  the  most  probable  value  of  the  m-th  source-level  density  and 
g  =  10  log  e  =  4.343. 

As  noted  previously,  the  number  of  ships  in  each  annulus,  nh,  is  also  an  rv 
(Poisson-distributed)  with  parameter  bh,  the  mean  number  of  ships  in  the  annulus. 
Since  each  annulus  has  the  same  width  Ar,  bh  is  proportional  to  the  radius  rh.  The 
medium  modifies  each  noise  source  by  the  factor  TEh,  so  that  the  total  noise  power 
from  annulus  h  becomes 


xh 


TEh  y*  • 


The  sensor  provides  to  the  processor  an  Incoherent  noise  sum  over  all 

annuli: 


z 


Xh- 


Rather  than  convolving  the  pdfs  of  the  rvs  xh  to  determine  that  of  z,  It  Is 
desirable  to  obtain  and  multiply  together  the  characteristic  functions  (cf)s  for  the  rvs 
In  order  to  determine  the  ambient-noise  power  distribution.  The  process  of  summing 
the  noise  contributions  from  the  various  annuli  Is  a  "generalized  Poisson  process"  as 
described  by  Papouils  (Ret.  4,  pago  676).  Tho  cf  for  x  Is  then  given  by: 


invoking  the  following  property  of  cfs 

<tx  «  $y(0)d) , 


where  x  ■  yd.  The  cf  for  z  is 


H 

*i<“>  •  n  ’ 
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and  for  y  it  is 


y  ^  [1  ♦  (C0TEh/c)2f 


2.3  NUMERICAL  PROCEDURES 

The  Inversion  of  $z(60)  to  form  f(z)  is  accomplished  by  a  discrete  Fourier 

transform  (DFT)  using  a  procedure  devised  by  Bird6  and  implemented  at  DREA  by 

Walker6.  Bird  has  shown  that  the  DFT  can  be  readily  performed  provided  that  the 
input  densities  are  zero  above  an  arbitrary  limit  Y.  He  then  deduces  appropriate 
values  for  Y,  and  for  J,  the  number  of  terms  in  the  finite  approximation  to  the  DFT  for 
the  general  case. 

Walker  obtains6 


)  slnc(- 


irlz 


)  e 


-jrr/z/Y 


(D 


for  z  <  Y,  where  slnc(z)  *  sln(z)/z.  (It  is  assumed  that  f(z)  *  0  if  z  >  Y.)  In  Eqn.  (1) 

$z(-y-)  *  $z(o» . 

In  the  case  of  a  Poisson  density  it  Is  necessary  to  truncate  the  probability  density 
function  (pdf)  at  a  value  Y  to  fulfill  the  above-mentioned  condition.  (In  FLAN,  Y  was 
chosen  to  be  equal  to  the  mean  noise  level  plus  1 5  dB.)  Furthermore,  the  infinite  sum 
over  l  in  Eqn.  (1)  can  not  be  accomplished  In  practice.  Bird  provides  arguments 
suggesting  that  a  practical  rule  for  truncation  is  to  determine  J,  the  maximum  value  of 
/,  from  the  following  expression: 


cggauio-. 

21tl 


The  expression  used  to  evaluate  F(z)  in  FLAN  is  Walker's  Eqn.  (A  10): 
J 


Hz) 
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3  ANALYSIS  PROCEDURES  AND  PREDICTIONS 


3.1  INPUT  PARAMETERS 

The  program  FLAN  was  run  for  a  variety  of  plsusible  input-value 
combinations  to  test  their  influence  upon  the  cdf,  mean  and  standard  deviation.  The 
input  variables  are  listed  in  Table  1. 

TABLE  1.  Input  parameters  for  noise  model  FLAN 

Q  Density  of  shipping  (number/km2) 

Rml(>  Minimum  analysis  range  (km) 

Rmax  Maximum  analysis  range  (km) 

N  Mean  total  number  of  ships  (If  specified,  governs  maximum 

analysis  range) 

Transmission-loss  range  dependence  (generated  externally 
by  an  acoustical  model  or  defined  by  the  following 
parameters: 

Rc  -  crossover  range  (from  spherical  to  cylindrical 

spreading  loss) 

or  -  attenuation  coefficient  (dB/km) 

A r  Width  of  annulus  (usually  about  2  km) 

Source-level  data: 

M  -  number  of  density  functions 

SL^  -  most  probable  source  level  of  each  density  function 

Wm  -  weighting  factor  for  each  density  function. 


3.2  PREDICTIONS— DEEP  WATER 

First  we  shall  examine  the  dependence  of  the  statistics  upon  shipping 
density,  with  other  factors  held  constant  at  typical  values,  and  transmission  loss 
assumed  to  increase  linearly  with  range.  Four  cdfs  covering  the  range  of  light  to 
heavy  shipping  densities  are  shown  in  Fig.  3.  Several  points  can  be  made  with 
respect  to  these  cdfs: 

(1)  the  median  noise  level  is  approximately  proportional  to  10  log(shlpping 

density):  as  the  density  of  shipping  rises,  so  does  the  mean  ambient 
noise  power; 

(2)  the  slopes  of  the  cdfs  (related  to  the  sd)  Increase  with  decreasing 

shipping  density  Q:  fewer  ships  provide  less  chance  of  averaging 
out  the  differences  in  source  level,  transmission  loss  and  ship 
numbers; 
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(3)  the  cdf  is  nearly  log  normal  over  2  to  3  standard  deviations,  at  least 
for  the  higher  values  of  Q.  In  a  small  fraction  of  the  samples  and  at 
short  ranges  a  very  few  ships  may  be  exposed  to  significantly 
differ  3nt  propagation  conditions,  such  as  the  nearest  convergence 
zone,  or  the  direct-path  region  at  a  range  less  than  R0  .  This  causes 
the  curves  to  break  near  the  90th  percentile.  Since  the  source 
levels  for  all  ships  are  drawn  from  the  same  pdf,  all  the  cdfs  of 
imbient  noise  will  tend  towards  the  same  value  at  the  high-noise 
limit. 

Several  examples  showing  this  break  in  the  cdf  are  provided  in  Fig.  4  for 
two  types  of  propagation:  on  the  left  for  convergence  zones,  and  on  the  right  for  a 
smoothly  increasing  dependence  of  TL  upon  range.  The  variation  of  cdf  shape  with 
changing  Rmln  points  to  nearby  ships  as  the  source  of  the  curvature. 

In  deep  water  the  mean  and  standard  deviation  of  the  noise  field  appear 
to  depend  most  strongly  upon  the  mean  number  of  ships  N.  This  dependence  Is 
examined  with  the  aid  of  FLAN  predictions  In  Fig.  5.  The  mean  ship  population 
contributing  to  the  noise  field  Is  governed  principally  by  Q,  the  shipping  density,  and 
by  the  upper  range  limit  Rmax.  All  other  model  parameters  are  held  constant  at 
representative  values.  It  is  clear  that  the  major  influence  on  sd  is  the  value  of  N.  For 

example,  a  change  in  ship  density  Q  from  1.0  to  0.2  ships/ 104  km2  is  seen  to  change 
sd  by  only  0.26  dB,  provided  that  a  corresponding  change  is  made  in  R^  to  hold  N 
constant.  It  is  Interesting  to  note  that  the  ambient-noise  sd  curves  In  Fig.  6  follow  a 
trend  similar  to  the  sd  of  total  ship  count,  where  the  latter  is  assumed  to  obey  Poisson 
statistics. 


A  major  concern  in  ambient-noise  studies  is  the  difficulty  of  specifying 
source-level  distributions  with  any  degree  of  confidence.  The  Influence  of  several 
plausible  source-level  pdfs  upon  the  sd  of  ambient-noise  level  is  displayed  In  Fig.  6. 
The  curves  relating  ambient  noise  sd  to  source-level  sd  were  calculated  for  a  single 
source-level  pdf  shape,  namely  that  of  Case  2  In  the  upper  left  corner  of  the  diagram. 
Larger  standard  deviations  are  obtained  by  stretching  the  pdf  to  cover  a  larger  range 
of  source  levels.  The  three  source-level  pdfs  illustrated  (Cases  1  to  3)  were  chosen 
to  assess  the  effect  of  extreme  changes  In  the  density  function  shape  upon  the  sd  of 

ambient  noise.  The  results  for  a  Q  of  1  ship/104  km2  and  a  source-level  sd  of  8  dB 
are  indicated  by  crosses.  In  all  cases  tested,  quite  drastic  changes  in  shape  or  in  sd 
of  the  source-level  pdfs  were  necessary  in  order  to  produce  even  modest  changes  in 
the  ambient-noise  sd. 

Next  we  present  a  cursory  examination  of  the  effect  of  transmission-loss 
trends  upon  the  statistics  of  interest.  Several  plausible  deep-water  range 
dependencies  of  TL  are  illustrated  in  Fig.  7.  Two  curves,  labeled  Labrador  Sea  and 
Mld-Atlanttc  were  generated  by  a  ray-trace  program  for  specific  ocean  areas.  The 
other  three  curves  reflect  various  combinations  of  R0  and  a.  Predictions  of  cdf  shape 
have  been  derived  for  tho  four  lowest  curves  in  Fig.  7  and  are  displayed  in  Fig.  8. 
(The  Labrador  Sea  case  will  be  deferred  to  Fig.  0.)  The  variations  in  shape  and  slope 
of  the  cdfs  are  seen  to  be  quite  weak,  considering  the  large  variations  In  R0  and  in  or. 
When  ot  Is  greatest,  the  influence  of  the  many  distant  ships  Is  diminished  relative  to 
those  nearby.  This  Increases  the  ambient-noise  variability,  i.e.  steepens  the  cdf 
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slope.  The  effect  is  evident  in  Fig.  8.  Little  dependence  of  the  cdf  on  R0  is  observed, 
however,  except  above  the  90th  percentile.  For  the  convergence-zone  propagation 
(mid-Atlantic  TL  profile),  both  the  slope  of  TL  and  of  the  cdf  are  most  closely  aligned 
with  the  geometric  cases  for  which  a  =  0.01  dB/km. 

Referring  again  to  Fig.  7,  we  observe  a  large  offset  between  TL  in  the 
Labrador  Sea  and  that  typical  of  the  mid-Atlantic.  If  all  other  factors  are  assumed 
equal,  then  we  would  expect  the  mean  ambient-noise  level  to  be  many  decibels  higher 
in  the  Labrador  Sea  environment.  The  cdf  predictions  for  the  two  locations  are 

represented  by  the  top  and  bottom  curves  for  equal  densities  of  1  ship/104  km2  in 
Fig.  9.  The  sd  of  1 .55  dB  for  the  lowest,  curve  is  quite  typical  of  measurements  near 

Bermuda7.  However,  DREA  data  (represented  by  dots  in  Fig.  9)  collected  in  the 
Labrador  Sea  imply  a  much  lower  ambient-noise  level  and  steeper  cdf  slope  than  that 
represented  by  the  top  curve.  The  good  fit  provided  by  the  middle  curve  was 
obtained  simply  by  reducing  the  shipping  density  by  a  factor  of  5,  and  the  ocean- 
basin  size  by  30%.  Both  changes  are  quite  plausible.  How  much  the  source-level 
densities  differ  in  practice,  and  how  significant  are  different  shipping  and  fishing 
patterns  is,  of  course,  much  more  difficult  to  assess. 


3.3  PREDICTIONS— SHALLOW  WATER 

The  shallow-water  acoustical  environment  requires  some  changes  in  scale 
to  be  made  to  Input  parameters  of  the  FLAW  model.  Transmission  loss  at  short  ranges 
tends  to  be  less  than  In  deep  water,  but  often  is  subject  to  a  relatively  large 
attenuation  as  range  Increases.  These  differences  in  TL  can  be  accommodated  by 
reducing  R0  and  increasing  a,  relative  to  their  deep-water  values.  (DREA  has  been 
cataloging  much  of  Canada's  eastern  shallow-water  area  in  terms  of  those 
parameters.) 

It  is  also  anticipated  that  should  be  much  smaller  in  the  confined 
continental  sholf  areas  than  In  the  deep  ocoan  basins.  A  value  of  RttUX  *  200  km  was 
chosen  for  this  study. 

Source  levels  of  fishing  vessels  generally  are  somewhat  smaller  titan  those 
for  merchant  ships.  However,  during  fishing  operations  thoir  radlatod  noiso  levels  can 
be  extremely  high  and  quite  variable. 

The  preceding  factors  all  tend  to  increase  the  slope  of  cdfs  relative  to 
their  typical  deep-water  values.  A  comparison  of  two  cases  for  which  data  are 
available  is  made  In  Fig.  10.  The  shullow-wntor  data  were  obtained  by  DREA  in  on 
area  characterized  by  heavy  shipping;  the  deep  water  data  were  obtained  off 
Canada's  east  coast.  Historical  shipping  densities  wore  used  in  the  model  predictions, 
togethor  with  the  source  level  density  function  illustrated  as  Cose  2  In  Fig.  6.  Rm„ 
took  on  values  of  750  km  and  200  km  for,  respectively,  the  deep  and  shallow-water 
eases.  The  only  parameter  adjusted  to  provide  the  model  fits  shown  in  Fig.  10  wore 
the  mean  values  of  the  two  source-level  pdfs.  The  main  factor  contributing  to  the 
steepness  of  the  shallow-water  cdf  Is  the  elimination  of  contributions  normally  made 
by  the  many  ships  beyond  200  km  range  in  deep-water  areas. 
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The  effect  of  variations  in  attenuation  coefficient,  a,  upon  the  cdf  is 
displayed  in  Fig.  1 1 .  Since  a  at  low  frequencies  can  be  an  order  of  magnitude  greater 
in  shallow  water  than  in  deep,  its  effect  upon  the  cdf  can  be  significant.  As  in  the 
deep-water  case  (Fig.  8),  the  cdf  becomes  steeper  as  a  increases.  A  low-noise  limit 
must  be  reached,  however,  as  the  shipping  noise  drops  to  zero,  leaving  only  wind¬ 
generated  noise  in  its  stead. 


4  ARRAY  DIRECTIVITY 

The  directional  sensitivity  (directivity)  of  an  array  of  sensors  can 
significantly  influence  the  iow-frequency  ambient-noise  statistics.  For  example,  a 
narrow  horizontal  beam  will  reduce  the  number  of  contributing  ships  from  an  isotropic 

shipping  distribution  by  the  ratio  of  beamwldth  to  360°  (or  1 80°  in  the  case  of  a  line 
array).  Not  only  is  the  mean  noise  power  reduced  by  this  factor,  but  the  cdf  slope, 
and  consequently  the  sd,  are  increased  by  an  amount  that  can  be  assessed  by 
reference  to  Figs.  3  and  5,  which  show  the  influence  of  mean  ship  count  upon  cdf 
shape  and  upon  sd,  respectively.  The  effect  of  sideiobe  structure  requires  special 
attention,  and  wilt  not  be  addressed  in  this  paper. 


5  SUMMARY  AND  CONCLUSIONS 

A  parametric  amblent-nolse  model  FLAN  has  been  developed  to  assess  the 
impact  of  variations  in  ship  source  levels  and  numbers,  of  transmission  loss,  and  of 
sensor  directivity  upon  the  first-order  statistics  of  ambient  noise  at  low  frequencies, 
A  special  effort  has  been  made  to  minimize  the  number  of  input  parameters  and  the 
effort  required  to  specify  their  values,  possibly  at  the  expense  of  generality.  This 
was  done  (1)  to  match  the  model's  Input  requirements  to  the  many  experimental 
situations  where  little  is  known  about  the  ships  contributing  to  the  noise  field  and  the 
prevailing  sound-transmission  conditions,  and  (2)  to  make  it  easy  to  vary  each 
parameter  in  a  systematic  manner  to  assess  its  impact  upon  the  statistics  of  ambient 
noise. 


The  principal  conclusions  reached  In  the  study  are: 

(1)  a  wide  variety  of  experimental  data  can  be  fit  by  such  a  model, 

(2)  good  model  fits  to  experimental  cdfs  do  not  Imply  an  unambiguous  set 

of  Input-parameter  values, 

(3)  under  many  conditions  the  cdf  of  ambient  noise  is  predicted  to  be  log 

normal  over  two  to  four  standard  deviations, 

(4)  the  ad  of  ambient  noise  Is  principally  dependent  upon  the  number  of 

ships  contributing  to  the  noise  field, 

(6)  in  the  case  of  directional  arrays  the  number  of  ships  contributing 
strongly  to  the  noise  field  is  diminished,  so  that  the  sd  tends  to 
become  larger. 
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(6)  the  sd  of  ambient  noise  is  weakly  dependent  upon  the  source-level 

pdf  of  contributing  ships,  upon  RmJn,  and  upon  R0, 

(7)  the  shape  of  the  cdf  (as  opposed  to  its  slope)  draws  away  from  log¬ 

normal  behavior  as  Rm!n  and  Q  decrease,  and  as  a  increases, 

(8)  shallow-water  cdfs  tend  to  be  much  steeper  in  slope  than  those  for 

deep  water,  primarily  because  of  the  smaller  region  over  which  ships 

can  contribute  to  ambient  noise. 

Extensions  of  the  model  to  include  explicitly  the  effect  of  shipping  lanes, 
convoys,  azimuthal  dependence  of  transmission  toss,  and  TL  fluctuations  can  be 
envisaged,  in  order  to  increase  its  ability  to  cope  with  realistic  environments.  The 
cost  of  such  improvements  will  be  additional  complexity  in  the  input  specification. 
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FRASER:  Properties  of  low-frequency 
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DETECTION  MODELS  AMO  TARGET  INFORMATION  PROCESSING 

by 

Dr.  James  W.  Bond 
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San  Diego,  Ca.,  USA 


ABSTRACT 


The  passive  sonar  equation  leads  in  a  natural  manner  to  sensor  detection 
performance  contours.  Useful  information  can  be  extracted  from  observed 
fractional  holding  time  on  non-cooperative  targets  by  use  of  these  contours. 
The  derived  information  could  aid  in  multiple  sensor  contact  data  associa¬ 
tion  and  the  localization  and  tracking  of  poorly  held  targets.  Target 
source  level  and  sensor  processor  recognition  differential  need  not  be  input 
to  the  prediction  model  to  prioritize  candidate  contact  data  associations. 

A  surprising  example  shows  that  use  of  predictions  to  estimate  the  range  to 
initially  detected  targets  is  not  likely  to  work.  An  idealized  but  reason¬ 
able  example  shows  that  the  probability  of  a  target  remaining  undetected  at 
further  ranges  and  first  detected  at  a  given  range  can  be  Independent  of 
the  range  of  the  target  from  the  sensor.  The  theory  of  stochastic  differ¬ 
ential  equations  applied  to  tracking  a  target  may  alleviate  the  difficulty. 
The  paper  closes  with  recommendations  for  the  development  of  model  capabili¬ 
ties  to  provide  more  useful  outputs  for  Information  processing. 


INTRODUCTION 

Detection  models  for  passive  acoustic  sensors  have  been  used  for  operational 
research,  sensor  trade-off  studies,  and  planning  sensor  deployments.  Detec¬ 
tion  models  have  been  little  used  fo*  real-time  Information  processing.  The 
purpose  of  this  paper  is  to  describe  several  applications  of  models  to  real¬ 
time  information  processing.  The  use  of  detection  models  to  aid  in  contact 
data  association  and  target  localization  and  tracking  will  be  nv*  ted. 

To  provide  a  context  for  the  discussion,  postulate  a  processing  ^enu 
designed  to  handle  non-cooperative  targets  through  the  reduction  and  analy¬ 
sis  of  passive  acoustic  sensor  data.  The  target  data  consists  of  detections 
of  signals,  whose  frequency,  bearing  relative  to  the  sensor,  and 
signal  level  relative  to  ambient  noise  are  measured.  The  environmental  date 
consists  of  directional  ambient  noise  and  clutter  levels  estimated  within 
the  processing  center  as  well  as  externally  provided  data  to  the  center. 

Targets  held  consistently  by  several  sensors  pose  no  intrinsic  information 
handling  problems  while  targets  held  infrequently  pose  very  difficult 
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handling  problems.  For  infrequently  held  targets,  detection  models  appear 
useful  in  two  ways:  (1)  to  extract  meaning  from  detections  and  non-detections 
and  (2)  to  supplement  measurement  data  on  a  target  with  predicted  data. 
Application:,  using  models  in  each  of  these  ways  will  be  discussed. 

In  section  one  an  operationally  observable  measure,  fractional  holding  time, 
is  defined.  Previous  results  obtained  by  the  author  comparing  opera¬ 
tional  fractional  holding  time  with  predicted  seasonal  detection  performance 
contours  are  briefly  reviewed.  These  operations  analysis  results  provide,  in 
particular,  a  characterization  of  the  variability  of  fractional  holding  time. 
The  variability  of  the  observable  is  especially  important  because  it  provides 
the  basis  for  an  assessment  of  the  operational  usefulness  of  detection  per¬ 
formance  predictions. 

In  section  two,  a  general  approach  to  using  detection  models  to  aid  in  con¬ 
tact  data  association  is  identified.  The  model  is  used  to  refine  geographic 
consistency  algorithms  now  commonly  used  for  data  association  by  accounting 
for  the  probabilities  of  observing  given  fractional  holding  times  condi¬ 
tioned  on  the  association.  The  predictions  required  need  not  use  an  input 
value  of  target  source  level,  thus  avoiding  the  often  very  difficult  problem 
of  estimating  this  parameter. 

In  section  three  the  use  of  detection  models  for  single  sensor  localization 
is  first  discussed.  The  inference  of  target  range  from  fractional  holding 
time  is  investigated.  An  example  indicates  that  a  most  likely  range  for  an 
initial  contact  need  not  exist.  A  more  general  approach  to  target  tracking 
appears  to  be  required.  The  use  of  stochastic  differential  equations  seems 
to  offer  some  chance  of  overcaning  these  difficulties  with  the  added  advan¬ 
tage  that  use  can  be  made  of  sensor  non-holding  information  (so-called 
negative  information)  as  well  as  better  use  of  sensor  holding  information 
(positive  information). 

Section  four  contains  a  brief  discussion  of  the  directions  that  passive 
acoustic  sonar  modeling  might  take  to  better  aid  in  information  processing. 

It  is  suggested  that  the  predictions  should  more  fully  utilize  measurement 
data  available  in  the  processing  center  and  should  aim  to  predict  additional 
observable  features  of  detected  signals. 


1  SINGLE  SENSOR  DETECTION  PREDICTIONS 

The  basis  for  passive  acoustic  sensor  detection  performance  predictions  is 
the  passive  sonar  equation  (in  decibels) 

SEt  «  SL  -  RD  -  PI  +  et,  t  a  0  .  (1) 

From  the  viewpoint  of  the  user  of  the  predictions,  source  level  (SL)  summar¬ 
izes  the  relevant  target  characteristics,  while  recognition  differential  (RD) 
summarizes  the  relevant  signal  processor  characteristics.  The  performance 
index  (PI)  is  a  function  of  range,  seasonal,  and  directional  dependent  trans¬ 
mission  loss  (TL),  season  and  directional  dependent  ambient  noise  (AN),  and 
directional  dependent  array  gain  (AG)  in  the  following  manner: 
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PI  =  TL  +  AN  -  AG  . 

The  performance  index  therefore  summarizes  all  of  the  non-target  and  non¬ 
processor  relevant  properties  of  the  environment  and  sensor  whose  performance 
is  being  predicted.  Signal  excess  SE^  varies  as  a  function  of  time  t. 

One  possible  representation  of  the  stochastic  behavior  of  SEj-  is  the  A-o 
jump  model  which  is  a  simple  step  function  of  time.  The  family  of  step  func¬ 
tions  are  conditioned  by  assuming 

(1)  ej.  is  normally  distributed  with  zero  mean  and  variance  o2 

(2)  cs  and  et  are  different  with  probability  1  -  exp[A(s-t)]  and 
equal  with  probability  exp[-A(s-t)],  A  >  0,  s  >  t. 

These  criteria  define  a  commonly  used  stochastic  process  for  SEt  which  has  a 
normal  distribution  with  variance  o2  and  an  exponential  autocorrelation  func¬ 
tion  with  a  relaxation  time  constant  of  1/A. 

A  natural  manner  to  summarize  expected  performance  over  a  period  of  time  T 
is  to  predict  the  probability  p  that  the  instantaneous  signal  excess  SE^  is 
non-negative.  For  example,  if  the  expected  value  of  signal  excess  were  zero, 
then  SEt  would  be  expected  to  be  non-negative  half  of  the  time  and  p  would 
be  i.  furthermore,  it  is  assumed  that  the  operational  measure  of  holding  per¬ 
formance,  the  fraction  of  time  that  contact  is  held,  fractional  holding  time 
(FHT)  is  an  unbiased  statistical  estimate  of  p.  Values  of  FHT  will  there¬ 
fore  be  the  observables  which  will  be  compared  to  predicted  values  of  p. 

From  (1)  it  follows  that 

p  =  Prob{SEt  >  0}  =  Prob(t t  >  PI  -  SL  +  RD} 


with 


O(x)  =  /xro  exp[-iy2]  dy  . 

(2) 

Another  formulation  of  equation  (2)  is 

$*l(l  -  p)  *  (PI  -  SL  +  RD)/a. 

(3) 

Equation  (3)  exhibits  4>"*(1  -  p)  as  linearly  related  to  PI  with  the  slope 
of  the  line  inversely  proportional  to  the  standard  deviation  a  of  signal 
excess  and  with  intercept  (SL-  RD)/o. 

The  observable  to  be  compared  with  the  predictions  is  fractional  holding 
time,  FHT,  defined  as  a  sum  of  holding  intervals  of  a  particular  target 
acoustic  source  over  the  sum  of  opportunity  time  in  a  region  bounded  by  two 
PI  contours  1$  dB  apart  as  illustrated  in  Figure  1.  The  FHT  for  such  a 
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region  will  be  compared  with  the  probability  p  for  the  average  value  of  the 
two  PI  values.  The  expected  variability  of  the  observable  FHT  is  a  function 
of  the  true  probability  of  non-negative  SEt  in  the  PI  bin  and  the  total  tar¬ 
get  opportunity  time  T  in  the  bin.  Let  f(t)  =  1  if  the  source  is  detected 
and  0  if  a  source  is  not  detected  by  the  sensor.  Then,  from  FHT  = 

1/T  /q  f(t)  dt  and  the  fact  that  e^.,  t  ?  0  is  a  zero  mean  A -a  stationary 
process  it  follows 

E[FHT]  =  p  (4) 

VAR[FHT]  =  2p(l  -  p)[(e"AT-l+AT)/(AT)2]  .  (5) 

The  term  involving  AT  converges  very  rapidly  to  1/AT  as  AT  •>  °°  and  there¬ 
fore  equation  (5)  simplifies  to 

Var[FHT]  -  2p(l  -  p)  1/AT  .  (6) 


The  variation  in  SE  induced  by  target  position  changes  and  target  source 
level  changes  as  well  as  sound  propagation  and  ambient  noise  variation  can 
be  estimated  by  regression  of  FHT  against  PI.  For  this  regression  each 
observation  4>_1(1-FHT)  should  be  weighted  according  to  its  expected  variance. 
The  formulas  used  for  the  results  presented  in  this  paper  are  derived  from 
equation  (6)  by  use  of  the  first  two  terms  of  the  Taylor  expression  of  <I>"1 
evaluated  at  1  -  p  and  resulted  in 


VarrMl-FHT)  = 


(7) 


where  1  -  p  =  <t(x)  and  x  =  E [FHT 3  using  known  properties  of  f. 

Data  was  pooled  for  several  targets  with  no  a  priori  assunptions  as  to  source 
level  or  recognition  differential.  The  standard  deviation  of  signal  excess 
was  also  assumed  to  be  unknown  but  dependent  on  season  but  independent  of 
target.  Figure  2  presents  an  example  of  the  regression  results  obtained. 

Data  is  presented  for  seven  targets  deployed  during  the  winter  season.  The 
left  scale  is  in  terms  of  the  transformed  variable  <1>“1(1-FHT)  and  the  right- 
hand  scale  is  natural.  The  overall  a  results  obtained  are  summarized  in 
Table  1.  The  noteworthy  implication  of  these  results  is  that  summer  esti¬ 
mates  of  n  are  from  4  to  7  db  higher  than  winter  estimates.  This  result 
lends  credibility  to  the  estimates  because  greater  variability  would  be 
expected  during  the  summer  than  the  winter  due  to  propagation  effects. 

The  fact  that  propagation  paths  are  more  sharply  focused  during  the 
summer  would  lead  to  greater  variation  of  received  signal  level  as  a 
function  of  target  position. 

The  variability  in  signal  excess  implied  by  the  values  of  o  in  Table  1 
exceeds  that  estimated  due  to  statistical  variability  by  equation  (7).  Some 
of  this  Increased  variability  is  induced  by  the  use  of  fixed  transmission 
loss  estimates  as  a  function  of  angle  for  each  sensor  and  season  and  the  use 
of  an  average  direction  dependent  ambient  noise  for  each  sensor. 


SACLANTCEN  CP-32 


19-4 


BONO:  Detection  moaeis 


TABLE  1.  EXAMPLE  OF  a  ESTIMATES  OBTAINED  FROM  A  REGRESSION 
ANALYSIS  OF  FHT  ON  REAL  TARGETS. 


Summer  Season  Winter  Season 

SENSOR  1 

★ 

Number  of  Targets 

5 

7 

Degrees  of  Freedom 

26 

41 

Point  Estimate  of  a 

17.1 

10.2 

90%  Fiducial  Bounds  of  o 

12.3 

8.2 

27.9 

13.6 

SENSOR  2 

Number  of  Targets 

4 

7 

Degrees  of  Freedom 

29 

14 

Point  Estimate  of  a 

10.1 

5.9 

90%  Fiducial  Bounds  of  o 

8.1 

4.4 

13.5 

8.7 

♦Case  presented  in  Figure  2. 


Figure  3  shows  an  example  of  the  variability  observed  for  a  nominal  case, 
not  the  best  case  nor  the  worst  case.  The  solid  lines  are  90%  confidence 
bounds  on  the  regression  line  and  the  dotted  lines  are  90%  confidence  bounds 
on  the  observations.  In  this  example,  even  though  the  regression  line  is 
estimated  with  reasonable  confidence,  the  variability  of  the  observations 
is  around  7  dB.  Predictions  which  use  real-time  estimates  of  transmission 
should  reduce  this  loss  or  ambient  noise  variability.  However,  at  this  time 
an  analysis  to  assess  the  variability  of  observed  fractional  holding  time 
relative  to  real-time  predictions  has  not  been  undertaken  to  the  author's 
knowledge. 

The  operational  analysis  results  imply  two  conclusions  relative  to  the  use 
of  models  for  real-time  information  processing: 

(1)  predictions  for  real-time  use  should  be  based  on  real-time  data; 
and 

(2)  although  regression  analysis  can  estimate  o  and  SL-RD  for  a 
large  data  base,,  these  estimates  will  normally  not  be  available 
for  real-time  applications. 


19- b 
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2  USE  OF  PREDICTIONS  FOR  CONTACT  DATA  ASSOCIATION 

The  fundamental  problem  in  data  association  is  to  calculate  the  probability 
that  a  given  set  of  observations  occurs.  The  calculation  can  be  factored 
into  the  calculation  of  two  probabilities: 


(1)  the  probability  that  the  observed  line  characteristics,  e.g., 
frequency,  bandwidth,  stability,  would  be  associated  with  the 
acoustic  sources  of  a  single  target;  and 

(2)  the  probability  that  the  available  operator  heuristic  informa¬ 
tion  would  have  been  generated  if  they  were  on  the  same  target. 

If  detections  are  associated  for  different  sensors,  then  an  additional  factor 
becomes  important: 

(3)  the  probability  of  the  observed  detections  occurring  on  specific 
bearings  for  some  sensors  and  not  occurring  on  specific  bearings 
for  other  sensors  for  each  target  acoustic  source. 

There  appears  to  be  a  real  possibility  of  exploiting  detection  models  to 
improve  the  estimation  of  the  third  factor  and  little  possibility  of  their 
use  in  the  estimation  of  the  other  factors.  Therefore,  we  focus  our  atten¬ 
tion  on  the  two  sensor  contact  association  problem.  It  is  convenient  to 
introduce  the  following  notation  to  aid  in  a  concise  formulation  of  the 
problem: 

The  observations  consist  of  fractional  holding  times  FHTj  (fj)  for  a  speci¬ 
fied  opportunity  time  T,  for  sensors  1  =  1,  2,  M,  and  acoustic  sources 
with  nominal  frequencies  fj,  f2»  ...»  ffl. 

The  sensor  Indices  1  a  1  and  2  are  chosen  so  that  the  contact  association 
consists  of  FHTj  (fj)  and  FHTg  (f^). 

The  FHTj  (fj)  are  all  supposed  to  be  appropriately  calculated,  i.e.,  the 
holding  time  on  a  sensor  i  for  target  source  with  frequency  fj  has  been  on 
a  single  target.  (The  reader  will  note  that  It  would  be  straightforward  to 
extend  the  results  presented  to  the  case  when  these  single  sensor,  single 
source  associations  had  been  made  correctly  with  probability  P. ,  over  the 
opportunity  time  T. ) 

To  Interpret  the  observed  fractional  holding  times  it  is  necessary  to  asso¬ 
ciate  PI  values  to  them.  A  value  of  T  of  24  hours  or  larger  appears  to  be 
reasonable.  To  associate  a  PI  value  to  an  observed  FHT  define  an  opportunity 
region  as  in  Figure  4  compatible  over  time  with  the  bearing  information 
observed  on  the  sources  fj  and  f^  for  sensors  1  and  2,  respectively,  during 
the  opportunity  Interval  dT,  Such  a  region  would  have  lengths  on  the 
order  of  2  VT  where  V  Is  the  average  speed-of-advance  of  the  target.  View 
the  region  defined  in  log  R  coordinates  and  take  the  centroid  (or  a  point 
approximating  It)  as  the  position  for  which  PI  predictions  are  compared 
with  the  observed  FHT. 
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The  PI  predictions  for  the  centroid  should  be  as  accurate  as  possible  and 
therefore  should  be  real-time  predictions  using  sensor  measured  directional 
ambient  noise  and  the  best  available  sound  velocity  depth  profiles  for  trans¬ 
mission  loss  predictions. 

Equation  (3)  exhibits  the  probability  of  detection  p  as  a  function  of  PI, 
with  parameters  o,  RD,  and  SL.  The  value  of  o  can  be  empirically  estimated 
by  sensor  based  on  an  analysis  of  historical  data.  The  value  of  RD  can  be 
updated  in  real  time  from  the  outputs  of  passive  acoustic  sensors  through 
the  use  of  an  automatic  detection  algorithm  to  calibrate  operator  responses. 

A  priori  estimates  of  SL  are  generally  not  available,  so  a  posteriori  esti¬ 
mates  will  be  used.  Toward  this  end,  view  p  as  a  function  of  SL.  Suppose 
that  the  opportunity  time  T  is  long  enough  and  the  target  maneuvers  such 
that  the  average  value  of  SL  observed  by  the  sensors  during  the  oppor¬ 
tunity  time  T  is  independent  of  sensor.  Then  equation  (3)  implies  that 
for  sensors  1,  2,  ...,  M  the  probabilities  pj(SL),  p2 (SL) ,  ...,  p.. (SL ) 
are  derivable  from  the  sensor  PI  values  Plj,  Pip,  ...,  PIM  for  the  cen¬ 
troid  of  the  opportunity  region  as  a  function  of  a  single  unknown  SL. 

A  natural  measure  of  consistency  between  p-j (SL)  and  FHT  for  a  given  acoustic 
source  and  the  M  sensor's  is  the  weighted  mean-square  error  E(SL)  defined 

by 


M 

E (SL )  =  £  (FHT.-p, (SL))2/(Var[FHT.]  +  Var[p. (SL)])  .  (8) 

1=1  1  1  1  1 


The  smaller  E(SL)  can  be  made  the  more  consistent  are  the  observations. 
Furthermore,,  the  value  of  SL  which  minimizes  E(SL)  provides  the  a  posteriori 
estimate  of  SL  as  previously  mentioned. 

The  value  of  Var[FHT]  used  in  equation  (8)  can  be  determined  using  equation 
(5)  with  the  value  of  X  estimated  empirically  from  an  analysis  of  historical 
operational  data.  (Such  an  analysis  is  now  in  progress  at  NOSC  by  the 
author.) 


The  value  of  Var[p]  can  be  estimated  by  differential  approximation  (refer  to 
Figure  5).  It  follows  that 


4P . « (f ) 


(9) 


where  SE  *  -PI  +  SL  -  RD  and  o  is  the  standard  deviation  of  signal  excess 
e-x2/2 

with  <J>(x)  =  by  the  fundamental  theorem  of  calculus.  Assume  the  value 

/2u 

of  Ap  to  match  the  observed  holding  time  as  a  function  of  time  is  uniformly 
distributed.  It  then  follows  that 


Var[p]  =  ^  Ap2 


(10) 


The  desired  approximation  for  Var[p]  is  obtained  by  substitution  of  (9)  into 

(10). 
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The  inclusion  of  the  variances  in  equation  (8)  appears  to  be  important.  It 
helps  ensure  correct  behavior  of  (8)  for  small  p  when  FHT  =  0.  This  is 
because  it  can  be  shown  that 

p/(Var[FHT]  +  Var[p])  -»  constant  as  p  ->  0. 

The  expression  in  (8)  could  clearly  be  used  to  rank  different  associations 
for  the  same  M  sensors.  It  is  not  so  clear  how  to  compare  E(SL)  values 
when  the  sensors  used  for  the  calculations  differ.  An  alternative  approach 
to  the  calculation  of  consistency,  albeit  more  complicated,  would  be  to  cal¬ 
culate  the  probability  of  the  observations  for  an  optimum  choice  of  SL.  For 
this  approach  the  random  variable 

(FHT  -  p(SL))  /  /Var[FHT]  +  Var[p] 

is  assumed  to  be  normally  distributed  with  unit  variance.  Then  the  probabil¬ 
ity  of  observing  Y  satisfying  | Y-p(SL) [  | FHT-p(SL)  1  is  given  by 

l-2'K-|FHT-p(SL)i  /  ^ar[FHT]  +  Var[p])  (11) 

and  the  probability  of  the  observations  FHT^,  ...,  FHT^  is  given  by 
M 

P(SL)  -  £  (l-2‘t,(-|FHT1-p.(SL)|  /  v^TFHf.]  +  Var[p,  (SL)]) )  (12) 

i=l  11  11 

for  a  choice  of  SL  leading  to  the  normal  distributions  of  the  random  vari¬ 
ables.  This  choice  of  SL  is  assumed  to  be  given  by  the  value  which  maxi¬ 
mizes  P(SL). 

For  either  approach,  it  is  suggested  that  SL  be  estimated  iteratively.  A 
good  first  choice  would  be  a  sensor  for  which  the  data  association  was  most 
straightforward,  as  long  as  it  was  not  held  all  the  time.  The  initial  choice 

for  SL  can  then  be  obtained  by  solving  equation  (3)  for  SL  when  the  observed 

fractional  holding  time  is  substituted  for  p. 

The  value  of  E(SL)  or  P(SL)  thus  obtained  would  be  used  in  the  ranking  of 

contact  data  associations.  The  estimates  of  SL  would  of  themselves  be  of 

great  value  since  target  source  level  of  a  non-cooperative  target  is  not 
observable  using  the  usual  non-calibrated  sensors  in  an  incompletely  charac¬ 
terized  environment.  The  method  provides  c  recursive  procedure  for  estimat¬ 
ing  SL  for  all  acoustic  sources  of  all  targets  being  tracked. 

At  the  present  time  this  methodology  is  being  investigated  by  studying  the 
SL  estimates  obtained  from  processing  of  historical  data  associated  to  oper¬ 
ational  targets  with  known  source  level. 
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3  USE  OF  PREDICTIONS  FOR  TARGET  LOCALIZATION  AND  TRACKING 

Another  important  potential  use  of  models  is  to  provide  range  estimates 
when  cross-fixes  are  not  available  as  illustrated  in  Figure  6.  Due  to  uncer¬ 
tainties  in  the  predictions  and  the  expected  variability  in  fractional  hold¬ 
ing  time  it  appears  unlikely  that  range  estimates  obtained  from  comparing 
predicted  probabilities  as  a  function  of  range  with  observed  holding  time 
would  provide  useful  information  when  cross-fixes  are  available  during  the 
time  frame  T  associated  with  the  fractional  holding  time  observation.  It 
is  natural  then  to  consider  the  case  when  an  isolated  detection  has  occurred 
by  a  single  sensor. 

Suppose  that  a  sensor  first  detects  a  target.  Figure  7  sunmarizes  the  per¬ 
formance  of  the  sensor.  Suppose  further  that  the  target  is  moving  toward 
the  sensor  so  that  it  moves  through  regions  of  poorer  coverage  toward 
regions  of  better  coverage.  Then  the  probability  of  being  first  detected  p 
in  region  K  is  the  product  of  the  probability  of  detection  for  region  K  and 
the  probability  of  non-detection  for  the  regions  1,  2,  ....  k-1,  i.e., 

^k  s  Pk(l"Pk_j )(^“Pk»2^  *  *  *  *  (13) 

Several  examples  indicated  that  Pk  tended  to  vary  slowly  with  K.  A  fellow 
worker  at  NOSC,  Stefen  Hui ,  foundKthe  following  important  example: 

Let 


Pk  *  Pj/d-tk-lJpj) 


(14) 


with  p.  *  1/N  to  assure  that  the  pt  as  defined  satisfy  0  ?  p.  •  1 .  Then 

Pi  *  pi  •  h  "  ■  pn- 

To  see  that  this  is  the  case,  proceed  by  mathematical  induction.  P.  *  p,. 
Suppose  Pk  =  Pj.  Then, 


pk+l  =  {Wl"pk)  •** 


fk+1 

Pk 


(1-Pk) 


» 


and  it  suffices  to  verify  that  Pk+i(l"Pk)/Pk  =  1. 
Using  (14)  it  follows  that 

fk+l  l-(k-l)p.  /  p,  \ 

pk  ^'Pk1  ’ '  T-'kp,"  •  l-(k-npj  j  ■ 


(15) 


by  cancellation  after  placing  the  second  factor  under  a  common  denominator. 
Furthermore,  whenever  (15)  holds,  i.e.,  whenever  pk+^  =  pk/(l-pk),  it  follows 


SACLANTCEN  CP-32 


19-9 


BOND:  Detection  models 


that  Pk+,  =  Pk.  For  example,  if  pk  =  1/n,  then  P|c+,  •  i  •  ^ 

The  example  is  not  farfetched  and  suggests  that  in  practice  it  will  be 
difficult  to  estimate  range  from  predicted  probabilities  of  detection  for 
a  sensor.  Indeed,  the  example  shows  that  for  a  target  of  known  track  (radial 
toward  sensor)  and  known  probabilities  of  detection  for  the  sensor  as  a  func¬ 
tion  of  range,  but  unknown  time  of  opportunity  for  the  sensor,  there  need  be 
no  range  information  associated  with  an  initial  detection  event.  To  try  to 
overcome  this  difficulty,  the  past  history  of  the  target  needs  to  be  con¬ 
sidered,  i.e.,  the  use  of  the  model  for  tracking  rather  than  localization 
should  be  considered. 

The  contact  data  association  calculation  suggests  that  SL  estimates  may  not 
be  necessary  for  the  tracking  application.  The  ocean  could  be  partitioned 
into  regions  with  some  boundaries  associated  for  the  sensors  holding  and 
others  defined  from  environmental  considerations.  For  each  region  a  calcu¬ 
lation  of  the  probability  of  observing  the  given  FHT-j  for  sensors  i  =  1,  2, 
....  M  could  be  made  and  the  calculation  of  P (SL )  defined  by  equation  (12) 
made  as  described  in  Section  2.  The  target  track  could  be  assuned  to  have 
random  course/speed  changes.  The  suggested  technical  approach  would  then 
attempt  to  apply  the  theory  of  stochastic  differential  equations  to  the 
problem. 

Stochastic  differential  equations  have  been  successfully  used  for  modeling 
target  motion  and  to  represent  a  more  general  approach  than  Kalman  filtering 
(reference  1).  The  interesting  observation,  I  believe,  is  that  a  model  using 
stochastic  differential  equations  might  provide  a  vehicle  for  using  negative 
information  as  well  as  positive  information.  Here,  using  negative  informa¬ 
tion  simply  means  using  the  information  that  a  sensor  did  not  gain  contact 
on  a  target. 

Consider  a  charged  particle  moving  through  a  field  of  charge  particles  of  the 
same  and  of  opposite  charge.  The  particle  is  attracted  by  charges  of  oppo¬ 
site  polarity  and  repelled  by  charges  of  like  polarity.  Furthermore,  while 
moving  through  the  field  the  particle  experiences  random  external  forces 
which  tend  to  perturb  its  track. 

Suppose  the  charged  particle  is  the  target.  The  like  charged  particles  are 
non-detecting  sensors  whose  repulsion  forces  are  described  by  the  predicted 
probability  of  non-detection  as  a  function  of  the  range  and  bearing  of  the 
target  from  the  sensor.  The  oppositely  charged  particles  are  detecting 
sensors  whose  attraction  force  is  also  described  by  predicted  probabilities 
of  detection.  The  random  forces  are  other  factors  leading  to  course  and 
speed  changes  of  the  target.  The  analogy  breaks  down  because  of  bearing 
Information  which  requires  that  the  target  be  along  given  bearings  at  given 
times.  The  problem  is  worth  pursuing  and  it  is  our  Intent  to  study  the 
problem  in  the  future. 

The  goal  of  using  stochastic  differential  equations  with  model  inputs  to 
allow  use  of  negative  Information  as  well  as  better  use  of  holding  Informa¬ 
tion  is  not  to  obtain  extremely  accurate  localization  information  but  to  pro¬ 
vide  Information  useful  for  planning  the  allocation  of  mobile  sensors  or 
re-processing  of  data  recorded  in  real  time. 
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4  THE  NATURE  OF  PREDICTIONS  DESIRED  FOR  INFORMATION  PROCESSING 

Performance  predictions  to  be  used  in  information  processing  can  always  be 
based  on  measured  or  recently  measured  ambient  noise  and  clutter. 

When  a  detection  has  occurred,  signal-to-noise  estimates  are  available  as  a 
function  of  time  during  detection  periods.  Autocorrelation  could  be  per¬ 
formed  on  the  received  signal  during  these  periods.  It  might  be  possible 
to  glean  multi-path  date  from  the  resulting  ambiguity  surfaces. 

It  is  rare  that  predictions  for  a  sensor  not  holding  contact  will  be 
required  unless  some  other  sensor  is  holding  contact.  For  the  case  when 
the  sensors  and  targets  are  all  in  a  homogeneous  environment,  it  seems  rea¬ 
sonable  to  try  to  use  what  was  observed  by  the  detecting  sensor  to  generate 
better  predictions  for  the  non-detecting  sensor.  To  my  knowledge,  such 
conditional  predictions  have  not  been  investigated. 

It  is  my  feeling  that  high  operational  payoff  might  result  if  modeling  empha^ 
sis  shifted  from  predictions  based  on  historical  or  synopic  data  to  predic¬ 
tions  based  on  information  derivable  from  directional  noise  measurements  for 
the  sensor  and  related  available  measurements  for  sensors  whose  outputs  can 
be  spectrun  analyzed  or  coherently  processed. 

The  durations  of  holding  and  non-holding  intervals,  signal-to-noise  ratios 
during  holding,  bearing  behavior,  all  may  contain  useful  clues  concerning 
propagation  and  hence  target  range  or  the  probability  that  detections  occur 
or  do  not  occur  on  several  sensors  simultaneously  or  near  simultaneously. 

The  prediction  of  observables  more  general  than  probabilities  of  detection 
should  be  the  direction  that  future  acoustic  modeling  takes  if  it  is  to  be 
a  major  tool  in  information  processing  of  passive  acoustic  sensor  data. 
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THE  PREDICTION  OF  TEMPORAL  STATISTICS  OF  DIRECTIONAL  AMBIENT  SHIPPING  NOISE 
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ABSTRACT 

A  technique  for  predicting  the  probability  density  function  (PDF)  of  the 
ambient-noise  power  received  by  a  narrow-band  system  in  a  modelled  shipping 
and  wind  noise  field  is  described.  The  densities  of  the  shipping  noise  are 
obtained  from  a  mathematical  model  for  spectral- 1 ine  components.  The 
importance  of  the  probability  of  occurrence  of  combinations  of  spectral- 
line  components  in  the  frequency  band  of  interest  is  emphasized.  It  is 
shown  that  the  uncertainty  in  receiving  a  particular  combination  of  lines 
in  the  frequency  band  of  interest  leads  to  the  definition  of  the 
unconditional  PDF  of  the  noise.  Its  prediction  consists  of  the  addition  of 
numerically  convolved  densities  in  conjunction  with  the  probability  of 
occurrance  of  a  line  combination.  The  results  show  that  the  shape  of  the 
resulting  PDF  depends  on  the  processing  bandwidth  and  the  statistics  of  the 
wind/fcackground  noise  level.  Suggestions  are  made  to  extend  the  technique 
to  use  it  for  closer  shipping. 


INTRODUCTION 

The  object  of  this  study  is  the  prediction  of  the  statistics  of  directional 
ambient  shipping  noise  with  a  mathematical -numerical  technique.  This 
reported  work  is  an  extension  of  the  already  existing  deterministic  model 
RANOI  II  <1>  to  which  random  theory  is  added. 

The  technique  for  predicting  the  statistics  is  based  on  the  statistical 
model  published  by  Dyer  <2>  and  gives  a  simple  analytical  expression  for 
the  probability  density  function  (PDF)  of  the  shipping-noise  power  as  a 
function  of  a  combination  of  received  spectral  line  components.  The  total 
short-time  averaged  noise  power  expected  to  be  received  in  the  frequency 
band  of  interest  depends  on  the  processing  bandwidth  relative  to  the 
assumed  distribution  of  the  lines  in  the  noise  spectrum.  For  a  "narrow"- 
band  system  -  narrow  in  the  sense  that  the  bandwidth  is  smaller  than  the 
average  line  spacing  in  a  single-ship  spectrum  **  the  number  of  lines 
present  in  the  band  will  be  on  the  average  less  than  the  number  of  sources. 
The  position  of  the  lines  in  the  spectrum  is  assumed  to  be  random; 
consequently  the  appearance  of  a  combination  of  lines  in  the  frequency 
band  of  interest  will  be  random.  As  all  these  combinations  are  mutually 
exclusive  (t  can  be  shown  that  the  unconditional  PDF  of  the  noise  due  to 
all  possible  combinations  of  spectral  lines  is  a  weighted  sum  of  the 
conditional  POFs  of  the  individual  combi  nations;  with  their  probability  of 
^  1  occurrence  as  weighting  factors.  Results  of  the  prediction  of  the 
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statistics  of  shipping  (plus  wind  noise)  from  a  modelled  shipping  noise 
field  from  RANDI  II  are  presented  and  discussed.  RANDI  II  and  the 
prediction  technique  are  developed  to  be  used  if  the  beam  of  a  sensor  is 
steered  in  the  direction  of  a  shipping  noise  field  that  is  kept  static 
during  the  modelled  observation  time.  This  situation  may  be  achieved  for 
distant  shipping.  Experimental  results  from  closer  shipping  show  a  highly 
variable  shipping-noise  field  due  to  the  moving  of  shipping  in  and  out  of 
the  beam.  This  is  also  observed  in  a  simulated  shipping-noise  field  from 
RANDI  II.  To  solve  this  problem  an  extension  of  the  technique  is 
suggested. 

Section  1  describes  the  modelled  physical  situation  and  gives  typical 
results  using  Dyer's  Statistical  Model.  Section  2  proposes  a  spectral 
model  and  describes  the  uncertainty  about  the  spectral  line  reception, 
which  is  the  majo*'  reason  for  defining  the  unconditional  PDF.  A  block 
diagram  shows  how  Dyer's  Statistical  Model,  the  spectral  model,  and 
RANDI  II  relate  to  its  numerical  calculation.  Results  are  presented  in 
Section  3  and  Section  4  discusses  the  moving  shipping  problem. 


1  MODELLED  PHYSICAL  SITUATION 

The  medium  and  system  modelled  by  the  prediction  system  are  shown  in 
Fig.  1.  The  noise  due  to  shippfng  propagates  through  multipaths  to  the 
sensor.  The  shipping  field  in  the  model  is  assumed  to  be  static  in  the 
sense  that  the  ships  do  not.  move  during  the  measurement  observation  time. 
The  sensor  is  an  array  of  hydrophones.  The  signals  received  in  the  beams 
are  passed  through  a  band  filter.  The  filter  output  is  squared,  short- 
time-averaged  anr  yields  the  beam-power  output.  Figure  1  also  shows  the 
most  Important  models  used  in  the  prediction  system.  Dyer's  statistics  of 
distant  shipping  noise  <Z>  models  distant  shipping  noise  analytically  as 
phase-random  line  components  arriving  from  a  number  of  Independent  sources. 
The  model  is  applicable  for  cases  where  distinct  assumptions  about  the 
random  propagation  mechanism  <3>  and  the  processing  system  are  satisfied. 
In  order  to  allow  the  addition  of  power  intensities  in  most  spectral-line 
combinations,  the  averaging  time-bandwidth  product  must  be  much  larger  than 
unity  (see  App.  A). 

Dyer's  PDFs  and  mathematical  formulas  for  different  spectral  line 
receptions  are  shown  In  Fig.  2,  In  which  Is  the  long-time  average 

intensity  of  each  line  and  l  Is  the  number  of  lines.  The  standard 
deviation  (SDEV)  of  the  level  of  one  line  is  constant,  about  5.6  dB.  If 
more  lines  are  combined,  the  PDF  becomes  narrower  and  the  standard 
deviation  lower.  This  Is  caused  by  the  logarithmic  transformation  of  the 
power  Intensities.  The  prediction  technique  uses  an  algorithm  that 
calculates  the  PDF  for  a  combination  of  lines  with  different  long-time 
average  Intensities  (see  App.  B). 
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2  THE  SPECTRAL  MODEL 


In  the  spectral  model  it  is  assumed  that  the  radiated  noise  of  merchant 
ships  consists  of  spectral-line  components  with  a  constant  line  spacing 
<2,4>.  We  consider  only  those  pass  bandwidths  that  are  smaller  than  the 
line  spacing  (line-spacing/bandwidth  ratio  greater  than  1);  no  more  than 
one  line  per  ship  is  expected  to  be  present  in  the  band,  as  shown  in 
Fig.  3.  To  characterize  the  reception  of  a  distinct  combination  of 
spectral  lines,  assuming  N  merchant  ships  in  the  beam  of  the  array,  we  need 
to  know  the  probability  of  receiving  this  distinct  number  of  lines  in  the 

band.  If  we  assume  that  both  the  position  of  the  lines  in  the  spectrum  and 

the  number  of  lines  in  the  band  of  interest  are  random,  we  define  the 
conditional  probability  of  SL  spectral  lines  given  N  ships:  P  {£/N}.  As 
the  lines  are  received  from  independent  sources,  we  can  use  the  binomial 
distribution  theory  to  calculate  this  probability  (see  App.  C). 

An  indication  for  the  reception  of  different  numbers  of  spectral  lines  in 
the  band  is  obtained  by  observing  the  results  from  SACLANTCEN's  towed-array 
experiments  <4>.  Figure  4  gives  the  beam-noise  level  output  as  a  function 
of  azimuth  angle  for  three  adjacent  frequency  bands  with  widths  of  2.5  Hz. 
The  mean  and  standard  deviation  of  the  noise  from  5G  data  samples  are 

plotted.  The  centre  frequencies  of  the  bands  are  148,  151.5  and  153  Hz. 

Although  in  general  the  mean  values  from  the  three  bands  differ  little 
within  one  beam,  the  standard  deviations  differ  in  some  beams  by  more  than 
3  dB.  Referring  to  Dyer's  PDFs  in  Fig.  2,  this  difference  is  presumably 
caused  by  a  different  combination  of  spectral- line  components  in  the  bands 
of  these  beams.  Because  it  is  unknown  which  combination  of  lines,  out  of 
all  the  possible  combinations,  will  be  present  in  the  band,  we  have  defined 
the  unconditional  PDF  of  the  noise.  The  formula  Is  shown  in  Fig.  5.  This 
PDF  fs  the  sum  of  all  the  conditional  PDFs  of  every  possible  combination  of 
lines,  weighted  by  the  probabilities  of  receiving  a  particular  combination 
of  lines  in  the  band.  Appendix  D  proves  that  the  weighted  summation  of 
PDFs  is  allowed  because  the  line  combinations,  with  which  the  PDFs  are 
connected,  occur  mutually  and  exclusively. 

The  block  diagram,  Fig.  6,  shows  the  relation  of  the  numerical  calculation 
of  the  unconditional  PDF  by  means  of  the  different  models.  Three  flows  of 
Information  are  extracted  from  the  RANDI  II  model  and  directed  to  the 
submodels: 

a)  The  power  Intensities  of  the  line  arrivals  of  N  ships  go  to 
Dyer's  statistical  model, 

b)  The  power  Intensities  of  the  wind/background  noise  and  the 
processing  time-bandwidth  product  go  to  the  wind  statistics 
model , 

c)  The  number  of  ships  and  processing  bandwidth  go  to  the  spectral 
model. 


The  values  of  f0  corresponding  to  the  PDFs  of  all  combinations  of  t 
*  k 


The  resulting  POF,  f. 


lines,  are  summed  and  normalized, 
with  the  PDF  of  the  wind  noise 
weighted  and  summed  over  all  possible  numbers  of  lines. 


f  ;  this  gives  the  PDF, 


is  convolved 
f,  which  is 


A 
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3  RESULTS  FROM  THE  PREDICTION  MODEL 

The  prediction  model  has  been  run  with  shipping  data  acquired  by  maritime 
patrol  aircraft  during  one  of  SACLANTCEN's  towed-array  experiments.  The 
data  set  is  given  in  Appendix  0  <4>.  A  typical  result  of  the  unconditional 
PDF  of  the  noise  from  the  shipping-noise  field  modelled  in  RANDI  II  with 
the  above  shipping  data  is  shown  in  Fig.  7.  This  figure  shows: 

The  wind  noise  level  and  the  noise  levels  of  the  7  spectral  lines 
from  7  ships. 

The  conditional  probabilities  P  {2/7}  (2  =0,..7) 

The  weighted  PDF  curves  for  different  combinations  of  lines  2. 

(2  -•=  0, . . 7). 

-  The  unconditional  PDF  curve. 

In  this  example  the  shape  of  the  resulting  PDF  is  bimodal.  This  is  because 
the  wind  noise  (and  If  the  wind  is  not  present,  the  background  noise)  is 
always  present  in  every  combination  of  spectral  lines.  Every  weighted  PDF 
curve  is  either  that  of  the  wind/background  noise  alone  or  the  result  of  a 
convolution  of  PDFs  of  wind  noise  and  line  noise.  It  is  clear  that  the 
influence  of  the  wind/background  noise  statistics  on  the  unconditional  PDF 
of  ambient  noise  will  be  substantial,  especially  for  higher  frequencies. 

Figures  8,  9,  10  and  11  show  the  results  of  predicting  the  unconditional 
PDF  of  noise  from  the  modelled  shipping  noise  field  for  different  beam 
outputs  and  various  time-bandwidth  products  (TBP)  and  line-spacing/ 
bandwidth  ratios  (LBR). 

The  modelled  beam  power  response  level  vs  azimuth  angle  and  the  uncondi¬ 
tional  PDFs  with  probabilities  of  line  reception  are  plotted  for  some 
beams.  It  needs  to  be  emphasized  that  no  conditions  are  made  about  the 
combination  of  received  lines  in  the  band  of  each  beam.  It  is  therefore 
not  possible  to  compare  these  results  directly  with  the  temporal  statistics 
of  real  data.  The  TBP  of  the  processing  system  Influences  the  statistics 
of  the  wind/background  noise  Intensity.  A  high  TBP  gives  a  narrow  PDF  of 
the  wind  noise.  For  narrow  bandwidths  the  LBR  will  be  high.  Figures  8  to 
11  demonstrate  that  the  Influence  of  the  wind-noise  level  relative  to  the 
total  shipping  noise  level  is  substantial,  especially  in  beams  with  few 
ships.  In  most  of  the  beams  the  shape  of  the  PDFs  is  bimodal.  The  range  of 
standard  deviation  Is  large. 

Vie  can  draw  the  following  general  conclusions  from  the  above  results: 

1)  The  unconditional  PDF  of  the  noise  level  is  non-gauss Ian  for  most 
of  the  beams. 


2)  In  the  case  of  narrow  bands  the  influence  of  the  wind-background 
noise  will  be  substantial  in  beams  with  few  ships,  because  of 
the  high  probabilities  of  receiving  no  lines  in  the  band. 

3)  low-level  wind-background  noise  leads  to  higher  standard 
deviations  than  high-level  wind-background  noise. 
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4)  Variation  of  the  bandwidth  influences  the  amplitude  of  the 
resulting  PDF  more  than  its  standard  deviation. 

5)  Wind-background  noise  with  a  narrow  PDF  dominates  the 
unconditional  PDF,  particularly  for  narrow  bandwidths. 


4  THE  MOVING  SHIPPING  PROBLEM 

The  above  proposed  technique  is  developed  to  predict  the  statistics  of 
ambient  shipping  noise  from  distant  shipping.  If  the  shipping  is  closer, 
the  influence  of  moving  shipping  will  have  substantial  effect  on  the 
statistics  of  the  ambient  shipping  noise.  The  varying  shipping  noise  field 
can  be  observed  in  the  experimental  data  and  can  be  modelled  with  RANDI  II. 
Figures  12a  and  b  plot  the  beam  noise  level  outputs  for  three  series  of 
successive  and  overlapping  noise  samples  from  towed-array  measurements. 
Three  values  of  mean  and  standard  deviations  of  the  noise  are  plotted 
against  azimuth  angle  for  150  and  502  Hz.  Beam  outputs  with  large 
variation  in  the  three  mean  values  have,  in  general,  high  standard 
deviations,  which  might  be  explained  by  ships  moving  in  and  out  of  these 
beams  during  the  measurement  period.  Another  indication  for  moving 
shipping  is  simulated  with  RANDI  II.  Figure  13a  shows  the  shipping  noise 
field  calculated  from  the  shipping  surveillance  data  (see  App.  D  <4>) 
collected  by  maritime  patrol  aircraft  during  the  measurements  shown  in 
Fig.  12.  The  spikes  are  shipping-noise  levels  in  one-degree  sectors.  The 
curve  is  the  noise  field  convolved  with  an  unambiguous  beam  of  5°. 
Figure  13a  shows  the  field  at  the  start  of  an  observation  time  period  and 
Fig.  13b  shows  the  field  at  the  end  of  the  period.  A  highly  variable  noise 
field  is  observed  due  to  moving  shipping.  The  impact  on  the  corresponding 
modelled  beam  noise  responses  is  shown  in  Figs.  14a  and  b. 

In  order  to  handle  moving  shipping  with  a  prediction  technique  a  dynamic 
shipping-noise  field  should  be  used.  The  statistics  of  the  noise  could  be 
calculated  out  of  noise  samples  from  a  varying  shipping-noise  field.  This 
solution  would  be  beyond  the  scope  of  this  study,  In  which  an  analytical 
calculation  fs  required.  An  extension  of  the  prediction  technique  is 
suggested  with  the  use  of  unconditional  statistics  (Fig.  15).  We  propose 
the  unconditional  POF  of  the  noise  due  to  all  possible  combinations  of  M 
ships  in  the  beam  during  the  observation  time  period.  This  PDF  is  the 
weighted  sum  of  the  unconditional  PDFs  due  to  all  possible  combinations  of 
spectral  lines  in  the  band,  given  a  particular  combination  of  N  ships  In 
the  beam  of  the  sensor.  The  weighting  factors  are  the  probabilities  of 
occurrence  of  a  particular  combination  of  ships  in  the  beam.  With  this 
extension,  which  is  under  development,  the  predicted  PDFs  are  expected  to 
be  much  wider,  with  much  higher  standard  deviations.  This  will  be  caused 
by  the  extra  uncertainty  about  the  shipping  covered  by  the  beam. 
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APPENDIX  A 

JUSTIFICATION  OF  THE  USE  OF  POWER  ADDITION  IN 
THE  PROPOSED  PREDICTION  TECHNIQUE 


If  the  output  of  the  band-pass  filter  in  Fig.  16  is  represented  by  the 
signal 


V(t)  =  1  Akcos(u>Rt  +  0k) 


as  a  sum  of  line  signals  with  amplitude  Ak,  frequency  fk  -  wk/2n,  and  phase 
0k,  it  can  be  shown  that  the  squared  and  short-time  averaged  (interval  I 
seconds)  signal  v(t)  becomes: 


P(t)  =  \ 


&  i  A.  A  Ul/2 

1  ±  —  [  du  cos  |  (u>. 

k=l  m=l  1  t- 1/2  K 

i! 

I  A?  +  off-diagonal  terms 
k=l  K 


uj  )u  +  0.  -  0  I 
nr  k  m' 


These  off-diagonal  terms  are  S  0  if 


min 

k#m 


w. 


and  P(t)  is  a  sum  of  line-power  intensities. 

If  we  define  the  minimum  relative  line  spacing  a  as 


A  _  min  f  . 
A  "  k*ra  Tk  ‘  rm 


we  get  the  condition  I  >  1/A 


If  we  use  a  time-bandwidth  product  much  larger  than  1,  then 


IW 


(Eq.  A.l) 


For  a  bandwidth  that  is  not  too  small  and  a  few  spectral  lines  received  in 
the  band,  it  is  allowable  to  use  the  addition  of  spectral  line  power  levels 
if  the  condition  A.l  is  fulfilled.  If  we  use  a  line- spacing/bandwidth 
ratio  of  5  and  a  maximum  possible  reception  of  7  spectral  lines  in  the  pass 
band,  a  relative  line  spacing  less  than  the  minimum  A  will  have  a 
probability  of  occurrence  of  less  than  10*. 


J&* 
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APPENDIX  B 


The  derivation  of  the  PDF  for  unequal  long-time  average  intensities  is 
shown  in  <2>.  An  analysis  of  the  characteristic  functions  of  PDFs  is  used. 
An  evaluation  by  the  calculus  of  residues  and  transformation  to  the 
logarithmic  domain  yields 


£  _i_i 

f(x)  =  l  M£  P£  e*p[x  -  exp(x/pm)J 
m-i  11 


where 


-  A 11 ' 


with  u  the  long-time  average  intensity  of  the  mth  line  and  £  the 
number 


of  lines  that  are  combined.  The  property  IP  =  1  is  used  in  the  pre¬ 
nyl  * 

diction  technique  as  a  quality  measure  for  the  numerical  calculation. 


APPENDIX  C 


We  assume  N  lines  in  a  spectral  line  space  divided  into  m  bands,  where  m  is 
the  line-spacing/bandwidth  ratio.  If  one  line  is  randomly  present  in  the 
space,  the  probability  that  this  line  is  present  in  a  marked  band  is  1/m 
and  the  probability  that  it  is  not  present  in  this  band  is  (m-l)/m. 


Since  all  the  N  lines  are  independent  and  not  labelled,  we  use  the  binomial 
distribution  to  calculate  the  probability  that  £  lines  out  of  N  are  present 
in  the  band  of  interest: 


(!)  (=£) N-4  (0 


! 

4 


SACLANTCEN  CP-32 


20-7 


GROEN:  Temporal  statistics  of  shipping  noise 


APPENDIX  D 

PROOF  WEIGHTED  SUMMATION  OF  DENSITIES 


Suppose  v.  is  the  event  that  the  value  of  the  random  variable  (r.v)  X  is 
greater  titan  or  equal  to  x;  and  less  than  or  equal  to  x..+p  the  probability 

of  event  v.  is: 

P^}  =  P{X.  <  X  <  xi+1}  =  J1+1  f(x)dx 

xi 

with  f(x)  being  the  PDF  of  r.v.  X. 

Further,  suppose  that  is  the  event  that  the  combination  of  spectral 
lines  is  received  in  the  band  of  interest,  with  probability  P { c ^ . 

If  v^  is  the  consequence  of  c^,  then  the  probability  of  event  VjCj  is 

PWjC.}  =  P{v./Cj}P{Cj}  (Eq.  D.l) 

The  events  in  the  set  v.Cp  v.c2,...,  are  mutually  exclusive  because 
Cp  c2,...,  ck  are  mutually  exclusive.  The  subscript  kT  is  defined  as  the 
total  number  of  combinations  of  lines. 

From  above  and  Eq.  D.l  we  obtain 


P{v 


or 


tCl  *  VIC2  *  '•  *  Pl^/CjJPfCji 


Ic  X 

P{V1(c1  ♦  C2  +  ...  ♦  Ck  )}  a  IT  ;1+1  fc  (X)  dx  P{Cj} 

1  jsl  X,  j 


or 


/,1'1  fc  ♦  ,  ♦  +  ,  (X)  dx  =  iT  /,tl  f  (x)  dx  P|c,i  (Eq.  D.2) 

x,  cltc2*"  +ckT  1=1  x,  CJ  J 


with  f.  (xj  being  the  probability  density  of  X  doe  to  the  f™  combination 

CJ 
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of  lines,  and  f  (x)  being  the  probability  density  of  X  due 

c  1  c  2  ••  c  ky 

to  "sum"  of  all  the  possible  combinations  of  lines. 

In  numerical  calculations  for  Ax  sufficiently  small,  we  obtain 
P{vi)  =  P{ x,j  <  X  <  x^  +  Ax}  S  fCx^Ax 


From  the  above  and  Eq.  D.2  we  conclude  that 


or 


+  C2  + 


(x.)Ax  =  I  f  (x.)  Ax  P{c.} 
kT  4-1  cj  '  J 


+  c. 


<x1> 


I  f,  (xf)PUJ 
J«1  j  J 


(Eq.  0.3) 


Hence  a  sample  of  the  total  PDF  of  the  noise-power  level  X  is  the  weighted 
sum  of  the  PDFs  due  to  all  possible  combinations  of  spectral  lines  with  the 
probability  of  occurrence  of  the  different  combinations  of  lines  as 
weighting  factors. 
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AMBIENT  NOISE  IN  THE  NORTHEASTERN  PACIFIC  OCEAN 

by 
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Science  Applications,  Inc. 

Canoga  Park,  Ca.,  USA 


ABSTRACT 


One-hour  averages  of  omnidirectional  ambient  noise  measurements  at  60  Hz  and 
165  Hz  are  analyzed  for  two  nearby,  deep  ocean  sites  in  the  Northeastern 
Pacific  during  February  and  March  1981.  Site  A  is  a  high  noise  site  and  is 
located  near  major  east-west  shipping  lanes  and  near  major  Pacific  storm 
paths.  Site  B  is  a  lower  noise  site  and  is  located  approximately  450  n.m. 
from  Site  A  away  from  major  shipping  lanes  and  Pacific  storm  paths.  The 
site  and  frequency  dependence  of  ambient  noise  is  found  to  be  highly  variable 
with  shippiing  noise  being  totally  dominant  at  60  Hz  at  Site  A  and  storm 
noise  being  totally  dominant  at  165  Hz  at  Site  B.  Both  shipping  and  storm 
noise  can  dominate  the  165  Hz  Site  A  or  60  Hz  Site  B  noise  levels  depending 
on  weather  conditions. 

Storm  noise  has  a  possible  Indirect  effect  on  shipping  noise,  since  very  low 
noise  periods  occur  In  between  storms,  especially  when  the  storm  passes  near¬ 
by  the  site.  A  limited  number  of  WMO  ships  sampled  Indicate  that  ships  slow 
down  or  stay  In  port  during  storms.  Average  omnidirectional  noise  levels  at 
60  Hz  at  both  sites  and  at  165  Hz  at  Site  A  were  3  to  4  dB  lower  during  a 
stormy  week  than  they  were  during  a  relatively  calm  week. 

Array  noise  gain  measurements  at  Site  B  Indicate  that  the  coherence  of  noise 
curing  stormy  periods  Is  much  less  than  It  Is  during  calm  periods.  Gener¬ 
ally,  this  Implies  that  the  Increase  in  beam  level  at  165  Hz,  caused  by 
storm  noise,  will  be  significantly  less  than  corresponding  increase  In  omni¬ 
directional  noise  levels.  At  60  Hz  at  Site  B,  this  Implies  that  beam  levels 
will  decrease,  In  general,  during  storms,  especially  If  the  storm  system  Is 
nearby. 


1.0  INTRODUCTION 

There  are  many  sources  of  ambient  noise1  In  the  ocean  and  it  is  shown  In 
this  paper  that  different  noise  sources  (shipping  and  storms  In  this  case) 
can  dominate  the  noise  characteristics  of  two  nearby  sites.  For  the  anal¬ 
ysis  presented  here-  it  is  Important  to  differentiate  between  short  term 
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(minutes)  noise  fluctuations  and  changes  in  the  long  term  (hours)  noise 
intensity  level.  The  two  major  factors  affecting  long  term  noise  intensity 
levels  are  the  distribution  of  noise  sources  relative  to  the  measurement 
site  and  transmission  loss  from  each  noise  source  to  the  measurement  site. 
The  data  analyzed  in  this  paper  illustrate  the  measurement  site  dependence 
and  frequency  dependence  of  long  term  averaged  ambient  noise.  Bottom  mount¬ 
ed  hydrophone  data  were  collected  from  two  sites  and  two  frequencies  in  deep 
water  in  the  Northeastern  Pacific  Ocean  during  February  and  March  1981.  The 
data  used  in  this  analysis  are  averaged  over  1  hour  at  60  Hz  and  165  Hz. 
Short  term  noise  fluctuations  will  not  be  addressed  in  this  paper. 

Two  dominant  sources  of  noise  in  the  Northeastern  Pacific  Ocean  during  the 
winter  are  shipping  and  storms.  Site  A  is  located  near  major  Pacific  ship¬ 
ping  lanes  and  is  in  the  path  of  many  large  Pacific  storms.  Site  B  is 
located  approximately  450  miles  from  Site  A  and  is  located  well  south  of 
major  shipping  lanes  and  Pacific  storms.  Consequently,  Site  A  is  much 
noisier  than  Site  B  and  the  60  Hz  noise  levels  at  Site  A  were  totally  domin¬ 
ated  by  shipping  noise.  On  the  other  extreme,  the  165  Hz  levels  at  Site  B 
were  totally  dominated  by  storm  noise.  The  objective  of  this  paper  is  to 
show  that  low  frequency  ambient  noise  can  be  extremely  site  and  frequency 
dependent  and  to  show  the  large  changes  in  array  noise  gain  that  can  occur 
in  a  period  of  1  day. 

To  differentiate  between  storm  and  nearby  shipping  events  in  the  ambient 
noise  time  series,  the  following  criteria  were  used:  Storm  events  were 
defined  as  increases  in  noise  level  lasting  from  2  to  7  days  at  10  Hz  and 
165  Hz  with  little  or  no  change  in  noise  level  at  60  Hz.  Nearby  shipping 
events,  on  the  other  hand,  are  defined  as  increases  in  ambient  noise  level 
at  all  frequencies  from  10  Hz  to  165  Hz  lasting  1  or  2  hours.  Distant  ship¬ 
ping  noise  background  has  a  maximum  near  60  Hz,  and  since  the  shipping 
spectra  is  relatively  flat,  storm  events  have  a  greater  effect  on  10  Hz  and 
165  Hz  noise  levels  than  on  60  Hz  levels.  The  10  Hz  noise  levels  were  used 
only  to  define  storm  events  and  no  analyses  of  these  data  are  presented  in 
this  paper.  Every  storm  event  defined  in  the  manner  described  above  was 
accompanied  by  a  major  Pacific  storm  and  it  is  unlikely,  In  the  author's 
opinion,  that  shipping  can  produce  the  storm  noise  characteristics  defined 
above  (and  vice  versa). 

A  subsequent  paper  will  cover  the  subjects  of  distant  versus  local  storm 
noise  and  modeling  ambient  noise  based  on  Wilson's  wind  source  level 
curves2""  and  historical  shipping  densities.  This  paper  covers  ambient 
noise  characteristics  of  Sites  A  and  B»  emphasizing  the  frequency  and  site 
dependence  of  ambient  noise. 

2.0  AMBIENT  NOISE  CHARACTERISTICS 

The  data  base  for  this  analysis  is  60  Hz  and  165  Hz  omnidirectional  data 
from  two  deep  water,  Northeastern  Pacific,  bottom  mounted  hydrophones.  A 
1-hour  (dB)  averaging  time  was  used  In  order  to  analyze  long  term  changes 
In  ambient  noise  time  series  at  each  site.  Wind  speed  information  was  taken 
from  the  Fleet  Numerical  Oceanographic  Center  (FN0C)  surface  wind  analyses 
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performed  every  6  hours  (Zulu  time).  The  accuracy  of  the  FNOC  surface  wind 
analyses  for  the  mid-latitudes  was  evaluated  in  the  NORPAC  Pole  Experiment5. 


2. 1  Site  A  -  The  High  Noise  Site 

Figures  1  and  2  show  the  60  Hz  and  165  Hz  omnidirectional  ambient  noise  data 
at  Site  A  for  the  month  of  February  1981.  Data  points  were  plotted  every  12 
hours  (Zulu  time)  except  on  12/13  February  1981,  when  no  data  were  available. 
Four  major  Pacific  storms  passed  to  the  north  of  Site  A  during  6  to  8  Feb¬ 
ruary,  16  February,  18  to  20  February  and  21  to  27  February  1981,  while  a 
small  local  squall  passed  overhead  Site  A  on  10  February  1981. 

Except  for  the  high  levels  on  10,  18-19  and  24  February  1981,  Figure  1  shows 
little  effect  of  the  passage  of  the  Pacific  storms  at  60  Hz.  It  can  be  con¬ 
cluded  that  the  ambient  noise  at  Site  A  is  dominated  by  shipping  the  vast 

majority  of  the  time.  The  60  Hz  levels  of  85  dB/pPa//Hz1*,  shown  in  Figure  1, 

are  high  enough  (relative  to  the  165  Hz  levels  at  Sites  A  and  B)  to  limit 
storm  effects  to  2  to  4  dB  increases  lasting  periods  of  1  day  or  less. 

The  5  to  7  dB  level  increases  at  165  Hz  shown  in  Figure  2  during  storm  per¬ 
iods  are  much  more  pronounced  with  the  high  levels  lasting  up  to  5  days  in 

duration.  The  fact  that  these  high  levels  correlate  very  well  with  the 
presence  of  large  Pacific  storms  is  a  direct  result  of  noise  produced  by 
storms.  The  indirect  effect  of  storms  is  also  potentially  significant. 

Figure  2  suggests  that  very  low  levels  (e.g.,  9,  11-12  and  21-22  February 
1981)  can  occur  at  165  Hz  between  storms.  Figure  1  shows  the  same  effect  at 
60  Hz,  but  to  a  much  lesser  extent.  The  question  that  arises  is  "Does  the 
presence  of  storms  reduce  shipping  noise  by  causing  ships  to  slow  down  or 
change  course?" 

The  author  has  not  done  an  extensive  study  of  ship  course  or  speed  changes 
during  storm  periods,  but  several  WMO  ships  were  tracked  that  were  in  the 
Northeastern  Pacific  during  storm  periods  on  6  to  8  and  10  February  1981. 

The  source  of  WMO  ship  position,  course  and  speed  data  are  the  WMO  ship 
reports  received  by  FNOC,  Montery.  Figure  3  shows  the  tracks  of  five  WMO 
ships  (with  reported  speed  time  series  shown  in  the  Insert)  from  5  to  11 
February  1981.  These  ships  were  selected  because  they  reported  frequently 
and  were  located  in  the  east-west  shipping  lane  that  crossed  the  storm 
center.  Although  the  WMO  ship  list  does  not  Include  all  ships  at  sea,  sev¬ 
eral  characteristics  of  WMO  ship  operations  observed  during  or  after  storms 
may  be  typical  of  most  ships.  In  examining  the  WMO  ship  locations,  courses 
and  speeds,  the  following  observations  are  made: 

o  All  WMO  ships  during  this  time  period  slowed  down  within  (e.g., 
JBSP)  or  near  (e.g.  H8DE)  the  major  storm  system. 

o  Not  all  WMO  ships  slowed  down  at  the  same  time,  since  location 
relative  to  storm  systems  was  the  major  factor  in  speed  changes. 

o  WHO  ships  that  did  not  cross  the  storm  system  (e.g.,  D5HD) 
tended  to  maintain  constant  course  and  speed. 
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o  Significant  long  term  changes  in  WMO  ship  course  were  not 

observed.  The  preferred  method  of  getting  around  a  storm  seems 
to  be  slowing  down  and  "waiting"  for  calm  weather. 

o  The  vast  majority  of  WMO  ships  during  this  time  period  were  on 
easterly  courses  with  very  few  WMO  ships  leaving  west  coast 
ports  on  westerly  trans-Pacific  routes. 

The  above  observations  were  made  on  a  very  limited  time  period  for  WMO  ships 
only  and,  therefore,  a  much  more  extensive  analysis  is  needed  to  say  these 
are  typical  of  ship  course  and  speed  changes.  However,  all  the  observations 
imply  less  ship  noise  during  and  between  storms.  This  could  be  a  possible 
explanation  for  the  low  noise  levels  observed  in  Figure  2.  After  a  storm 
has  passed,  it  could  take  approximately  a  day  for  shipping  densities  and 
speeds  to  increase  to  normal  levels. 

The  week  from  5  to  11  February  1981  was  a  typical  storm  period  for  the 
Northeastern  Pacific,  with  one  normal  and  one  small  storm  during  the  week. 
The  period  from  20  to  27  March  1981  was  an  extended,  unusually  calm  weather 
period  for  the  Northeastern  Pacific  during  winter,  with  one  storm  beginning 
on  the  last  day  of  the  week.  The  period  from  20  to  26  March  1931  was  ship¬ 
ping-dominated  and  a  comparison  of  the  1-hour  averaged  noise  statistics 
between  the  2  weeks  (5  to  11  February  1981/20  to  27  March  1981)  is  shown  in 
Figure  4.  On  the  average,  noise  levels  during  the  calm  March  week  are  3  dB 
and  4  dB  higher  at  65  Hz  and  165  Hz,  respectively,  than  corresponding  noise 
levels  during  the  stormy  February  week.  Again,  the  possibility  of  storms 
lowering  the  average  noise  level  between  60  Hz  and  165  Hz  is  suggested  for 
high  noise  sites  like  Site  A.  Certainly,  Figure  4  shows  there  are  many 
more  very  low  noise  periods  during  the  stormy  week.  It  is  interesting  to 
note  the  blmodal  distribution  of  noise  levels  at  165  Hz  during  the  stormy 
week  shown  In  Figure  4.  Physically,  this  occurs  because  the  noise  levels 
during  stormy  and  calm  periods  are  typically  6  dB  apart  and  there  is  rela¬ 
tively  little  time  in  transition  between  calm  and  stormy  periods. 

Although  more  Intensive  analysis  of  the  effect  of  storms  on  shipping  noise 
Is  needed,  the  data  presented  here  suggests  the  possibility  that,  for  high 
noise  areas  like  Site  A,  storms  may  produce  lower  overall  noise  levels 
between  60  Hz  and  165  Hz  by  reducing  ship  noise. 

2.2  Site  B  -  The  Low  Noise  Site 

As  discussed  previously.  Site  B  is  located  approximately  450  n.m.  from  Site 
A  and  is  well  south  of  the  major  east-west  shipping  lanes  and  major  Pacifc 
storms.  Consequently,  noise  levels  at  60  Hz  and  165  Hz  at  Site  B  are,  on 
the  average,  approximately  11  dB  and  14  dB  quieter  during  calm  weather 
periods  than  corresponding  levels  at  Site  A.  Figures  5  and  6  show  the  omni¬ 
directional  levels  for  Site  B  during  February  1981  at  60  Hz  and  165  Hz, 
respectively.  The  data  are  1-hour  averages  plotted  every  12  hours  Zulu 
time. 
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In  comparing  the  data  in  Figures  5  and  6  with  the  data  in  Figures  1  and  2, 
some  differences  between  Sites  A  and  B  (other  than  average  level)  are  evi¬ 
dent.  The  165  Hz  levels  at  Site  B  are  totally  dominated  by  storm  noise  with 
level  increases  from  9  dB  to  13  dB  during  every  February  1981  storm  period 
described  previously.  During  the  calm  period  between  storms,  the  165  Hz 
levels  at  Site  B  are  very  low,  approaching  59  to  60  dB/pPa/ZHz^  as  a  "noise 
floor."  At  these  very  low  noise  levels,  shipping  noise  is  not  a  factor 
since  a  10  or  15  knot  local  wind2*3  is  sufficient  to  dominate  these  noise 
floor  levels. (  A  large  distant  storm  can  raise  the  noise  level  above 
70  dB/pPa//Hz'2  and  Figure  6  hows  that  the  omnidirectional  165  Hz  levels  at 
Site  B  act  as  a  "weather  barometer"  for  Pacific  storms.  Figure  7  shows  the 
time  series  at  165  Hz  at  Site  B  for  the  week  of  13  to  21  February  1981, 
which  illustrates  the  major  noise  sources  (local  shipping,  local  storm, 
distant  storm  and  calm  or  distant  shipping  background). 

The  60  Hz  levels  at  Site  B  shown  in  Figure  5  vary  much  less  (2  to  4  dB)  with 
storm  events  than  the  165  Hz  levels.  This  indicates  that  the  60  Hz  levels 
are  dominated  more  by  shipping  noise  than  storm  noise  most  of  the  time. 
However,  the  60  Hz  levels  at  Site  B  have  more  variability  than  the  60  Hz 
levels  at  Site  A  (Figure  1)  which  were  totally  dominated  by  shipping. 

It  appears  that  the  60  Hz  levels  at  Site  B  are  dominated  by  relatively 
nearby  ships,  since  the  lowest  levels  in  February  1981  occurred  on  11  Feb- 
rary  1981,  just  after  the  small  local  storm  on  10  February  1981.  The  other 
more  distant  storms  in  February  1981  did  not  result  in  low  levels  between 
storms. 

Figure  8  shows  a  comparison  of  the  1-hour  averaged  statistics  at  Site  B  for 
the  same  storey  week  (5  to  11  February  1981)  and  the  same  calm  week  (20  to 
27  March  1981)  used  in  Figure  4  for  Site  A.  The  average  60  Hz  levels  at 
Site  B  increased  by  approximately  3  dB  during  the  calm  week,  possibly  be¬ 
cause  shipping  densities  and  speeds  increased  to  normal  levels  when  no 
storms  were  present.  However,  the  average  165  Hz  levels  at  Site  B  decreased 
by  5  to  6  dB  during  the  calm  week  due  to  decreased  wind  noise.  The  high 
noise  levels  in  the  bimodal  distribution  shown  in  Figure  8  for  165  Hz  at 
Site  B  during  the  calm  week  are  due  to  a  storm  the  last  day  of  that  week. 

3.0  ARRAY  GAIN 

Analysis  of  the  omnidirectional  data  in  the  previous  section  illustrated  site 
and  frequency  variability  of  ambient  noise.  However,  the  coherence  of  the 
noise  field  is  an  extremely  important  factor  in  assessing  array  performance 
and  omnidirectional  data  is  not  impacted  by  noise  coherence.  Array  noise 
gain  (ANG)  measurements  were  collected  for  a  near  broadside  beam  for  both 
Sites  A  and  B  during  part  of  the  February/March  1981  data  collection  period. 
To  eliminate  measurement  system  effects,  theoretical  array  noise  gain, 
ANGtheo»  was  subtracted  from  measured  array  noise  gain  ANG^eas  as  follows: 

HiE  ‘  mmn  •ANGth<io  (1> 
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The  objective  of  this  section  is  to  illustrate  the  highly  variable  noise 
source  dependence  of  ANG  by  eliminating  system  gain  against  an  ideal  noise 
field.  The  results  discussed  for  ANG  depend  on  the  coherence  of  the  dom¬ 
inant  noise  source  only  and  apply  to  any  array  system. 

If  the  noise  field  were  isotropic  and  totally  incoherent,  ANG"  would  be  zero 
for  all  beams  and  array  gain  would  be  constant  and  equal  the  theoretical 
array  gain  value.  For  Site  A,  and  especially  Site  B,  ARE  was  highly  vari¬ 
able  and  dependent  on  the  noise  source.  Physically,  ANG  is  expected  to  be 
nearly  euqal  to  ANGt(le0  when  the  dominant  noise  source  is  a  large,  distant 

Pacific  storm.  Since  a  typical  Pacific  storm  has  a  fetch  (greater  than  20 
knot  winds)  on  the  order  of  2000  n.m.  (2/3  the  size  of  the  Continental 
United  States),  it  should  appear  more  incoherent  and  isotropic  than  a  noise 
field  dominated  by  nearby  surface  ships.  Figure  9  shows  ANG  at  Site  B 
plotted  for  several  weeks  in  February/March  1981  with  the  165  Hz  omnidirec¬ 
tional  levels  at  Site  B  used  as  an  index  or  "storm  barometer."  ARE'  at 
165  Hz  is  anti  correlated  with  the  165  Hz  levels  probably  because  distant 
storm  noise  is  uncorrelated.  The  noise  field  during  calm  periods  may  be 
dominated  by  a  relatively  few  nearby  ships  and,  consequently  would  be  more 
highly  correlated.  Since  Site  B  is  a  low  noise  site,  it  is  probably  not 
Impacted  greatly  by  distant  shipping  noise.  The  60  Hz  ANG  shipping  values 
shown  in  Figure  9  support  this  supposition.  The  storm  on  16  to  18  March 
1981  passed  very  close  to  Site  B  and  the  60  Hz  ARE  values  were  also  anti- 
correlated  with  the  165  Hz  levels.  Other  more  distant  storms  did  not 
result  in  as  great  a  decrease  in  ARE.  For  example,  the  60  Hz  ANE  values 
did  not  decrease  as  much  during  the  distant  storm  on  1  to  3  April  1981  as 
it  did  during  the  more  local  storm  of  16  to  18  March  1981. 

The  ARE  values  at  60  Hz  and  165  Hz  may  have  different  performance  Implica¬ 
tions  for  an  array  of  elements  when  the  changes  In  ARE  plus  the  omnidirec¬ 
tional  level  are  considered.  At  165  Hz,  the  results  Imply  that  beam  level 
Increases  during  storm  periods  will  not  be  nearly  as  large  as  the  increases 
in  omnidirectional  noise  levels.  In  other  words,  the  coherence  of  storm 
noise  relative  to  the  coherence  of  the  noise  field  during  calm  weather 
periods  results  in  beam  level  increases  that  are  much  less  than  the  omni¬ 
directional  level  increases.  At  60  Hz,  however,  the  situation  is  different 
since  the  60  Hz  omnidirectional  noise  levels  at  Site  B  remain  constant  or 
increase  very  little  during  storm  periods.  Since  ARE  generally  decreased 
more  than  the  omnidirectional  level  increase,  the  storm  indirectly  resulted 
in  lower  beam  levels  at  60  Hz  at  Site  B.  Changes  in  ship  speeds  and  densi¬ 
ties  are  possible  explanations  for  lower  beam  levels. 
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FIG.  J  SELECTED  WMO  SHIP  TRACKS 


SITE  A  60HJ-STORMY 


SITE  A  165HZ-STORMY 


SITE  A  SOHz-CALM  SITE  A  185HZ-CALM 


FIS.  4  ONE-HOUR  AVERAGE  NOISE  STATISTICS  FOR  SITE  A 
FOR  ACALN  AND  STORMY  WEEK 
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Time  February  1981 

FIG.  7  SITE  B  -  165  Hz  TIME  SERIES 


SITE  B  60HZ--STGRMY  SITE  B  165HZ-STORMY 


SITE  B  60HZ-CALM  SITE  B  165HZ-CALM 


SPECTRUM  LEVEL  IN 

FIG.  8  ONE-HOUR  AVERAGE  NOISE  STATISTICS  FOR  SITE  S 
FOR  A  CALM  AND  STORMY  MEEK 
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SITE  DEPENDENCE  OF  WIND-DOMINATED  AMBIENT  NOISE  IN  SHALLOW  WATER 

by 

S.N.  Wolf  and  F.  Ingenito 
Naval  Research  Laboratory 
Washington,  D.C.  20375 


ABSTRACT 


A  survey  of  ambient  noise  data  in  shallow  water  taken  under  wind-dominated 
conditions  shows  substantial  differences  in  spectrum  level  (often  greater 
than  10  dB)  under  the  same  windspeed  and  sea-state  conditions.  This  area- 
dependent  effect,  which  is  evident  even  in  long-term  averaged  data,  is 
caused  by  differences  in  ocean  bottom  properties,  water  depth,  and  sound- 
speed  profile.  Comparisons  between  those  data  sets  which  are  accompanied 
by  suitable  environmental  information  and  a  wind-noise  model  developed 
jointly  at  NRL  and  SACLANTCEN  are  shown.  The  areas  in  which  the  data  were 
taken  cover  a  broad  range  of  bottom  types,  water  depths,  and  sound-speed 
profiles.  Calculated  differences  In  noise  level  agree  well  with  the 
measured  differences  and  it  is  concluded  that  the  model  is  capable  of 
predicting  the  effect  of  the  shallow  water  environment  on  the  spectrum 
level  of  wind-dominated  ambient  noise.  We  also  conclude  that  commonly 
used  universal  ambient  noise  curves  which  do  not  take  into  account 
environmental  differences  are  inferior  to  the  model  predictions  in 
shallow  water. 


Introduction 

In  this  paper  we  present  the  results  of  some  recent  studies  of  the 
dependence  of  wind-driven  ambient  noise  level  In  shallow  water  on  environ¬ 
mental  parameters  other  than  windspeed  and  sea  state.  Among  these  parameters 
are  water  depth,  sound  speed  profile  and  ocean  bottom  composition.  The  paper 
will  begin  with  the  presentation  of  some  measurements  which  suggest  that 
an  observable  difference  in  wind-driven  ambient  noise  level  exists  for 
different  sites  under  the  same  windspeed  conditions.  The  model  used  In 
the  calculation  of  the  effect  of  the  environment  on  noise  level  Is 
described.  We  then  present  the  results  of  the  application  of  this  model 
to  several  data  sets  In  order  to  account  for  observed  level  differences. 
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Shallow  water  wind  driven  noise  measurements  made  at  a  number  of 
locations  [1-9]  are  shown  in  Figs.  1  and  2.  Since  the  measurements,  in 
general,  were  not  made  under  exactly  the  same  windspeed  conditions,  the  data 
are  assembled  by  reported  sea  state,  or  by  sea  state  associated  with  the 
reported  windspeed  by  using  Vine  and  Volkmann's  nomogram  [10].  The  arbitrary 
division  of  a  continuum  of  possible  windspeeds  into  a  small  number  of  sea 
states,  the  subjective  nature  of  estimating  sea  states,  and  the  introduction 
of  measurement  error  are  expected  to  introduce  some  differences  in  the 
measured  spectra.  The  data  show  spectrum  level  differences,  which  in  some 
cases  exceed  10  dB.  These  differences  are  certainly  too  large  to  attribute 
to  usual  uncertainty  of  measurements  of  acoustic  noise  level.  If  we  use 
Piggot's  [2]  measurements  as  a  guide,  we  find  that  a  difference  of  10  dB  in 
ambient  noise  level  corresponds  approximately  to  a  tripling  of  the  wind- 
speed,  or  a  change  from  sea  state  2  to  sea  state  5.  Although  the  comparison 
of  shallow  water  ambient  noise  levels  by  grouping  the  data  by  associated 
sea  state  is  not  "exact",  differences  in  level  appear  which  cannot  be 
explained  by  usual  errors  in  measurement  or  sea  state  estimation. 

Another  comparison  of  wind-generated  ambient  noise  [1]  is  shown  in 
Fig.  3.  In  this  illustration  we  show  noise  spectra  measured  at  two  sites,  one 
having  a  sand  bottom,  the  other  having  a  silt  bottom  along  with  Piggot's  [2] 
measurements.  The  present  measurements  at  frequencies  above  100  Hz  show 
spectral  shapes  similar  to  those  obtained  by  Piggot  and  were  free  from  lines 
or  other  indications  of  shipping  noise,  and  thus  appear  to  be  wind-driven 
noise.  We  see  that  differences  in  level  (~5  dB)  in  the  two  sets  of  measure¬ 
ments  over  the  sand  bottom,  are  similar  to  the  differences  in  Piggot's 
measurements  at  the  corresponding  wind  speeds,  but  the  absolute  levels 
made  over  the  sand  bottom  are  somewhat  higher.  The  measurements  made  over 
the  silt  bottom  are  somewhat  lower  than  those  reported  by  Piggot.  In 
addition,  we  can  see  that  the  measurements  made  at  windspeed  7  m/sec  over 
the  silt  bottom  are  2  dB  lower  than  those  made  over  the  sand  bottom  at 
windspeed  3-5  m/sec.  These  data  provide  further  support  for  the  existence 
of  a  site-dependence  of  wind-driven  noise.  We  now  turn  to  the  problem  of 
quantitatively  accounting  for  the  observed  differences  In  level  by  examining 
the  ocean  environment. 


Model  Description 

The  wind-noise  model  used  In  this  study  was  developed  jointly  by  NRL 
and  SACLANTCEN  and  has  been  described  In  detail  elsewhere  [11].  The  wind- 
noise  sources  are  modeled  as  monopoles  distributed  over  an  infinite  plane 
slightly  below  the  surface  and  parallel  to  It  (see  Fig.  4),  For  the 
purposes  of  this  study  the  sources  are  assumed  to  be  uncorrelated. 

Because  of  the  pressure  release  surface  the  sources  are  equivalent  to 
dipoles  at  the  surface.  The  water  column  has  depth-dependent  sound-speed 
and  overlies  a  bottom  consisting  of  a  sedimentary  layer  with  depth-dependent 
sound-speed  above  a  homogeous  fluid  or  solid  half-space.  Absorption  Is 
present  In  all  three  layers.  In  addition,  attenuation  due  to  roughness 
at  the  surface  and  at  the  Interface  between  the  water  column  and  the 
sedimentary  layer  Is  Included.  When  the  field  of  an  individual  noise 
source  Is  represented  In  terms  of  the  normal  modes  of  the  system  the 
total  intensity  at  depth  z  obtained  by  integrating  over  the  entire 
source  plane  Is 
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In  Eq.  (1)  I ( z , z ' )  is  the  intensity  at  depth  z  for  sources  located  at 
depth  z',  the  un's  are  the  mode  amplitude  functions,  the  kn's  are  the 
complex  modal  wavenumbers,  and  q  is  a  source  term. 


It  can  be  seen  from  Eq.  (1)  that  the  n  =  m  terms  dominate  the  sum. 
The  modal  wavenumber,  kn,  can  be  written  as 


kn  =  Kn  +  i  ctn 


(2) 


where  a  n  is  the  mode  attenuation  coefficient.  Neglecting  the  n  /  m  terms 
and  assuming  that  an  <<  Kn,  Eq.  (1)  reduces  to 


I(z.z') 


(3) 


Since  an  is  proportional  to  the  attenuation  in  the  system,  Eq.  (3) 
states  that  the  intensity  is  inversely  proportional  to  the  attenuation. 

Note  that  an  is  the  total  mode  attenuation  coefficient  and  is  the  sum  of 
several  mode-dependent  terms  representing  the  contributions  of  the  various 
attenuation  mechanisms.  In  most  cases  absorption  in  the  sediment  dominates 
but  loss  due  to  surface  roughness  and  absorption  in  the  water  column  can 
be  important  at  higher  frequencies. 


Results  and  Discussion 


In  order  to  use  the  wind-noise  model  described  above  to  calculate 
environmental  effects  on  the  measured  noise  level  it  Is  necessary  to  have 
certain  accompanying  environmental  data.  These  data  are  sound  speed  In 
the  water  column  as  a  function  of  depth  and  the  density,  sound  speed  an'1 
attentuatlon  coefficient  of  each  bottom  layer.  A  number  of  sufficiently 
complete  data  sets,  consisting  of  noise  spectrum  levels  and  environmental 
data,  have  been  collected  [1,5,12,13].  Figures  5  and  6  show  measured 
noise  levels  at  the  frequencies  500  Hz  and  1000  Hz,  respectively,  vs 
wlndspeed.  In  most  cases  the  wind  showed  short-term  fluctuations  over 
a  range  of  wlndspeeds  during  the  course  of  a  noise  level  measurement  so 
the  points  are  plotted  at  the  mean  wlndspeed.  Considerable  scatter  is 
observed  in  the  data.  Again  the  scatter  is  too  large  to  be  attributed  to 
Inaccuracies  In  the  measurements  and  uncertainties  in  the  wlndspeed.  The 
trend  of  the  measurements  indicates  that  the  measured  noise  level  varies 
approximately  as  the  square  of  the  wlndspeed  at  both  frequencies.  Thus 
the  data  is  consistent  with  earlier  observations  reported  In  the 
literature  [2]. 
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To  determine  the  environmental  effects  we  recast  Eq.  (1)  as  follows: 


Here  10  log  I  is  the  measured  noise  level.  The  second  term  on  the  right- 
hand  side  of  Eq.  (4)  can  be  calculated  from  the  model,  giving  the  source 
level,  10  log  q2.  Note  that  the  calculated  term  also  depends  on  windspeed 
through  the  attenuation  caused  by  surface  roughness.  The  source  level  will 
depend  on  windspeed  but  will  be  independent  on  the  other  environmental 
conditions.  The  success  of  the  explanation  of  the  scatter  of  the  levels 
in  Figs.  5  and  6  may  be  inferred  from  the  degree  to  which  the  independent 
determinations  of  source  level  produce  points  which  lie  on  a  common  curve 
when  plotted  vs  windspeed. 

The  source  levels  for  the  data  of  Figs.  5  and  6  are  plotted  in 
Figs.  7  and  8.  It  is  immediately  evident  that  the  scatter  has  been 
reduced  considerably  by  the  removal  of  the  environmental  effects.  We 
note  that  because  the  calculated  absorption  coefficients  depend  on  surface 
roughness,  and  hence  on  windspeed,  the  calculated  source  level  depends  on 
a  different  power  of  the  windspeed  than  does  the  noise  level.  At  500  Hz 
the  source  level  increases  with  the  fifth  power  of  the  windspeed;  at 
1000  Hz  the  source  level  depends  on  the  sixth  power  of  the  windspeed. 


Conclusion 

Measurements  of  wind-driven  ambient  noise  from  a  number  of  sites  under 
similar  windspeeds  show  a  site-to-site  variability  of  noise  level  which 
can  not  be  readily  attributed  to  errors  in  acoustic  and  environmental 
measurements  alone.  Some  of  the  measurement  sets  were  accompanied  by 
sufficient  water  column  and  ocean  bottom  information  to  permit  estimates  of 
the  effects  of  the  acoustic  environment  on  noise  level.  These  calculated 
site-dependent  terms  account  for  the  observed  differences  among  the  sites. 
We  conclude  that  a  more  accurate  prediclton  of  wind-noise  level  can  be 
made  If  water  column  and  ocean  bottom  information  is  available  than  can  be 
obtained  from  "universal"  ambient  noise-vs-windspeed  curves. 
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FIG.  1  MEASURED  AMBIENT  NOISE  SPECTRUM  LEVELS  UNDER 
SEA  STATE  2  CONDITIONS 


FIG.  2  MEASURED  AMBIENT  NOISE  SPECTRUM  LEVELS  UNDER 
SEA  STATE  3  CONDITIONS 
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ABSTRACT 


Low  frequency  noise  in  most  ocean  areas  is  dominated  by  shipping  or 
industrial  noise.  The  composite  noise  field  is  the  result  of  radiated 
energy  from  each  discrete  noise  source  propagating  over  relatively  large 
distances  to  the  receiving  site.  The  ocean  bottom  in  general,  and 
topographic  features  In  particular,  result  in  vertically  directional  noise 
fields  in  addition  to  the  expected  azimuthal  directionality.  The  noise 
field  near  the  surface  displays  broadband  coherent  interaction  between  the 
direct  and  surface  image  energy,  the  details  of  which  can  be  traced  in  part 
to  the  intervening  topography.  At  deeper  depths  the  level  versus  depth 
curves  also  display  the  impact  of  the  topographic  blockage.  A  measured 
data  set  displaying  near  surface  to  near  bottom  low  frequency  noise  levels 
is  analyzed  with  respect  to  range  dependent  environmental  factors  and 
shipping  densities. 


INTRODUCTION 


Recently,  the  depth  dependence  of  the  low  frequency  ocean  ambient  noise 
field  has  been  a  subject  of  some  Interest1.  Depth  profiles  of  ambient 
noise  were  obtained  from  two  AUTOBUOY  dives  to  depths  of  2800  m  which  were 
conducted  at  widely  separated  sites  under  moderate  sea  conditions  as  part 
of  the  NORLANT  *72  experiment*.  In  both  cases  no  ship  traffic  was  observed 
In  the  general  area  of  the  dive*  The  deployment  vessel,  USNS  SANDS 
(T-AGOR-6),  was  at  quiet  ship  condition  during  most  of  each  dive  (using 
only  Its  gas  turbine  on  the  02  level  to  power  the  ship's  facilities)  at  a 
distance  of  5  nm  from  the  launch  point  for  the  first  dive  on  22  July  and 
2  nm  from  the  launchpolnt  for  the  second  dive  on  31  July  1972.  Aural  and 
graphical  monitoring  of  the  recorded  data  Indicated  that  no  obvious  nearby 
ship  noise  was  present.  There  were,  however,  waterborne  transients 
occurring  at  repetition  Intervals  of  approximately  10,  18,  and  24  seconds. 
These  transient?  had  low  frequency  characteristics  and  time  durations  which 
Indicate  that  their  sources  were  some  distance  away  and  these  have  been 
Identified  as  sparkers  and  air-guns  being  used  for  sub-bottom  profiling  In 
the  Labrador  Basin  and  over  the  Grand  Banks.  The  effect  of  these  impulses 
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on  the  ambient  level  was  eliminated  by  processing  only  in  the  interval 
between  pulses.  Thus,  the  data  reported  here  represent  only  ambient  noise 
generated  by  the  sea  surface  and  distant  shipping. 

The  resulting  data  display  different  noise  level  versus  depth  relationships 
at  the  two  sites.  Furthermore,  at  one  site  data  were  acquired  to  very 
shallow  depths  and  near-surface  interference  was  observed  in  the  resulting 
noise  versus  depth  patterns.  Two  noise  models,  DANES3  and  FANM4,  were 
tested  against  the  data  to  determine  their  utility  in  describing  the 
observations.  In  addition,  the  range  dependent  Parabolic  Equation5 
transmission  loss  model  was  used  to  assess  the  effects  of  seafloor 
bathymetry  on  noise  versus  depth  resulting  from  distant  shipping. 


AUTOBUOY  DESCRIPTION 

AllTOBUOY 6  is  a  self-contained,  untethered,  programmable,  free-diving, 
acoustic  data  gathering  system  which  is  capable  of  maximum  depths  to 
6,000  m.  It  is  generally  programmed  to  descend  quickly  to  a  selected  maximum 
depth,  to  release  a  descent  weight  in  order  to  achieve  neutral  buoyancy  and 
then  to  hover  sequentially  at  four  shallower  depths.  At  each  depth,  several 
types  of  data  are  recorded  on  the  seven  channel  tape  recorder  located  in 
the  instrumentation  pressure  vessel.  These  data  include  calibration  signals 
and  the  acoustic  signals  from  each  of  two  closely  spaced  (*W0  m) 
hydrophones  together  with  depth,  time,  ascent  velocity  and  temperature 
signals.  Depth  is  controlled  during  the  dive  by  adjusting  the  buoyancy 

around  a  neutral  value  in  response  to  depth  and  velocity  error  signals. 
This  is  accomplished  by  the  alternate  valving  to  sea  of  light  and  heavy 
fluids.  While  hovering,  AUTOBUOY  Is  free  to  move  with  the  water  mass. 
Thus,  no  local  noise  is  generated  by  water  flow  past  the  hydrophones.  The 
fluid  valving  system  is  operated  by  solenoids  which  open  valves  such  that 
fluid  flows  out  by  gravitational  force.  No  noise  Is  detectable  from  this 
operation.  The  depth  and  velocity  signals  also  determine  start  times  of  an 
analog  tape  recorder  which  is  preprogramed  for  a  specific  sample  size 
(typically  20  minutes).  When  recording  is  completed,  the  tape  is  stopped 

and  AUTOBUOY  ascends  to  the  next  programmed  depth*  When  the  last  depth 
station  is  completed,  the  system  rises  to  the  surface  where  an  attached 
recovery  system,  consisting  of  an  rf  beacon  and  a  flashing  light,  is 
activated. 


DIVE  DESCRIPTIONS 

The  dive  sites  for  NORLANT  *72*  are  shown  in  Figure  1.  Dive  l  occurred  in 
the  northeastern  section  of  the  Labrador  Sea  between  Greenland  and  the 
Reykjanes  Ridge.  Oive  2  occurred  in  the  deeper  area  of  the  basin  350  m  due 
south  of  the  southern  tip  of  Greenland.  For  each  dive,  AUTOBUOY  was 
programmed  to  descend  to  2900  m  and  then  to  hover  at  2300,  1700,  1200  and 
600  m  depths  while  recording  20  minutes  of  ambient  noise  at  each  depth.  At 
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the  shallowest  depth  station,  recording  continued  until  the  tape  ran  out. 
Figures  2a  and  3a  pictorial ly  display  the  actual  dive  patterns  for  each 
event.  During  the  first  dive,  hovering  was  achieved  only  at  the  two 
deepest  depths  while  during  the  second  dive  hovering  occurred  at  all  but 
the  shallowest  depth.  When  hovering  did  not  occur,  ascent  rates  were 
approximately  1.2  m/sec  (about  0.25  kt).  At  this  slow  speed,  no  observable 
mechanical  or  flow  noise  is  introduced  into  the  hydrophones.  Failure  to 
hover  during  a  depth  station  is  usually,  the  result  of  valving  too  much 
heavy  fluid  at  the  deeper  depths  and  not  having  sufficient  reserve  for  the 
shallower  depths.  Also  indicated  on  Figures  2a  and  3a  are  the  ship  and 
wind  speed  conditions  during  the  dives. 


SPECTRAL  ESTIMATES  AND  DEPTH  DEPENDENT  RESULTS 

Figures  2b  and  3b  show  the  resulting  spectra  from  the  hydrophone  acoustic 
outputs  for  the  two  dives.  Two  hydrophones,  vertically  separated  by  *40  m, 
yielded  virtually  identical  results.  All  data  were  processed  in  1/3-octave 
bands  centered  at  the  indicated  frequencies.  Data  were  processed  twice: 
once  on  board  the  SANDS  through  log-amplifiers  onto  a  strip  chart  recorder 
and  later  ashore  using  a  digital  automatic  data  processing  system  (ADP). 
In  both  cases  all  transients  (soismlcs,  airguns,  etc.,)  were  eliminated 
from  the  data*  Calibration  signals  over  the  entire  band  were  injected  at 
the  hydrophone  preamplifiers  in  the  laboratory  prior  to  deployment  to 
produce  a  calibration  tape.  Also,  a  broadband  calibration  signal  was 
automatically  Injected  at  the  recorder  inputs  during  the  dive  at  the  start 
of  each  20  minute  segment.  The  shipboard  and  ADP  results  were  essentially 
Identical.  The  principal  difference  between  the  spectra  for  the  two  dives 
Is  their  depth  dependence.  At  the  deepest  depth,  the  spectrum  of  the  22 
July  data  has  a  variation  in  level  with  frequency  of  only  7  dB  over  the 
range  25  Hz  to  800  Hz.  As  depth  decreases  the  spectra  slowly  approach  a 
more  familiar  shape.  The  leveling  off  of  the  spectra  between  200-1000  Hz 
Is  a  characteristic  noted  by  others7*1  and  Is  attributed  to  the  overlap 
between  ship  generated  noise  and  wind  generated  noise.  In  order  for  this  to 
be  the  explanation  of  the  22  duly  spectra  depth  dependence,  a  significant 
reduction  In  ship  noise  contribution  would  be  required  for  the  deep  data* 

Figures  2c  and  3c  show  the  same  data  plotted  to  display  the  noise  depth 
dependence  at  the  two  sites*  Each  data  set  represents  measurements  over 
approximately  6  hours  with  depths  measured  sequentially,  rather  than 
concurrently.  Nevertheless,  In  the  absence  of  Individual  nearby  noise 
sources,  the  data  should  be  a  valid  representation  of  persistent  background 
noise  as  a  function  of  depth  for  the  measured  locations  and  for  the 
conditions  that  exist  In  the  sumer  months.  Horizontal  noise  pattern 
measurements  in  the  northeast  Pacific  by  R.  Hagstaff 4 ,  for  Instance, 
demonstrate  repeatable  results  from  year  to  year  for  the  same  season  In  the 
shipping  noise  dominated  portion  of  the  spectrum* 

Several  Interesting  features  of  the  noise  depth  dependent j  can  be  noted 
from  the  data  displayed  in  Figures  2  and  3.  The  most  striking  is  the 
strong  depth  dependence  of  the  low  frequency  noise  below  critical  depth 
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(see  the  sound  speed  profile  on  Figure  2c)  at  the  location  of  Dive  1.  This 
low  frequency  noise  reduction  below  critical  depth  is  consistent  with 
observations  in  other  ocean  basins1.  It  results  in  part  because  most  cf 

the  energy  from  the  distant  noise  sources  is  carried  by  low  order  normal 

modes  of  the  water  column  and  for  these  modes  the  amplitude  decays 
exponentially  with  depth  below  the  critical  depth.  The  rate  of  decay  of 
the  noise  level  with  depth  below  critical  is  determined  by  the  energy 

distribution  among  the  normal  modes.  This  distribution  is  governed  by 
several  factors  including:  the  distribution  of  energy  among  the  normal 

modes  of  the  water  column  of  the  source  region(s),  changes  in  the 
distribution  of  noise  energy  among  modes  as  the  energy  travels  from  sources 
to  receiver  through  a  water  column  with  horizontally  changing  sound  speed 
versus  depth  structure  (and  thus  changing  normal  mode  eigenfunctions)*, 
^varying  modal  attenuation  along  the  propagation  path  resulting  from  several 
factors  including  bottom  interaction  of  the  modes.  We  will  return  to  a 
discussion  of  some  of  these  factors  below  when  comparisons  are  made  with 
numeri cal  model  noise  computations. 

Other  notable  features  of  the  noise  depth  profiles  include  the  virtually 
constant  level  with  depth  for  frequencies  above  300  Hz  at  the  site  of  Dive 
1  (Figure  2c)  and  the  much  weaker  depth  dependence  throughout  the  water 
column  at  the  site  of  Dive  2  for  all  frequencies  as  compared  with  the  low 
frequencies  at  the  site  of  Dive  1. 


DEPTH  DEPENDENCE  ANALYSIS 

Two  computerized  ambient  noise  numerical  models  were  used  to  produce 
simulations  of  the  noise  fields  for  the  sites  of  Dives  1  and  2.  The  models 
used  are  DANES3  and  FANM 4  .  Transmission  loss  models  are  useful  for 

analyzing  some  features  of  the  noise  field  depth  dependence.  For  this 
purpose,  simulations  were  run  using  the  Parabolic  Equation  (PE)  computer 
numerical  transmission  loss  model  of  Tappert5.  In  this  section,  the  models 
and  model  Inputs  are  briefly  described.  Subsequently,  model  simulations 
are  compared  with  the  NORLANT  *72  AUTOBUOY  data  and  the  comparisons  are 
discussed. 

The  FANM  ambient  noise  model  is  essentially  a  computer  Implementation  of 
the  Talham  model10.  It  treats  water  column  parameters  as  range  Independent 
(single  sound  speed  profile),  but  allows  bottom  depth  and  bottom  less  to 
vary  with  range}  thus,  energy  from  distant  ships  can  be  cut  off  by 
intervening  bathymetric  features.  DANES  utilizes  a  fully  range  dependent 
environment  with  propagation  described  bv  a  range  averaged  transmission 
loss  for  15°  resolution  tracks.  Data  base  for  sound  velocity  profiles  (by 
season  and  area),  bottom  loss  provinces,  bottom  depths  and  shipping 
distribution  (by  month)  are  accessed  by  DANES  as  indicated  by  input 
information  (site  location,  time  of  year,  weather).  Sources  are 

characterized  as  a  function  of  frequency,  by  a  source  at  depth  of  6  m; 
source  spectra  are  modified  ROSS-ALVAREZ  and  ship  distribution  Was  obtained 
from  the  RMS  data  base.  Similar  data  bases  for  bottom  depth,  bottom  loss 
and  shipping  distribution  were  accessed  by  FANM.  Figure  4  summarizes  the 
ship  count  data  base.  This  data  base  contains  ship  counts  in  14  by  1° 
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( tat. /Long. )  squares.  The  data  are  shown  in  Figure  4  using  revised 
divisions  which  more  clearly  indicate  the  relationship  of  ship  counts,  the 
bottom  bathymetry  and  the  measurement  sites-  Aerial  reconnaissance  during 
NORLANT  '72 2  identified  only  5  ships  in  a  5°  square  area  near  the 
measurement  sites.  No  ships  were  observed  from  the  SANDS  during  the 
measurements  at  the  site  of  Dive  1.  The  higher  density  shipping  to  the 
south  in  the  region  shown  in  Figure  4  is  expected  because  of  shipping  lanes 
in  this  area. 


In  Figure  5,  the  NORLANT  data  at  25  Hz  and  50  Hz  are  compared  with  the 
simulated  ambient  noise  depth  profile  from  the  DANES  model.  DANES  was 
executed  with  a  sound  speed  profile  for  the  Dive  2  site  which  was  measured 
during  NORLANT  *72  and  with  archival  sound  speed  profiles  in  other  portions 
of  the  area.  Other  data  bases  needed  (ship  distribution,  bottom  type, 
etc.),  were  accessed  as  appropriate  to  the  model.  Although  the  DANES 
simulated  noise  levels  are  reasonably  close  to  the  measurements  at  shallow 
depths,  the  details  of  the  depth  dependence  differ  between  data  and  model. 
FANM  simulated  results  (with  similar  input  parameters)  at  25  Hz  are  also 
shown  on  Figure  5.  They  agree  with  the  data  less  than  the  DANES  results  at 
shallow  depths  and  exhibit  a  greater  change  in  level  with  depth  than  either 
the  data  or  the  DANES  simulations.  It  has  been  shown11  that  range  averaged 
transmission  loss  simulations  can  be  valuable  in  analyzing  the  depth 
dependence  of  ambient  noise  fields.  Such  simulations  were  made  using  a 
Parabolic  Equation  5  computer  model  with  a  noise  source  depth  assumed  to  be 
6  m  and  with  range  varying  input  parameters  for  bottom  bathymetry  and  sound 
speed  profiles  as  shown  in  Fig  6b.  These  input  parameters  are 
representative  of  the  environment  to  the  south  of  the  Dive  2  site,  in  the 
direction  of  most  of  the  noise  sources  (ships).  Results  of  these 
transmission  loss  simulations,  averaged  over  a  source-to-receiver  range 
Interval  of  500  to  600  nm  are  also  shown  In  Figure  5.  The  simulations 
shown  are  from  a  low  loss  bottom  (0.2  dB/nm  sediment  attenuation).  It  is 
significant  that,  while  they  differ  In  terns  of  the  details  of  the  noise 
level  versus  depth  profile,  in  general,  all  of  the  simulation  results  agree 
with  the  data  in  describing  only  a  small  (at  most)  dependence  of  noise 
level  on  depth  over  the  sampled  depth  interval. 


A  more  Interesting  set  of  comparisons  is  shown  in  Figure  7,  wherein  NORLANT 
results  for  the  Give  1  site  are  shown  in  comparison  with  DANES  noise  and  PE 
average  transmission  loss  depth  profiles  at  25  Hz  and  50  Hz  and  with  FANM 
noise  depth  profiles  at  25  Hz.  There  is  clearly  significant  disagreement 
between  the  measurements  and  all  of  the  simulations.  The  monotonic 
decrease  of  measured  noise  levels  with  depth  is  significantly  greater  than 
either  noise  model  simulation  and  the  character  of  the  average  transmission 
loss  profile  disagrees  with  the  measurements  especially  in  the  sharp 
Increase  in  loss  (which  would  produce  a  similar  decrease  in  noise  level) 
below  about  2000  m.  this  situation  has  been  analyzed  further  by  a 
parametric  .  study  of  the  Inpact  of:  a)  the  range  dependent  nature 
(horizontal  variation)  of  the  sound  speed  profile  in  this  area;  and  b)  the 
topographic  blockage  introduced  by  the  Reyk  janes  Ridge  which  Is  between  the 
receiving  site  for  Dive  1  and  most  of  the  ship  noise  sources  (Figure  4).  A 
Parabolic  Equation  model  was  used  In  this  study;  the  basic  assumption  Is 
that  the  major  contributor  to  the  noise  field  at  the  Dive  i  site  is  distant 
shipping  as  local  shipping  was  non-existent.  The  majority  of  the  distant 
shipping  lies  on  the  far  side  of  the  Reykjanes  Ridge  at  distances  greater 
than  300  nm.  The  Ridge  is  of  varying  depth,  having  a  somewhat  average 
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monotonic  decline  from  the  surface  at  Iceland  in  the  northeast  to  the 
Labrador  Basin  in  the  southwest,  ending  approximately  due  south  of  the 
southern  tip  of  Greenland.  The  PE  model  runs  used  a  source  at  6  m  depth 
and  were  made  to  a  range  of  600  nm  from  the  receiver.  The  water  depth  at 
the  source  was  chosen  as  4500  m  rising  to  3300  m  at  the  start  of  the 
Reykjanes  Ridge,  which  was  chosen  as  the  midway  point.  The  base  of  the 
Ridge  was  assumed  to  be  200  nm  wide  with  the  Ridge  plateau  rising  to 
selected  depths  below  the  surface.  The  receiver  was  located  100  nm  beyond 
the  far  side  of  the  Ridge  base  (see  Figure  6a).  All  propagation  events 
used  here  assumed  a  flat  pressure  release  surface  and  bottom  character¬ 
istics  as  follows:  sediment  sound  speed  gradient  of  0.1  m  along  the  entire 
track;  0.2  dB/nm  attenuation  in  the  sediment  from  the  source  to  the  Ridge 
and  a  fully  reflecting  bottom  at  the  Ridge  and  beyond.  While  this 
characterization  is  not  realistic,  runs  performed  using  a  highly  absorbing 
bottom  along  the  entire  track  differed  in  transmission  loss  versus  depth 
characterizations  by  less  than  3  dB  overall  with  near  bottom  exponential 
roll-off  starting  at  about  100  m  shallower  depth.  The  sound  speed  was 
assumed  to  be  range  dependent  with  the  profile  at  the  receiving  site  chosen 
to  be  that  measured  during  the  exercise  and  all  others  chosen  from  archival 
data  for  midsunmer  in  the  direction  from  the  receiver  to  the  southeast. 


Figure  8  displays  transmission  loss  versus  depth  profiles  averaged  over  the 
500  to  600  nm  source-to-receiver  range  interval.  Examining  the  depth 
dependence  for  such  a  range  averaged  transmission  loss  is  equivalent  to 
examining  the  noise  field  depth  dependence  for  a  uniform  distribution  of 
noise  sources  over  the  500  to  600  nm  range  interval.  The  family  of  curves 
shown  was  obtained  by  repeating  the  simulation  for  a  variety  of  Ridge  crest 
depths  (shown  in  the  circle  by  each  curve).  Also  shown  on  Figure  8  is  the 
range  averaged  transmission  loss  depth  profile  resulting  from  a  simulation 
with  a  single  bottom  depth  (3300  m)  and  a  single  sound  speed  profile  (that 
measured  at  Site  1  during  the  acoustic  measurements).  The  first  thing  to 
note  from  Figure  8  is  the  significant  increase  In  the  portion  of  the  water 
column  depth  with  essentially  constant  average  transmission  loss  which 
results  from  range  dependent  simulations  (for  the  case  of  no  bathymetric 
Intrusion)  as  opposed  to  those  from  constant  environment  simulations. 
Thus,  In  the  absence  of  the  Reykjanes  Ridge,  one  might  anticipate  a  noise 
field  of  very  nearly  constant  level  with  depth  down  to  depths  on  the  order 
of  3  km.  As  one  would  anticipate,  the  attenuation  of  energy  carried  by 
higher  order  modes  resulting  from  the  bathymetric  intrusion  by  the  Ridge 
results  In  confinement  of  the  noise  energy  to  shallower  portions  of  the 
sound  channel.  Indeed,  as  the  top  of  the  (simulated)  Intrusion  becomes 
shallower,  the  resulting  distribution  of  energy  (Figure  8)  more  and  more 
nearly  represents  a  single  low  order  mode  of  the  sound  channel  at  the 
receiver  site.  Although  none  of  the  profiles  of  average  transmission  loss 
with  depth  in  Figure  8  exactly  matches  the  Site  1  data,  the  profiles  do 
make  the  case  that  the  bathymetric  Intrusion  represented  by  the  Ridge  will 
significantly  change  the  depth  distribution  of  energy  for  noise  sources 
across  the  Ridge  from  the  receiver.  Specifically;  the  Ridge  combined  with 
the  local  sound  speed  profile  variation  will  tend  to  confine  more  energy  to 
the  upper  portion  of  the  water  column.  This  tendency,  the  varying  depth  of 
the  Ridge  along  Its  path  between  sources  and  receiving  site,  and 
contributions  from  sources  not  blocked  by  the  Ridge  (including  at  least  a 
few  at  shorter  range)  conspire  to  produce  the  observed  noise  level  versus 
depth  profile. 
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AMBIENT  NOISE  BOUNDARY  EFFECTS  NEAR  THE  SEA  SURFACE 

During  Dive  1  (Figure  2)  AUTOBUOY  made  a  slow  ascent  from  approximately  360 
m  to  20  m  below  the  surface  (as  indicated  by  its  depth  gauge  signal)  where 
it  hovered  for  several  minutes  until  it  received  the  command  to  surface. 
During  this  ascent,  as  AUTOBUOY  neared  the  surface,  a  radical  change  in  the 
ambient  noise  level  at  low  frequency  was  observed  (see  Figure  9).  The 
level  change  during  the  last  340  m  of  depth  change  ranged  from  11.5  dB  at 
25  Hz  to  0.5  dB  at  200  Hz  for  the  upper  hydrophone  and  from  7  dB  at  25  Hz 
to  0.5  dB  at  125  Hz  for  the  deeper  hydrophone.  These  two  hydrophones,  as 
mentioned  earlier,  are  approximately  10  m  apart  so  that  when  the  upper 
hydrophone  is  at  20  m,  the  lower  hydrophone  is  at  30  m  depth.  Assuming  a 
directional  noise  field,  these  observations  can  be  interpreted  as  a 
manifestation  of  the  Lloyd  Mirror  Effect  which  basically  states  that  the 
resultant  output  of  a  hydrophone  near  a  boundary  can  be  described  as  the 
algebraic  sum  of  the  output  of  the  hydrophone  and  its  image  (the  reflection 
of  the  hydrophone  in  the  boundary  taking  into  account  the  change  in 
acoustic  impedance  at  the  boundary).  For  the  case  of  a  pressure  release 
interface,  the  boundary  effect  results  in  a  180°  phase  change  of  the 
acoustic  signal  when  reflected  from  the  surface.  Figure  10  displays 
effective  vertical  dipole  beam  patterns  for  a  receiver  at  20  m  below  the 
surface  for  two  different  frequencies j  viz.,  31  and  100  Hz.  If  the  receiver 
were  60  m  below  the  surface,  then  the  31  Hz  pattern  would  be  nearly 
identical  to  the  100  Hz  pattern  for  a  depth  of  20  m.  If  the  noise  field 
were  Isotropic,  integration  over  these  two  patterns  for  a  single  frequency 
but  different  receiver  depths  (say  for  31  Hz  at  20  and  60  m)  would  yield 
approximately  the  same  value  and  therefore  no  significant  depth  effect 
would  be  noted.  However,  If  noise  is  vertically  directional  as  would  be 
expected  if  the  dominant  sources  are  distant  ships,  then  a  significant 
decrease  can  be  expected  as  the  receiver  approaches  a  boundary. 
Furthermore,  the  shape  of  the  receiver  dipole  pattern  Is  not  changed  for 
the  horizontal  null  If  the  noise  Is  band  limited  rather  than  single 
frequency.  Note,  In  Figure  10,  that  the  dashed  curve  for  1/3-octave  noise 
Is  Identical  to  that  for  a  single  frequency  at  the  band  center  for  the 
horizontal  null  and  becomes  progressively  modified  for  subsequent  nulls. 


ANALYSIS  OF  NEAR  SURFACE  NOISE  DEPTH  DEPENDENCE 

The  pattern  function  for  the  dipole  receiver,  described  above.  Is  given  by: 

I  ■  4I0  s1na(2ir  |  sin  0) 

where  the  symbols  are  defined  In  Figure  11.  The  FANM  noise  model  computes 
the  noise  eneroy  arriving  at  a  receiver  In  several  vertical  arrival  angle 
bins.  A  special  version  of  FANM  was  generated  to  apply  the  above  pattern 
function  separately  to  each  vertical  arrival  angle  bln  (using  the  angle 
with  the  horizontal  of  the  center  of  the  bln  as  6)  and  then  sum  the 
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resulting  energy  in  the  bins  for  a  total  (omnidirectional)  noise  level  . 
This  FANM  version,  producing  what  is  called  a  coherent  sum,  was  used  to 
simulate  the  noise  versus  depth  profiles  for  the  upper  hydrophone  data  of 
Figure  9„  In  these  simulations,  sharp  decreases  in  level  occurred  at  low 
frequency  as  the  surface  was  approached.  This  is  in  contrast  to  a  normal 
(incoherent)  FANM  calculation  for  which  noise  level  either  remains  constant 
or  actually  increases  slightly  as  a  function  of  decreasing  depth  near  the 
sea  surface.  In  Figure  12,  calculations  with  the  new  (coherent)  FANM  model 
are  compared  with  the  NORLANT  '72  near  surface  data.  The  depth  dependence 
of  the  data  and  the  FANM  simulations  is  in  excellent  agreement. 

Recall  that  the  FANM  model  cuts  off  noise  energy  which  is  blocked  by 
bathymetric  features.  In  the  earlier  discussion  of  the  noise  depth 
profile,  it  was  postulated  that  the  noise  field  at  Site  1  resulted  from  a 
combination  of  noise  energy,  much  of  which  had  undergone  varying  degrees  of 
bathymetric  blockage  by  the  Reykjanes  Ridge.  This  was  based  in  part  on  the 
information  shown  in  Figure  8.  In  addition  to  a  redistribution  of  energy 
in  the  deep  portions  of  the  water  column,  the  bathymetric  blockage  modifies 
the  near  surface  interference  pattern.  Figure  13  shows  an  expanded  display 
of  the  near  surface  portion  of  the  parametric,  PE  transmission  loss  study 
of  Figure  8.  Modification  of  the  details  of  the  near  surface  interference 
pattern  by  the  bathymetric  feature  is  evident.  One  persistent  feature  of 
the  near  surface  interference  pattern  is  the  final  monotonic  decrease  in 
noise  energy  as  the  receiver  approaches  the  surface.  This  is  a  result  of 
the  pressure  release  surface  and  resulting  exponential  decay  of  all  modes 
near  the  surface12.  Reduction  of  the  interference  with  increasing  Ridge 
height  (decreasing  N)  is  consistent  with  the  earlier  description  of  higher 
order  mode  attenuation  as  the  height  of  the  bathymetric  intrusion 
increases.  Increased  dominance  of  lower  order  modes  can  be  interpreted  as 
a  reduction  of  ray  angle,  with  the  horizontal,  for  the  interfering  energy. 
From  the  equation  of  Figure  11,  as  the  height  of  the  bathymetric  intrusion 
increases,  the  attendant  decrease  of  should  result  in  a  shift  to  greater 
depth  of  various  features  of  the  interference  pattern,  such  as  the  depth  of 
occurrence  of  peaks  and  nulls.  Oust  this  behavior  is  exhibited  by  the 
results  shown  in  Figure  13.  For  example,  the  first  minimum  in  transmission 
loss  (maximum  of  noise  energy)  occurs  at  its  shallowest  depth  with  no 
bathymetric  interference  and  monotonically  increases  in  depth  as  the  Ridge 
height  increases.  The  near  surface  features  of  the  NORLANT  '72  data  are 
consistent  with  an  Interpretation  of  the  field  as  a  result  of  noise  energy 
which  has  undergone  varying  degrees  of  bathymetric  blockage  by  the 
Reykjanes  Ridge.  Note  that  the  depth  of  the  second  (and  higher  order) 
nulls  in  the  patterns  of  Figure  13  would  be  filled  in  by  the  1/3-octave 
processing  of  the  NORLANT  data,  although  the  first,  near  surface  null  would 
be  unaffected  (recall  Figure  10).  This  and  the  reduction  of  the 

Interference  null  by  sea  surface  roughness  are  probable  explanations  for 
the  absence  of  second  and  higher  order  nulls  in  the  measured  surface  noise 
field. 


CONCLUSION 

In  the  NORLANT  *72  data,  both  the  near  surface  Interference  pattern  and  the 
depth  profile  of  the  noise  field  for  the  deeper  portion  of  the  water  column 
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have  been  shown  to  be  consistent  with  a  postulated  description  of  the  npjse 
field  as  a  combination  of  noise  energy  fields  from  distant  sources  which 

have  experienced  varying  degrees  of  blockage  by  the  Reyk janes  Ridge. 
Numerically  simulated  range  averaged  transmission  loss  has  been  shown  to  be 
a  useful  tool  for  analyzing  the  depth  dependence  of  ambient  noise  fields. 
Such  simulations  have  illustrated  the  impact  of  bathymetric  blockage  on 
both  the  noise  versus  depth  profile  in  the  deep  portion  of  the  water  column 
and  on  the  near  surface  interference  pattern.  Although  they  exhibited 
reasonable  success  in  simulating  the  noise  depth  profile  in  the  simpler 
region  of  the  NORLANT  *72  Dive  2,  neither  the  DANES  nor  the  FANM  model  was 
successful  in  reproducing  the  noise  depth  profile  of  more  complicated 
origin  at  the  Dive  1  site. 
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THE  SEAMOUNT  AS  A  NOISE  BARRIER 


by 

Herman  Medwin,  Emily  Childs*  and  Edgar  A.  Jordan** 

Physics  Department,  Naval  Postgraduate  School,  Monterey,  CA  93940,  U.S.A. 


ABSTRACT 

When  the  path  from  a  localized  noise  source  is  intercepted  by  a  seamount  the 
shielded  noise  consists  of  a  diffracted  component  and,  if  the  seamount  is 
close  to  a  sufficiently  smooth  ocean  surface,  a  multiply-reflected  component 
as  well.  The  magnitude  and  frequency  filtering  of  the  blockage  caused  by 
the  seamount  is  here  determined  by  using  the  impulse  solution  of  diffraction 
by  a  rigid  wedge  as  a  module  for  the  computer  calculation  of  multiple 
diffraction  at  gross  changes  of  slope  along  the  seamount  and  at  Its  crest 
line.  These  computer  calculations  are  validated  by  comparison  with 
laboratory  measurements  using  a  scale  model  of  DICKENS  Seamount  which  has 
previously  been  shown  to  be  an  effective  technique  for  the  calculation  of 
topographical  diffraction  at  sea. 


BACKGROUND 

The  ambient  noise  measured  at  any  location  depends  on  the  noise  source  and 
the  transmission  path.  A  factor  that  may  significantly  affect  the 
propagation  Toss  is  the  Interception  of  the  noise  path  by  a  seamount; 
shadowing  losses  of  15  to  20  dB  are  commonly  measured. 

One  successful  technique  for  the  prediction  of  shadowing  losses  by  a 
seamount  is  the  use  of  a  2H  x  2M  laboratory  scale  model1  (Figs.  1,2) j  the 
scaled  frequency ^dependent  diffraction  losses  at  the  seamount  are  added  to 
the  losses  in  the  body  of  «.ne  sea  which  are  calculated  by  a  computer  model 
such  as  FACT.  There  is  good  agreement  with  ocean  data  (Figs.  3,4).  The 
seamount  is  assumed  to  be  rigid;  no  other  assumptions  art  needed.  The 
increase  of  shadowing  loss  with  frequency1/2  has  since  been  observed  in 
several  ocean  experiments  by  different  laboratories*.  Whether  this 
diffraction  component  is  dominant  depends  on  the  distance  between  the 
seamount  and  the  sea  surface  as  well  as  the  roughness  of  the  sea  surface^. 

‘Ocean  Acoustics  Associates 
“IT,  U.S.  Navy 


SACUNTCEN  CP-32  24-1 


MEDWIN  et  al:  Seamount  as  a  noise  barrier 


In  the  example  studied,  DICKENS  Seamount  with  a  slope  of  about  14° 
culminating  in  a  crest  approximately  500M  from  the  surface,  a  35  knot  wind 
resulted  in  an  ocean  surface  of  rms  height  1.6M  calculated  from  the 
Pierson-Moskowitz  wind  wave  spectrum.  In  this  case,  the  model  studies  show 
that  not  only  was  the  diffraction  component  the  f i rst,  it  was  also  the 
largest  component;  multiple  reflections  between  the  seamount  and  the  wind 
blown  sea  produced  a  much  weaker  signal.  Note  signal  after  203  ps  in  Fig.  5 

The  physical  basis  for  the  interpretation  of  the  diffraction  is  the  fact 
that  the  diffracted  energy  comes  principally  from  a  region  near  the  sound 
ray  track.  The  Crestline  of  the  seamount  at  the  sound  track  may  be 
approximately  most  simply  by  the  apex  of  a  flat  wedge;  the  rigorous  closed 
form  impulse  solution  of  a  rigid  wedge  is  given  by4  eqs.  (l)-(3), 

/>(/-)  =  (-.Spc/4naB){p}  [rr0sinhv  exp(-nv/Pw) , 
where 

;  _ Sin[(fl/r,„K>nft±fln)l  _ 

1  -2exp(-n7tfJeosl(fl/GJ(!i±tf  ±i’o)]  +  exp(-27r}'/0J  ’ 


(1) 

(2) 


Y 


-arc  cosh 


(3) 


The  term  (7r+0+0o)  is  written  for  simplicity;  the  curly  bracket 
consists  of  the  sum  of  four  terms  obtained  by  using  the  four  possible 
combinations  of  the  angles,  +0  and  +0O. 

Source  coordinates  are  (r0,  Q0,  0),  receiver  coordinates  are 

(r,  0,  Z).  The  angle  of  the  wedge  measured  in  the  fluid  region  Is  called 

0W. 

The  sound  diverges  spherically  from  the  point  source  and  cylindrical ly  from 
the  crestal  line  scatterer.  The  earliest,  strongest,  pressure  arrivals  vary 
as  (time) -1/2,  and  the  frequency  dependence  of  the  scatter  goes  as 
f  1/2.  For  these  reasons,  the  appropriate  model  scaling  parameter  Is  the 
far  field  diffraction  strength,  DS,  given  by 

DS  ■  20  logic  C(PD/P0)(r0/R0)(rA)l/2  ]  (4) 

where 

Pg  *  diffracted  pressure 

Pc  a  reference  pressure  at  range  R0 

R0  a  reference  range  (-1  M) 

A  =  wave  length 

r0  =  range,  source  to  diffracting  edge 
r  =  range,  diffracting  edge  to  receiver 

The  far  field  diffraction  strength  is  a  function  of  wedge  angle,  angle  of 
Incidence,  elevation  and  azimuthal  scatter  angles. 
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The  concept  of  diffraction  strength  is  used  throughout  this  work.  To 
convert  to  diffraction  loss  DL  in  dB  one  uses 

DL  =  20  1og10  C(r0/Ru)(r/X)l/2]  .  DS  (5) 

where  Ry  is  the  direct  range  from  source  to  receiver. 

The  equivalence  of  DS  in  the  laboratory,  where  dimensions  are  measured  in 
centimeters  and  frequencies  go  to  100kHz,  and  in  the  ocean  with  ranges  in  Km 
and  frequencies  in  hundreds  of  Hertz,  is  shown  in  Fig.  6. 

The  asymptotic  approach  to  the  far  field  diffraction  strength  is  presented 
as  a  function  of  r/  \  in  Fig.  7.  It  will  be  noted  that  the  plane  14°  wedge, 
the  simplest  approximation  to  the  seamount,  has  a  DS  some  lOdB  too  high;  and 
the  “contour  wedge"  which  is  contoured  to  the  seamount  profile  along  the 
sound  track,  and  which  is  the  physical  realization  of  the  two  dimensional 
parabolic  equation  solution,  is  also  several  dB  too  high. 


METHOD 


The  purpose  of  this  paper  is  to  demonstrate  the  extent  to  which  the 
diffraction  strength  of  the  seamount  may  be  determined  by  simple  computer 
models  based  on  plane  wedge  diffraction.  The  degree  of  validity  of  each 
computer  model  Is  measured  by  its  ability  to  match  the  results  of  the 
laboratory  scale  model,  whose  effectiveness  has  already  been  proven  by 
comparison  with  ocean  experiments.  Three  types  of  computer  models  are 
considered:  {1}  a  double  wedge  interpretation  of  the  two-dimensional  contour 
of  a  seamount;  (2)  a  segmented  plane  wedge  idealization  of  the  crest  of  the 
three-dimensional  seamount;  (3)  a  multiple  wedge  reepresentation  of  the 
three-dimensional  seamount. 

Model  1 : 

Impulse  calculations  of  double  diffraction  have  been  used  to  describe 
the  noise  shadowing  by  a  thick  barrier.  An  example  of  the  success  of  this 
technique  for  a  thick  plate  is  shown  in  Fig.  8.  The  two-dimensional  double 
diffraction  technique  is  now  used  to  approximate  the  seamount  as  a  wide 
barrier  with  two  changes  of  slope;  diffraction  at  the  first  edge  Is 
Interpreted  as  due  to  a  line  of  properly-timed  secondary  sources  which 
reradiate  over  the  second  edge.  Fig.  9  presents  the  two-dimensional  Track  6 
ray  path  contour  of  DICKENS  seamount  and  two  different  approximations  for 
the  double  diffraction  calculations.  In  approximation  one,  plane  surfaces 
are  assumed  from  the  source  to  position  a  where  the  sound  diffracts  through 
angle  U#  ■  186* *  from  a  to  position  c  where  the  diffraction  angle  Is 
D#  *  210°,  and  from  c  to  the  receiver.  In  approximation  two,  the 
assumed  plane  surfaces  are  from  source  to  b,  b  to  d,  and  d  to  receiver;  the 
two  wedges  thereby  defined  have  0*  *  212°  and  186°  at  b  and  d, 
respectively.  The  comparison  of  the  double  diffraction  computer  calculation 
with  laboratory  experiment  shows  good  to  excellent  agreement. 
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Model  2: 

A  nai've  assumption  will  now  be  evaluated,  namely  that  the  crest  of  the 
three-dimensional  seamount  is  where  most  of  the  major  changes  of  slope 
occur,  and  that  this  may  be  approximated  by  crestal  line  segments  of  finite 
length,  Fig.  10.  Each  segment  produces  a  response  that  is  a  truncation^ 
of  the  temporal  response  given  by  eq.  (1),  due  to  the  limited  time  of 
interaction  with  the  finite  segment.  The  diffraction  at  the  crest  is  then 
comprised  of  the  sum  of  the  plane  wedge  impulse  responses.  When  this  sum  of 
the  segmental  responses  for  DICKENS  Seamount  is  transformed  to  the  frequency 
domain  one  obtains  the  DS  shown  in  Fig.  10.  Only  the  lowest  frequencies  are 
seen  to  be  in  agreement  with  experiment,  as  one  should  expect  for  this  plane 
wedge  model  which  ignores  changes  of  slope  along  the  slopes  of  the  seamount. 
The  DS  of  this  crestal  segment  model  is  1  or  2  dB  better  than  that  of  the 
simple  plane  wedge,  0#  *  208°,  but  it  is  still  about  10  dB  too  large, 
compared  with  laboratory  measurements. 

Model  3: 

The  correct  calculation  should  consider  the  multiple  diffraction  at  all 
significant  changes  of  slope  at,  and  adjacent  to,  the  ray  path.  The 
simplifying  assumption  is  now  made  that  we  can  visually  identify  the  major 
changes  of  slope  near  TRACK  6  and  use  these  for  an  adequate  calculation  of 
DS.  Several  wedge  identifications  have  been  made  at  regions  of  large  change 
of  slope.  Fig.  11  shows  these  wedge  approximations  drawn  on  a  photograph  of 
the  critical  part  of  the  seamount  model.  We  find  (Fig.  12)  that  computer 
calculations  for  four  such  approximations  are  within  2,  3  dB  of  each  other 
and  about  2  or  3  dB  too  high  compared  to  the  scale  model  experiment.  The 
Identification  of  larger  numbers  of  changes  of  slope  would  be  needed  to 
improve  the  accuracy  of  the  shadowing  prediction. 


CONCLUSIONS 

Using  topographical  data  that  show  the  lines  of  significant  changes  of  slope, 
the  double  diffraction  technique  permits  accurate  computer  calculations  of 
the  diffraction  loss  of  poise  in  the  shadow  of  a  two-dimensional  seamount, 
and  approximate  predictions  for  the  three-dimensional  body. 
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FIG.  1  LABORATORY  MODEL  OF  DICKENS  SEAMOUNT  TO  SCALE  I  km  -  5  Inches 


FIG.  2  LABORATORY  MODEL  OF  DICKENS  SEAMOUNT  WITH  SCALED  MODEL  OCEAN  REPRE¬ 
SENTING  3$  kn  MIND  CALCULATED  FROM  PIRRSON-MOSKOWIT?.  MIND  NAVE 
SPECTRUM 
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OCEAN  RAY  DIAGRAM 
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LABORATORY  MODEL  OF  PROPAGATION  OVER  DICKENS  SEAMOUNT 


FIG.  3  RAY  DIAGRAM  OF  FORWARD  SCATTER  AND  DIFFRACTION  BY  DICKENS  SEAMOUNT 
ALONG  TRACK  6  OF  THE  OCEAN  EXPERIMENT 
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FIG.  4  PREDICTED  SHADOWING  LOSS  CAUSED  BY  DICKENS  SEAMOUNT  FOR  SOURCE-RECEIVER 
CASE  SHOWN  IN  FIG.  3.  Shaded  area  shews  error  bounds  of  scale  model 
calculation.  Circles  are  ocean  data  provided  by  N.R.  chapman  and 
G.R.  Bbbeson,  Defense  Research  Establishment ,  Pacific,  Victoria  B.C. , 
Canada.  The  theoretical  slope,  solid  line,  is  3  dB/ octavo 
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PRESSURE  DIFFRACTED  AND  MULTIPLY  REFLECTED  MIRROR  OCEAN  OVER  SEAMOUNT  MODEL 


PRESSURE  DIFFRACTED  AND  MULTIPLY  SCATTERED  TYPICAL  ROUGH  OCEAN  OVER  SEAMOUNT  MODEL 
♦°”i  I  I  1  I 


FIG.  5  LABORATORY  STUDY  OF  IMPULSE  DIFFRACTION  AND  MULTIPLE  REFLECTION  AT 
DICKENS  SEAMOUNT.  Top,  free  field  pressure  of  source ;  second  graph, 
pressure  diffracted  over  seamount  model;  third  graph,  pressure 
diffracted  and  multiply  reflected  by  a  mirror  smooth  ocean  surface ; 
bottom  graph,  pressure  diffracted  and  multiply  reflected  at  a  model 
sea  surface  for  a  35  kn  wind.  The  first  SO  microseconds  is  the 
diffracted  signal;  from  50  to  203  usee  is  scale  model  reflection 
scatter  tnaL  does  not  occur  at  sea;  from  203  usee  the  signal  is 
multiple  scatter  that  represents  the  ocean  situation 
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FIG.  6  DIFFRACTION  STRENGTH,  DS,  FOR  A  SYMMETRICAL  PLANE  RIGID  WEDGE  WITH 

0^  *  208*  and  l  *  0*  The  solid  line  gives  the  theoretical  value  for 

r  *  rQ  »  3.5  ku  in  water;  data  points  are  for  a  laboratory  model 

experiment  in  air  at  three  ranges,  25,  35  and  45  cm 
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FIG.  7  DIFFRACTION  STRENGTHS  FOR  THREE  MODELS  OF  DICKENS  SEAMOUNT.  Solid 
lines  from  top  to  bottom  are:  theory  for  plane  wedge  of  0^  =  208° ; 

asymptote  for  the  two-dimensional  contour  wedge  along  the  sound 
track ;  asymptote  for  the  three-dimensional  model  along  the  sound 
track .  Experimental  data  given  by  circles,  triangles  and  crosses/ 
respectively. 
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FIG.  8  SHADOWING  LOSSES  BY  DOUBLE  DIFFRACTION  CALCULATIONS ,  SOLID  LINES, 

COMPARED  WITH  EXPERIMENTAL  MEASUREMENTS  FOR  THICK  PLATES . 

The  range*  are  r  •  r  •  25  cm 
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TWO  DOUBLE  WEDGE  APPROXIMATIONS 
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FIG.  9  TWO  DOUBLE  DIFFRACTION  COMPUTER  CALCULATIONS  OF  DIFFRACTION  STRENGTH 

COMPARED  WITH  LABORATORY  EXPERIMENTAL  MEASUREMENTS  OF  THE  TWO- 

DIMENSIONAL  "CONTOUR  MODEL"  OF  SOUND  TRACK  6  OVER  DICKENS  SEAMOUNT. 

The  double  diffraction  barriers  are  plane  surface  fittings  to  the 

true  contour  (inset) .  Single  plane  wedge  theory  for  0  »  208° 

(14°  slopes)  and  0  «  212°  (16°  slopes)  are  shown  forw comparison 

ww 
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FIG.  10  PREDICTED  DS  FOR  CRESTLINE  SEGMENTAL  WEDGES  AS  SHOWN  IN  INSET.  f 

Laboratory  measurements  for  DS  for  scale  model  of  DICKENS  seamount  | 

shown  by  crosses.  The  dotted  line  is  the  solution  for  a  single  ] 

plane  wedge  Qy  *  206*. 
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FIG.  11  PHOTOGRAPH  OF  THE  DICKENS  SEAMOUNT  MODEL  WITH  IDENTIFICATION  OF 
WEDGE  CRESTS  USED  FOR  THE  DOUBLE  DIFFRACTION  CALCULATIONS  NEAR 
SOUND  TRACK  6 
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FIG.  12  MULTIPLE  DIFFRACTION  CALCULATIONS  OF  DS  FOR  WEDGES  THAT  ADJOIN 

THE  PRINCIPAL  RAY  PATH  CONTOUR  (TRACK  $)  COMPARED  WITH  LABORATORY 
MEASUREMENTS  FROM  THE  THREE  DIMENSIONAL  SCALE  MODEL  OF  DICKENS 
SEAMOUNT.  Referring  to  Fig .  11,  approximation  1  consists  of 
segments  A 1  and  B1  and  a  wedge  between  A2  and  A3  and  between  82 
and  83.  Approximation  2  is  comprised  of  the  double  wedge  Al,  A2 
and  double  wedge  Bl ,  82.  Approximation  3  uses  a  wedge  between 
A2,  A3  and  82,  B3  and  a  second  (unmarked)  wedge  closer  to  the 
receiver.  Approximation  4  is  a  triple  diffraction  which  is  composed 
of  wedges  Al,  81  and  A2,  82  and  A3,  83 
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